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Abstract 

Gauge theories of the Yang-Mills type are the single most important building block of 
the standard model of particle physics and beyond. They are an integral part of the strong 
and weak interactions, and in their Abelian version of electromagnetism. Since Yang-Mills 
theories are gauge theories their elementary particles, the gauge bosons, cannot be described 
without fixing a gauge. Therefore, to obtain their properties a quantized and gauge-fixed 
setting is necessary. 

Beyond perturbation theory, gauge-fixing in non-Abelian gauge theories is obstructed 
by the Gribov-Singer ambiguity, which requires the introduction of non-local constraints. 
The construction and implementation of a method-independent gauge-fixing prescription 
to resolve this ambiguity is the single most important first step to describe gauge bosons 
beyond perturbation theory. Proposals for such a procedure, generalizing the perturbative 
Landau gauge, are described here. Their implementation are discussed for two example 
methods, lattice gauge theory and the quantum equations of motion. 

After gauge-fixing, it is possible to study gauge bosons in detail. The most direct access is 
provided by their correlation functions. The corresponding two- and three-point correlation 
functions are presented at all energy scales. These give access to the properties of the gauge 
bosons, like their absence from the asymptotic physical state space, particle-like properties 
at high energies, and the running coupling. Furthermore, auxiliary degrees of freedom are 
introduced during gauge-fixing, and their properties are discussed as well. These results are 
presented for two, three, and four dimensions, and for various gauge algebras. 

Finally, the modifications of the properties of gauge bosons at finite temperature are 
presented. Evidence is provided that these reflect the phase structure of Yang- Mills theory. 
However, it is found that the phase transition is not deconfining the gauge bosons, although 
the bulk thermodynamical behavior is of a Stefan-Boltzmann type. The resolution of this 
apparent contradiction is also presented. In addition, this resolution provides an explicit 
and constructive solution to the Linde problem. 

Thus, the technical and conceptual framework presented here can be taken as a basis how 
to determine correlation functions in Yang-Mills theory, therefore opening up the avenue to 
investigate theories of direct practical relevance. The status of this effort will be briefly 
described, alongside with connections to other approaches to Yang-Mills theory beyond 
perturbation theory. 



Preprint submitted to Physics Reports 



December 26, 2011 



Keywords: Yang-Mills theory, Gauge-fixing, Dyson- Schwinger equations. Lattice gauge 
theory. Correlation functions. Landau gauge. Propagators, Vertices, Temperature, Phase 
diagram 



2 



Contents 



1 Introduction 

2 Quantizing Yang-Mills theory in Landau gauge 8 

2.1 The structure of gauge orbits 

2.2 Perturbative gauge-fixing 

2.3 Correlation functions 

2.3.1 Propagators and vertices 

2.3.2 The construction of gauge-invariant observables 

2.4 Non-perturbative gauge- fixing and Gribov copies 

2.5 Proposals for resolving the Gribov-Singer ambiguity 

2.5.1 Gribov regions 

2.5.2 Minimal Landau gauge 

2.5.3 Absolute Landau gauge 

2.5.4 Landau- -B gauges 

2.6 Gauge relations 

3 Methods 

3.1 Lattice gauge theory 

3.1.1 The action 

3.1.2 Fixing the gauge 

3.1.3 Correlation functions 

3.1.4 Lattice artifacts 

3.2 Quantum equations of motion 

3.2.1 Formulating the equations 

3.2.2 Renormalization group methods 

3.2.3 Uhraviolet 

3.2.4 Infrared 

3.2.5 Intermediate momenta 

3.2.6 Numerical methods 

3.2.7 Dyson-Schwinger equations on a lattice 

3.3 Slavnov- Taylor identities 

3.4 Renormalization 

4 Describing gluons at zero temperature |^ 

4.1 Correlation functions ul 

4.1.1 Propagators 

4.1.2 Vertices 

4.2 Schwinger functions, mass, and confinement 

4.3 Confinement scenarios 

4.3.1 Perturbative BRST [100 

4.3.2 Non-perturbative BRST and the Kugo-Ojima confinement scenario . Il02 



3 



4.3.3 Gribov-Zwanziger confinement scenario . 

4.3.4 Gribov-StingI scenario 

4.4 Relations to topological structures 

4.5 Beyond Yang-Mills theory 

4.5.1 Matter fields 

4.5.2 Composite objects 

Describing gluons at finite temperature 

5.1 Finite temperature 

5.1.1 Formulation 

5.1.2 Modifications of the methods 

5.2 Propagators 

5.2.1 Infrared 

5.2.2 Intermediate momenta and temperatures 

5.2.3 Infinite-temperature limit 

5.2.4 Schwinger functions and confinement . . 

5.3 Thermodynamic quantities 

5.3.1 Thermodynamic potential 

5.3.2 Order parameters 

5.4 Vertices 

Summary and discussion 

6.1 The state of the art, pitfalls, and odds and ends 

6.2 Summary and outlook 



4 



1. Introduction 



The best theoretical description of particle physics is currently arguably the standard- 
model of particle physics [l|, [ij , despite indirect as well as conceptual evidence for its incom- 
pleteness The most significant of the latter is the missing gravitational interactions. 
Otherwise, the standard model contains all known interactions: The strong and weak nu- 
clear force, electromagnetism, as well as the, at time of this writing still hypothetical 0, [sj 
Higgs interaction. 

The basic building block of the standard model, as well as of a multitude of extensions 
of the standard model are gauge theories These consists of a gauge sector, 

described by an Abelian or non-Abelian Yang-Mills theory, and a number of matter fields 
coupled to these gauge fields. The only exception to this rule is the Higgs interaction, which 
takes the form of a Yukawa theory . 

Thus, the entire standard model of particle physics rests on the description of gauge 
fields using Yang-Mills theories. It is this central element of particle physics with which 
the following will be concerned, dropping almost always all matter fields. The aim is to 
describe the elementary degrees of freedom of Yang-Mills theory, the gauge bosons, as well 
as their interaction. In this context the fact that Yang-Mills theory is a gauge theory will 
play a central role, as that makes the description of the elementary excitations necessarily 
gauge-dependent. This subtlety will be discussed in detail in the following. 

The starting point of any such discussion is the classical Lagrangian of Yang- Mills theory 

c = -\f^f!^ (1) 

i^;. = d^Al-d^A'^ + gf^^A'^At. (2) 

Herein the fields are the gauge fields describing the gauge bosons, called here for simplicity 
gluons. Of course, these could easily be selected to be the weak isospin bosons or the 
hypercharge gauge bosons. Furthermore, the parameters of the Lagrangian are the coupling 
constant g, and the numbers Z"^'^, which are the structure constants of the associated gauge 
algebra A. This gauge algebra can, in principle, be any semi-simple Lie algebra. In case of 
the standard model it is A =su(3) xsu(2) xu(l) jsf. The factor su(3) generates the strong 
interactions, and the factor su(2)xu(l) yields, after mixing, the weak and electromagnetic 
interactions. Herein, the discussion will be restricted to non-Abelian gauge algebras, as 
without matter fields Abelian gauge theories are trivial theories of non-interacting gauge 
bosons, as they miss any kind of interactions between the gauge bosons jij. 

For a semi-simple algebra the interactions of Yang- Mills theory factorize into the dynam- 
ics of Yang- Mills theories with independent simple algebras. They would only be linked by 
the presence of matter fields, which in the standard model are the leptons, the quarks, and, 
possibly, the Higgs field. In the standard model, these matter fields couple only minimally 
to the gauge fields, leading in case of the electroweak sector also to a hiding of the symmetry 
by the Higgs mechanism [2|. It thus suffices to consider here simple gauge algebras. 
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However, even without the presence of the matter fields, non-Abehan Yang-Mills theories 
are far from trivial. As long as the gauge algebra is non-Abelian the structure constants 
jabc g^j^g non- vanishing. Then, due to the quadratic part of the field-strength tensor ([2]), 
interactions between the gauge fields are introduced. These interactions are found to obey 
asymptotic freedom, i. e., they become weaker with increasing energy On the other 
hand, with decreasing energy they become stronger. As a consequence, perturbation theory 
can only capture the leading behavior at large energies, but fails completely when the energy 
of a process reaches the typical scale of the theory, which is denoted as Aym in the following. 
E. g., in the strong interactions, this scale is of the order of 1 GeV. With reaching this scale, 
all types of genuine non-perturbative effects are present, in particular the existence of bound 
states like glueballs, or the confinement of gluons joj. 

Since these effects are qualitatively present even in the absence of matter fields, this is the 
ratio to first investigate here the simpler case of neglecting matter fields, and concentrate on 
Yang-Mills theories with a simple Lie algebra. Understanding this system non-perturbatively 
then provides a firm basis for going back to the standard model. It is one possible version 
of this necessary foundation of understanding gauge theories with which the following will 
be concerned. 

At the heart of many obstacles to be encountered in this process is the fact that in 
Yang-Mills theory redundant degrees of freedom are introduced for the sake of having a 



local quantum field theory [10| . As a consequence, Yang- Mills theory is a gauge theory. One 
of the most striking features of the Lagrangian ([T]) is therefore its invariance under local 
gauge transformations, which take the infinitesimal form 

a; ^ Al + Df<l>, (3) 
= 5'^'d, + gffAl 

where 0" are some arbitrary functions. Therefore, the gauge fields themselves cannot be 
entities of the physical reality, as any observations should be independent of the chosen 
gaug^. Consequently, any particle-like excitations associated with the gauge fields, also 
cannot represent physical observable particles. In fact, the particle-field duality 13] turns 



out to be only a high-energy property [2|, and at low energies the fields lack a structure 
which can be easily interpreted as a particle-like excitation. The theory must therefore in 
some way remove them from the physical spectrum, an effect which cannot be captured 
by perturbation theory j^: They are said to be confined. Indeed, there appears to also 
exist no mechanism, which could turn the gauge-dependent gluons into gauge-independent 



ones 



10|, since otherwise it should be possible to experimentally detect isolated gluon-like 
states. Within experimental uncertainties, this is not the case [jj. These properties will be 
discussed extensively in section 14.21 

Since the confinement of gluons is not described by perturbation theory, which has asymp- 
totically observable gluon states it is necessarily non-perturbative. Furthermore, the 



^For the theory to be consistent it is required that no anomahes are present fill]. Whether this is the 
case for Yang-Mills theories beyond perturbation theory has not yet been proven, but no evidence exists 
which suggests the presence of anomalies, and their absence will therefore be assumed here. 
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behind-the- moon-problem illustrates that the mechanism is also nonlocal [12|. This con- 
finement of gluons is not yet fully understood. Here, possible scenarios will be described in 
section I4.3[ though important connections still remain to be made, see in particular section 
14.41 In addition, the underlying mechanism could very well be gauge-dependent, but even 
an understanding in one gauge would be substantial progress. However, the gauge-invariant 
statement of gluon confinement must be valid in any gauge. 

With this in mind, in the following the discussion will be essentially restricted to the 
Landau gauge, which will be introduced belovjl. However, the statement of what Landau 
gauge means is non-trivial beyond perturbation theory, as a consequence of the Gribov- 



Singer ambiguity [iTj, |l8|. This ambiguity and its consequence will be discussed in great 
detail in sections 12.41 and 12.51 One possibility is to split the perturbatively unique Landau 
gauge in a family of gauges beyond perturbation theory. This will be discussed in more 
detail in section 12.51 Occasionally, comments on the state of affairs in other gauges will also 
be made, but this is not the main focus here. Of course, eventually the aim can only be to 
understand the same mechanism in at least a sizable number of different gauges, similar to 
the concept of gauge-parameter independence of perturbation theory [2]. 

Irrespective of the particular choice, once a gauge is unambiguously fixed, the concept of 
gluons is well-defined, and questions relating to their properties can be posed. This will be 
done in section |H The methods used for this purpose are described beforehand in section 
[31 and are the equation of motions and lattice gauge theory. These will be used jointly 
to determine the correlation functions, a concept introduced in section 12.3. of gluons. 
From the two-point functions, the propagators, it is possible to infer properties of gluons. 
By determining the corresponding three-point functions it is possible to determine their 
relative interaction strength, the running gauge coupling, but also estimates of decay rates 
and fusion processes. Ultimately, this will give access to much more complicated quantities 
like scattering amphtudes in the future. 

These investigations can be extended to finite temperature, which will be done in chapter 
|5l This setting is of quite some importance when conside ring th e early universe, but also 



for its laboratory-based recreation in heavy-ion collisions |19l-l22|. Here, the change of the 
properties of gluons will be read off from the temperature-dependent propagators, which 
will be determined once more by the quantum equations of motion and lattice gauge theory 
simultaneously. It is found that gluons are not deconfined at any temperature, though they 
generate a thermodynamic potential which at high temperature coincides with that of a 
gas of non-interacting gluons. This apparent contradiction will be resolved, showing that 
long-range physics is sub-leading to hard effects for bulk thermodynamics. As a by-product. 



this will yield a constructive resolution of the Linde problem [23 . 

Finally, the state of the art, the required steps to make this a full formal procedure, and 
possible pitfalls will be summarized in section [61 Also the results will be set into perspective 
to other investigations and other theories, in particular QCD, over the course of this text, 
in particular in section 14.51 The final result of the presented set of methods and concepts is 



^Reviews, which also cover some parts of the foUowing, though with partially different eniphasis, are 



a framework for the description of gauge theories beyond perturbation theory at all energy 
scales. Though by far not yet a simple out-of-the-box solution given the complexity of the 
subject, it represents a versatile toolbox to describe the gauge- invariant, observable physics 
of gauge theories using explicitly the elementary degrees of freedom, the gauge bosons. It 
thus makes explicit contact between the basic structure of the theory and the consequences 
of it, which can be measured in experiment, irrespective of the involved scale, and without 
the addition of further parameters to the original theory, as is necessary in model or effective 
theory abstractions. The continued development of this framework and its application to 
various problems is therefore a very active field of research. Especially, it has already been 
applied to a wide range of theories, from Yang-Mills theory itself, as described here, over 
QCD [13|, llJ], Higgs ph ysic s |24| . Technicolor extensions of the standard model 
to quantum gravity 



27-21 



The starting point of all of these investigations is, of course, quantized Yang-Mills theory. 



2. Quantizing Yang-Mills theory in Landau gauge 

2.1. The structure of gauge orbits 

As noted, the Lagrangian ([1]) is invariant under the infinitesimal gauge transformation 
([3]). In fact, it is also invariant under the finite gauge transformation 

Afj, = TaA^ 

h = T,r 

where t"' are the generators of the gauge algebra in the fundamental representation. As a 
consequence, the set of fields connected by gauge transformations 

= {aP VM, (4) 

are all equivalent representations of a given, fixed space-time history A^ of the gluon field. 
Such a space-time history will be called a configuration in the following. The set Q depending 
on such a configuration is called its associated gauge orbit. The value of all objects which do 
not change under a gauge-transformation and thus are equal for all members of the gauge 
orbit, e. g. the action, are called gauge- invariant. On the other hand, objects which are 
changing under gauge-transformations will be called gauge- dependent. In general, they will 
depend on which element of the set Q is selected to calculate them. The chosen element of 
Q is called the representative of the gauge orbit. However, this is not necessary. Gauge- 
dependent quantities may still be invariant under the choice of elements from a subset of Q, 
or may change only between different subsets of Q. 

To determine the expectation value < Q > of a quantity Q{A^), depending on the field 
variables A^, the path-integral formalism can be useclf] j3j|. The value of < Q > is then 



■^From here on everything wih be given in Euchdean space-time in natural units. The reconstruction 

8 



given by 

<Q>= I VA^Q{A,)e-^''"^'', (5) 



where the measure "PA^ is normaUzed such that < 1 >= 1. The functional integral is over 
all gauge orbits of all configurations. It is then necessary to ensure that the integral is only 
sampling each configuration with the same weight using adequate normalization to obtain 
a well-defined result for Q. For gauge-invariant quantities, this is sufficient jo], [Hj . For a 
gauge-dependent quantity, like the gauge- fields, this is not sufficient. The integration over 



all representatives of a gauge-orbit yields that all such expectation values vanish [3J|. To 
obtain a non-zero value, it is necessary to define how to select a representative, or a weighted 
subset of representatives, of each gauge orbit, which is equivalent to choosing a coordinate 
system in terms of the field variables |2j. This has to be done by introducing a gauge-fixing 
function into the path integral weight, which cuts the domain of integration. Only when 
performing the integral over such a reduced gauge orbit, non-zero values for gauge- dependent 
observables can be obtained. 

How to perform this is not only interesting as a conceptual question in its own right, nor 
just because it is not possible to discuss gluon properties without fixing a gauge, but also for 
practical reasons. As will be discussed in section I2.3.2[ it is very convenient to use gauge- 
dependent quantities in intermediate stages to finally obtain gauge- invariant, observable 
quantities. 

2.2. Perturbative gauge-fixing 

The simplest way to implement a selection criterion is obtained in perturbation theory. 
Furthermore, the implementation in perturbation theory is also a viable starting point to 
generalize to the full theory. 

In perturbation theory, it is sufficient to implement a local condition to select a unique 
representative on the gauge orbit ji], [sij. One possibility to do so is by requiring 

d,A; = 0, (6) 

the Landau gauge condition. This gauge condition will be used exclusively here, i. e., 
whatever other conditions a field configuration satisfies in the following, it will always satisfy 
OH]) as well. 

The Landau gauge is a limiting case of the class of covariant gauges (if, which are well 
suited for perturbative purposes. Many other gauge conditions, like Coulomb gauge, max- 
imal Abelian gauges, interpolating gauges, axial gauges, Curci-Ferrari gauges, background- 
gauges, and others have also been used. However, Landau gauge has a number of distinct 
advantages. One is that it is manifestly covariant, and none of the technical complications 



of any Minkowski quantity can be done either by explicit Wick rotation [30] or by usage of the Schwinger 
function reconstruction [lOj. In a lattice formulation it can be proven that this is possible for all gauge- 
invariant quantities [3l]|. A proof for gauge- dependent quantities, necessary to make the theory well-defined 
in Euclidean space-time, is still lacking, but henceforth it will be assumed to be possible. See [sl] for more 
details. 
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associated with non-manifest covariance appear [36|. Furthermore, renormahzation is most 
simple in Landau gauge, as the degree of the divergence is the lowest of all possible, i. 
e., just logarithmically divergent and not quadratically divergent in four dimensions. As a 
consequence, all correlation functions in lower than four dimensions are finite. Finally, in 
Landau gauge, transverse tensor structures and longitudinal tensor structures of correlation 



functions, with respect to gluon momenta, are as decoupled as possible [37[. In particular, 
the minimum number of tensor structures is required when treating non-amputated cor- 
relation functions. All these properties make Landau gauge technically convenient. As a 
consequence, it is by all means the best-studied gauge beyond perturbation theory so far. 

The condition is sufficient to perturbatively single out one representative of each gauge 
orbit, up to global gauge transformations. This follows from the fact that no transformation 
functions 0° exist, which takes a gauge copy satisfying the Landau gauge by an infinitesimal 
transformation ([3]) into a gauge copy which also satisfies the Landau gauge condition ([6]). 
When expanding a finite gauge transformation in powers of the coupling constant, this also 
holds. 

The Landau gauge condition (E]) can be introduced as a restriction of the space of all 
gauge orbits to a hypersurface of representatives satisfying the Landau gauge condition 



by either use of the Faddeev- Popov procedure [2|, |38[, or by the more general anti-field 



formalism [11|, |35|. Both equivalent formalisms introduce two additional auxiliary anti- 
commuting, scalar fields c"" and c", the ghost and anti-ghost fields belonging to the adjoint 
representation of the gauge algebra. This permits to construct a local formulation of the 
restriction of the gauge orbit to a hypersurface. As a result, the path integral ([5]) takes the 
form 

<Q> = \im J VA^VcVcQ{A^,c,c)e-^'^'''^'> (7) 

in d dimensions. The introduction of a gauge parameter ^ is necessary in an intermediate 
step to make all possible inversions of operators well-defined. In implementations using 
lattice gauge theory, see section 13. it is actually not necessary to make explicit references 
to it, and also in calculations using functional methods, see section 13. 2[ it is possible to 
remove this parameter rather early on. It will therefore be mostly left implicit hereafter. It 
should be noted that for all such expressions implicitly necessary weight-factors have been 
included in the measure to ensure that still < 1 >= 1 is valid. 

Of course, given the presence of auxiliary fields it is possible to also construct quantities 
Q including them. From the appearance of the covariant derivative it follows that the ghosts 
interact in non-Abelian gauge theories with the gauge fields. Thus it is mandatory to include 
them in any calculations, which makes explicit use of the gauge-fixed path integral in the 
form ([7]). However, they do not represent physical particles, and as such will not appear 
in observable physical states, but only in intermediate states. Nonetheless, it is possible 



is_p 

to form bound states from them, called ghost balls, or as gluon-ghost hybrids 12]. Such 
bound-states are still unphysical, and thus still not be observable, even if they are color 
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singlets. 

After having obtained a perturbatively gauge-fixed setting, it is possible to determine 



correlation functions in perturbation theory, using standard methods [33[. For going beyond 
perturbation theory, it is quite useful to first introduce the concept of correlation functions 
in more detail, as they will play a role when attempting to gauge-fix beyond perturbation 
theory in the way proposed here. 

2. 3. Correlation functions 
2.3.1. Propagators and vertices 

The basic entities, which determine the complete partition function, are the Green's or 
correlation functions [33|- In their simplest version, these are polynomials of n fields 

< Aa,{Xi)...AaSXn) >= [ PA, (o^i ) . . . C" ^ '^'"^ (8) 



with Oj a multi-index encompassing the type of the field as well as color and Lorentz indices. 
The number n of fields involved will henceforth be called the order of the correlation function. 
From these the connected correlation functions can be derived, which are defined for the 
two-point functions as 

< A1A2 >c=< A1A2 > - < Ai>< A2 >, (9) 
and correspondingly for higher n-point functions Furthermore, the amputated correla- 



tion functions, or vertex functions, F* fulfill the relation |33 

n <'°' 

In particular, for the two-point functions, the propagatorqfl, the relation 

< Aaix)At,iy) >,= T''^^^-\x -y)= D^^^^ix - y), (11) 

holds. It has been assumed that the system is translationally invariant, which will be the 
situation encountered throughout. It should be noted that for Yang-Mills theory the one- 
point functions in all covariant gauges, and in particular in Landau gauge, always vanish, 
as otherwise Lorentz symmetry would be broken. Therefore, connected and full Green's 
functions coincide for two- and three-point correlation functions, and a difference is only 
encountered for higher correlation functions. 

The relation (El) p ermits the definition of the generating functional for the various cor- 
relation functions [l3|. The equations for the full correlation functions are obtained from 



6 



6 



+ ja]Z[j] = (12) 



SA, 

Z[j] = I I)A,e-^+^'^'"^'^^», (13) 



"^The inverse appears essential due to the amputation of the two-point correlation functions. 
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where the ja are sources to be set to zero at the end of the derivation. The derivative 
in the argument of the expression of the action imphes that all appearances of the fields 
have to be replaced by a derivative with respect to the sources, inheriting their indices. 
This equation yields the one-point correlation functions. Higher correlation functions are 
obtained by taking further functional derivatives of this equation with respect to the fields. 
The connected correlation functions are obtained from 



6S 



6A, 
W = \nZ 



5W_ 6_ 



+ Ja = 



(14) 



where W is also called the free energy. The vertex functions are determined from 



6T 6S 



5A^ 
-W 



Aa + 



(15) 



Ja- 



Thus, the Legendre transform F of W, the effective action, is the generating functional of 
vertex functions. In particular jij 



r^-i-^-n(^Xi,...,Xn) 

T^^^'{x,y) 



6Aai{xi)...6Aa„{x, 

D^-^^-\x-y) = 



(16) 



5^T 



6Aa{x)6Ai,{y) \Sja{x)6jb{y) 



6^W 



-1 



The order of the field indices Ai is relevant not only because of assignment of the arguments, 
but also if anti-commuting fields appear. 

From f ll6p follows the reconstruction of the original path-integral ( I13p as Q 

CO „ 

^[•^1 ^ S / d'^^l'-d'^^n < Aa^{xi)...AaSXn) > jai{Xl)---ja„{^n)- (17) 
n=0 

Similar reconstructions can be performed for W and F in terms of the connected correlation 
functions and the vertex functions, respectively. The equations for the vertex functions 
generated by derivations of (fT5|) with respect to the fields are known as Dyson- Schwinger 
equations 



13l . |33| , and are described in more detail in section 13.21 



As a general convention, in the following the function, or functions, which modify a 
propagator or vertex away from its tree-level value will be called a dressing function. E. g., 
for a propagator F, its dressing function 7 is defined as 



7 
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where the inversion is used to connect to the hterature [l3]j, while is the tree-level 
propagator. In general, for both propagators and vertices, additional tensor structures both 
in Lorentz and color indices can appear beyond tree-level. In this case more than one dressing 
function will be needed. They will be defined as dimensionless factor functions, and the new 
tensor structures are taken to not contain any information beyond those needed to specify its 
transformation properties and engineering dimensions. E. g., if the gluon 2-point vertex were 
transverse at tree-level, but acquired a non-trivial longitudinal dressing beyond tree-level, 
this would read ^ ^ 



Of course, this is not the case in Landau gauge jl, |4o[. In particular, the transverse gluon 
propagator dressing function will be denoted by Z throug houfl while the ghost propagator 
dressing function will be denoted by —G. As a general shorthand, D^^, will be used for the 
gluon propagator T^^~^, D for Z{p)/p'^, and Do for the ghost propagator F^"^"^, and color 
indices are kept explicit if not a factor of 5"^ has been factored out, i. e., D"-^{p) = 6°'^D{p). 

2.3.2. The construction of gauge-invariant observables 

By expanding the equation f|T5|) and its functional derivatives in the coupling (or another 



parameter), a perturbative set of equations for the correlation functions is obtained [33 
This permits straightforward computations of the perturbative behavior of the correlation 
functions. This process can be simplified by the use of Feynman rules j2| . 

In particular, the matrix elements of two-body decays are described by the three-point 
correlation functions, while two-to-two scattering cross sections matrix elements can be 



derived from the four-point correlation functions [2, |30| • The latter are the lowest correlation 
functions which can be combined to yield a momentum-dependent gauge-invariant quantity 
in Yang-Mills theory. Their computation is therefore the minimal requirement to obtain 
gauge-invariant information in general. In perturbation theory, their calculation can often 



be reduced to simpler calculations by, e. g., Cutosky rules [30[. Any such gauge-invariant 
result will necessarily always have all color indices in some way contracted, as on any non- 
contracted index a gauge-transformation could act, and thus would modify the result. 

Note that gauge-invariant physics also influences gauge-dependent correlation functions, 
for any order. E. g., a gluon can emit a virtual glueball, and reabsorb it later. Turning the 
argument around, there must exist indirect ways to infer gauge-invariant information from 
the gauge- dependent correlation functions. Therefore, in principle, it should be possible to 
obtain this information also from lower-order correlation functions. An example for this will 
be the determination of the phase transition temperature and order parameters from the 
two-point correlation functions in section 15.21 



^Note that in the context of (functional) renormahzation-group equations often the inverse of 7 is called 
a dressing function [soj. 

®Note that this function corresponds to the trace in both color and Lorentz indices of the full gluon 
propagator, multiplied by p^, thus being (— p)^°(p)) / {{d — 1)Na) with the dimension d of space-time 

and Na the size of the adjoint representation of the gauge algebra. It is thus always positive semi-definite. 
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It should be noted that, though the final result is gauge- invariant, it is often possible to 
choose a gauge which is especially suited to calculate a particular gauge-invariant quantity. 
This makes gauge-fixed calculations, though requiring more degrees of freedom, for many 
purposes quite attractive. E. g., almost all of perturbation theory is formulated in this way 
[2, 3^, though it is not necessary when only calculating gauge-invariant quantities jo], |41, 42 



The calculation of the hadronic spectrum in lattice gauge theory utilizes also gauge-fixed 



intermediate stages for a more efficient treatment [431 . 

Thus, the complete stage for perturbative calculations is now set. Unfortunately, many 
interesting effects, like bound-states, phase transitions, hiding or breaking of symmetries, 
confinement, etc. cannot be obtained in perturbation theory |i2|], and a treatment beyond 
perturbation theory is necessary. This is the main topic of this work. 

2.4- Non-perturbative gauge-fixing and Gribov copies 

To calculate gauge- dependent quantities non-perturbatively, in particular correlation 
functions, it is necessary to obtain a gauge-fixing which is valid beyond perturbation theory. 
However, beyond perturbation theory, gauge conditions like the Landau-gauge condition (|6]) 
have no longer a unique solution for a given configuration. There are several explicit exam- 



ples illustrating this fact [17|, I44j-|46| . Such independent solutions are called Gribov copies, 
and the associated ambiguity of the gauge condition is termed the Gribov-Singer ambiguity 
171 18| . This is a property which pertains to any non-Abelian gauge theory of Yang- Mills 

typqil. 

Unfortunately, it turns out that it is impossible to construct any kind of purely local 



gauge condition to single out a unique representative for a configuration |18l) I. This is a 
rather general result, which is based on the structure of the gauge orbits in non-Abelian 
gauge theories. As a consequence, it is only possible to specify a unique representative of a 
gauge-orbit when this specification is done non-locally. 

2. 5. Proposals for resolving the Gribov-Singer ambiguity 
2.5.1. Gribov regions 

In principle, it would be possible to average, in a well-defined way, over the Gribov copies 
to obtain also non-perturbatively a well-defined gauge, similar to, e. g., covariant gauges in 



perturbation theory. However, such Hirschfeld gauges [55|, 156[ induce significant cancella- 
tions, and no practical implementation has been constructed so far, but only conceptually 
developed [57-60]. Thus, the alternative is to select by some prescription a single represen- 
tative or a smaller subset for each gauge orbit, which therefore satisfies furthei|j constraints 
in addition to (El). 



'^In curved space-times, the properties of Gribov copies can change substantiahy, see e. g. [47H49j. Fur- 
thermore, Gribov copies also appear in gravity (46l. Isoj . Both findings are beyond the scope of this review. 

^This is only proven for covariant gauges [l8|. However, it has not yet been possible to find any local 
gauge condition which resolves the problem, and at best it can be traded in for a different problem of similar 



complexity, as is done, e. g, in direct Laplacian gauges 5ll-l53l|. For recent approaches to circumvent this 



problem, see e. g. |54{ . 

^The global gauge freedom will be left unfixed. Fixing it would lead, e. g., to Polyakov gauges [6lj or 
aligned gauges |62| . 
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Thus, further (non-local) constraints are required. To implement the constraints, a 
sequence of further conditions can be applied. In all cases, the first step taken is always to 
reduce the perturbative gauge freedom by implementing a local gauge condition, here the 
Landau gauge IQ. 



The remaining set of Gribov copies is called the residual gauge orbit [63||. Since the 



condition is perturbatively unique, only finite gauge transformations connect two different 



elements of the residual gauge orbit [6J]. This is trivially so, since any infinitesimal gauge 
transformation will move along the gauge orbit automatically out of the gauge-fixing hyper- 
surface implemented by the perturbative gauge-fixing by construction. 

The first restriction of the residual gauge orbit taken here is to constrain it to the so- 
called first Gribov region 17|. This Gribov region is defined by the requirement that the 



Faddeev-Popov operator M 



M"^ = -d^Dlf, (20) 
all of its eigenvalues are zero or positive. This region 



ab 



is strictly positive semi-definite, i. e 
can be shown to be bounded and convex 



65 



and the Faddeev-Popov operator has zero 
the so-called Gribov horizon. It can be 



eigenvalues only on the boundary of this regio 
shown that all gauge orbits pass at least once through the first Gribov region ^] , and hence 
no physical information is lost by restricting to it. The boundedness is a remarkable fact, as 
it implies that when calculating physical observables no arbitrarily large field fluctuations 
have to be taken into account. It contains the origin of fleld-space, and thus perturbation 
theory, as well. This follows from the fact that in the vacuum case fl2Up reduces to the 



positive semi-deflnite Laplacian [4^. Thus by restricting to the flrst Gribov region, ordinary 
perturbation theory is always included. 

Besides this flrst Gribov region, the remainder of the residual gauge orbit is a set of 
further Gribov regions. These are separated by further concentric Gribov horizons, each 
having more and more negative eigenvalues. The number of negative eigenvalues increases 
by one by passing the boundaries of these regions, but stays constant inside [l7|, |45| . It is 



expected that every residual gauge orbit passes through every Gribov region, though there 
is not yet an explicit proof of this. 



This restriction can be implemented using a ^-function [17|, |67|] in the perturbative gauge- 
flxed path integral ([7]) 



{-d,Df) 



lim 

n^(A, 



VA^VcVcQe i-d^Df) (A 



c)e 



fd-^xCg 



(21) 



where A, is the ith eigenvalue of the Faddeev-Popov operator fl20l) . Thus, only if all eigen- 
values are positive the ^-function contributes. To make contact to the implementation of 
this restriction in lattice gauge theory below in section [3.1.2[ the deflnition ^(0) = 1 has to 
be made for the step function. 



Apart from trivial zero modes due to constant eigenmodes, which will always be implicitly factored out. 
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Unfortunately, a unique, method-independent prescription how to effectively implement 
this restriction to the first Gribov region explicitly has not yet been constructed. There are, 
however, a number of possibilities, which have been explored. 

A proposal for how to implement this restriction using additional ghost fields, and thus in 
a similar way as in perturbation theory, has been made by approximating the 6'-function by 
a 5-function with the argument that in a high-dimensional space only the boundary contains 



an appreciable part of the volume [67|]. This generates the so-called Zwanziger Lagrangian. 
However, due to subtleties related to the definition of the step-function it is not yet proven 
that this is a valid procedure, though it has many interesting properties, and has been 



investigated in great detail, see e. g. |64 l67H79|. Furthermore, no Gribov copy, or any gauge 



copy in general, is preferred compared to another [80] . It would thus be completely legitimate 
to always chose the innermost Gribov copy for each gauge orbit. If (almost) all gauge orbits 
have a representative away from the Gribov horizon, this would yield distinctively different 
results for gauge-dependent quantities, e. g. the expectation value of the lowest Faddeev- 
Popov eigenvalue. Thus, such a replacement is already implementing a certain selection 
of Gribov copies, and thus corresponds to an extended gauge-fixing procedure. This is 
completely correct, provided (almost) all gauge orbits have Gribov copies on the Gribov 
horizon. Though not proven, this appears very likely. 

Another proposal [ill how to enforce this restriction explicitly stems from empiric ob- 
servations, and insights gained from two-dimensional Yang-Mills theory in Coulomb gauge 



82| . These suggest that a possible method-independent characterization of the first Gribov 



region would be the requirement that the ghost propagator 

= (c"c^) = (M'^^-i), (22) 

is negative-semi-definite in position space. Though at least the former statement is clearly 
a necessary condition, it has not yet been possible to show that this is a sufficient condition. 
Therefore, this has the status of a proposal which requires further investigations. However, 
for the methods employed here there are some method-dependent means how to restrict to 
the first Gribov region. This is described in section |3l 

After restricting to the first Gribov region, the remainders of the residual gauge orbits 



still possess a large number of Gribov copies j66|, l83l-l85| . This set will also be denoted as 



the residual gauge orbit in the following, to avoid the term residual of the residual gauge 
orbit. In fact, in an infinite volume this number is likely infinite, and in a finite volume V it 



appears to be a rapidly rising function of V [85| , possibly even proportional to exp V [81 



Actually, counting Gribov copies is in practice a non-trivial problem [87|], since two 
Gribov copies are different if and only if they differ at least at one space-time point after 
factorizing all possible global gauge transformations and all space-time transformations. 
This implies that for the decision whether two representatives of a gauge orbit are identical 
or Gribov copies, it is required to compare their field values at every space-time poiniF^. 



"'^^It appears that Gribov copies differ from each other over some large domain [88l|, so in practice aheady a 
coarse search can yield that two candidates are different. However, to ensure that they are the same requires 
a check of the whole space-time point by point. 
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Number of Gribov copies in four dimensions 




It is also in general non-trivial how to find all Gribov copies, so that they can be counted. 



Thus, except for very special systems j85j, only a lower limit can be posed on the number of 
Gribov copies. Examples for these are shown for two, three, and four space-time dimension 
on a finite lattic^ in figure [TJ It should be noted that, though two dimensions has trivial 
dynamics |91|] , gauge fixing has the same subtleties as in higher dimensions. Two-dimensional 
Yang-Mills theory is therefore an ideal laboratory to study these issues without the obscuring 
dynamics [o^. Remarkably, the number of Gribov copies not only increases with volume 
but in three and four dimensions also with increasing cutoff, while the latter seems not to 
be the case in two dimensions 



62, 87 



This is one of many hints encountered in this work 
that two-dimensional Yang-Mills theory has indeed some distinct features. 

Once more, it should be noted that one Gribov copy has no intrinsic difference compared 



^^Details of the lattice method are given in section [5TT| All units here and hereafter have been fixed by 
setting the string tension to (440 MeV)^, see for details |89ll90|. 
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to another Gribov copy, since they are physically equivalent. Thus any choice of a Gribov 
copy to represent the residual gauge orbit is equally acceptable. This is nicely illustrated by 
stochastic quantization, in which it is found that there is no stochastic force acting along 
a gauge orbit, and thus in the stochastic equilibration process no point on a gauge orbit 



is preferred [65|, |92|-|94|. This, of course, is just in disguise the problems encountered when 



defining the path integral, which require to introduce a gauge condition in the first place. 

2.5.2. Minimal Landau gauge 

As stated, the residual gauge orbits inside the first Gribov horizon possess further Gribov 
copies. It is therefore necessary to constraint the choice of Gribov copies further. There are 
two strategies mainly in use currently for that purpose. Both were originally motivated by 
studies on a finite lattice 



95, 96 



The first method is essentially a stochastic approach. In this case, instead of specifying 
conditions for selecting a Gribov copy, a random Gribov copy is chosen for each residual 
gauge orbit 0- This prescription, termed the minimal Landau gauge, therefore averages 
over Gribov-copy-dependent properties when calculating correlation functions. Assuming 
the choice to be ergodic, unbiased, and well-behaved, this implies that this prescription 
is equivalent to averaging over the residual gauge orbit 62, [87|. However, a constructive 



prescription how to make this choice in a path integral formulation is only developing [62|, |80 



87 



Precise definitions of this gauge therefore exist only as operational definitions in terms 

The second approach attempts to characterize 



of algorithms in lattice gauge theory [96 



Gribov copies and make a choice based on these characteristics. Two possibilities for this 
characterization will be presented in the next two sections. 

The central element of all operational definitions of the minimal Landau gauge is the 
fact that any Gribov copy in the first Gribov region maximizes the functional 8J, l97| 



F[A] 

{m) 



1 - 

1 - 



- j d'xAlAl 



(23) 



dpp'-'D;i{p) 



D 



ab 
fit/ 



{A^Al), 



on each configuration, where D'^^j^ is again the gluon propagator. This implies that this gauge 



minimizes the integrated weight of the gluon propagator. That this is indeed satisfying the 
minimal Landau gauge follows from the fact that the first derivative of fl23|l is the Landau 



gauge condition, and the Hessian is the Faddeev- Popov operator [97| . If any given algorithm 



finds one of all the maxima with equal probability, it would be a faithful representation of 
the distribution along the residual gauge orbit, and also be ergodic. 

This gauge has a further important potential. It has been argued that all Gribov copy 
selection procedures in the infinite- volume and continuum limit, the so-called thermody- 



namic limit, should yield the same result [65|], and thus therefore the one of the minimal 



Landau gauge. The basic idea behind this argument is essentially of thermodynamic nature: 
The configuration space is in this limit infinite-dimensional, thus almost all Gribov copies of 
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Raw distribution of F(A) in two dimensions 
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Figure 2: The distribution of ([23|) for a fixed dis- 
cretization of a = 0.22 fm, after subtracting one and 
rescaling by a factor of -1000 [6^. Results are shown 
for two (top- left panel), three (top-right panel), and 
four dimensions (bottom panel). Here and hereafter 
always 1000 configurations with 20 checked Gribov 
copies each have been used. 
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almost all residual gauge orbits should lie at the boundary of the first Gribov region. The 
important assumption is then that all Gribov copies at the boundary have similar correla- 
tion functions. Then the same argument as in thermodynamics holds, and the equilibrium 
behavior should emerge. Since the minimal Landau gauge is already averaging, this equilib- 
rium behavior should coincide with the one found in the minimal Landau gauge. However, 
this argument is only expected to hold for any finite polynomial in the gauge fields. E. g., 
the ghost propagator (12^ cannot be expressed as a finite series in polynomials of the fields, 
and therefore the argument does not need to hold for it. 

Take as an example the expression (123|) itself, which is a finite polynomial in the fields. 
Its distribution over the residual gauge orbit is shown in figure |2j It is clearly visible how 
the distribution of (123|) moves towards a more peaked distribution, and is actually well 
described by a Gaussian, with a volume-dependent width, shown in figure El Thus, it 
appears plausible that this quantity, and thus likely also the gluon propagator, will have the 
same value, irrespective of the representative chosen along the residual gauge orbit. 
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Figure 3: Lower limit of the width of the distribution 
of ([^5)) for a fixed discretization of a = 0.22 fm on a 
lattice [g^, after subtracting one and rescaling by a 
factor of -1000. 
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However, this statement does not take renormalization into account 87|. In particular, 
has to be regularized, to be meaningful in the continuum limit. It is not a-priori 
clear how this quantum effect affects the distribution, as it involves the regularization of 
a composite, non-local operator. From U(l) gauge theory it is known that this may cause 
significant problems joil ]. 

The consequences of changing the cut-off and performing only a multiplicative renormal- 
izatiorE^l are shown in figure HI While in two dimensions, where almost no Gribov copies 
are present, there is little change, the distribution appreciably changes in higher dimensions. 
Thus, regularization of this composite, non-local operator plays an important role in under- 
standing the thermodynamic limit, and how stringent a finite polynomial of the fields has 
to take the same value irrespective of the choice of Gribov copy. It is therefore better to 
to investigate instead the gluon propagator, and see whether it takes a unique value in the 
thermodynamic limit. This will be discussed at length below in section 14.1. 1[ 



2.5.3. Absolute Landau gauge 

An alternative way to choose a representative on the residual gauge orbit is the absolute 
Landau gauge 63|, [SJ, l97|, which makes a very definite choice rather than a random choice. 
This gauge choice is derived from the following observation. The functional fl23l) has, up to 
topological identifications 8J], a unique absolute maximunrj. The resulting set of absolute 
maxima, called the fundamental modular domain or region [97| , in analogy to conformal field 
theory, is by definition embedded in the first Gribov region, and includes the origin. It is less 
trivial to show that it is convex and bounded, and thus connected 65 . It can furthermore 



be shown that part of the boundary of the fundamental modular domain coincides with the 
Gribov horizon in the thermodynamic limit only 65|, |97(|. All possibly remaining degenerate 



^•^Notc that volume effects can also play a role here, and thus volumes as large as possible should be used, 
^''it should be noted that in some theories the relative and absolute maxima become arbitrarily close in 
the thermodynamic limit (os} . 
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Raw distribution of F(A) in two dimensions 



Raw distribution of F(A) in three dimensions 
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Raw distribution of F(A) in four dimensions 
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Figure 4: The change of ((23| when the cut-off is in- 
creased, and only a muhipUcative renormalization is 
performed [g^l- The top- left panel shows a volume of 
(5.7 fm)^ in two dimensions, the top-right panel a vol- 
ume of (3.1 fm)'^ in three dimensions, and the bottom 
panel a volume of (1.3 fm)^ in four dimensions. A one 
is subtracted before renormalization. 
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absolute minima are on the boundary. This boundary has actually a quite rough structure 



99|, including wedge singularities, and topological configurations, like e. g. instantons [44 



are located there. By construction, in this region the gluon propagator has its least integrated 
weight. A sketch of the situation is shown in figure O 

Based on this observation, the absolute Landau gauge is defined as selecting the Gribov 
copy which belongs to the fundamental modular domain 
realized by either checking the absolute minimiza tion of i 
of a suitable weight factor in the path integral 



100, 101 



102 



This condition can be 
explicitly or by the introduction 
In case the residual gauge orbit has 
more than one Gribov copy on the boundary of the fundamental modular domain, again a 



random choice is made [63|]. It should be noted that if the thermodynamic arguments made 
before were correct, the absolute Landau gauge and the minimal Landau gauge will coincide 
in the thermodynamic limit, at least for correlation functions being finite polynomials of the 
gauge fields. 



It should be noted once more that the relevant quantity fl23|l has to be renormalized. As 



21 




Figure 5: A sketch of the geometry of field configuration space. A gauge orbit is shown, which intersects 
the gauge-fixing hypersurface, defined by ([6]), six times in the visible part. Three of these Gribov copies are 
inside the first Gribov region. These three form the residual gauge orbit, after the restriction to the first 
Gribov region has been performed. As is visible, they are not connected by infinitesimal transformations. 
One of the Gribov copies is (necessarily) inside the fundamental modular region (FMR), and in this case a 
second Gribov copy is located very close to the boundary of the fundamental modular region, which itself 
touches the Gribov horizon. 

has been illustrated by figure HI this is likely not just a multiplicative renormalization. How- 
ever, an additive renormalization, which is in principle Gribov-copy-selection-dependent, 
can, in principle, make a local maximum a global maximum, and this is not yet fully un- 
derstood. In addition, finding the global minimum along the gauge orbit is in general a 
hard problem, since this requires to know all Gribov copies inside the first Gribov region. 
Therefore, actually reaching the absolute Landau gauge is in practice very complicated, see 
chapter [31 

2.5.4- Landau-B gauges 

A further proposal are the so-called Landau-i? gauges. They are motivated by two 
observations. 



One is that all information on a theory is contained in its correlation functions [33|. As 
a consequence, if two gauges differ, at least one correlation function has to differ at at least 
one momentum configuratioiJ^. E. g., the perturbative Landau gauge condition ([6]) can be 



^^To the extent that the correlation functions have to be analytic functions, this implies that they then 
differ for almost all momenta. 
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cast into a condition on the gluon propagator as (40|, |81| 



P.PuDt = {pM-p)p,AM) = 0, 



(24) 



for all momenta equally, and thus locally, and for each configuration individually. Thus the 
gauge condition in Fourier space already implies the vanishing of the longitudinal part of the 



gluon propagator. Similar conditions can also be posed for covariant gauges in general [4G 



It appears therefore possible that the non-local condition to select a particular representative 
of the residual gauge orbit should also be specifiable by a condition on one or more correlation 
functions. Since the conditions need only to be non-local in position space, they can be local 
in momentum space. 

The second observation comes from the solutions of the Dyson-Schwinger equatioii s 
(DSEs), as will be discussed below. There a family of solutions is found 37|, Il03l - ll06 



which are parametrized by the value of the ghost dressing function at zero momentum 



37 



or, equivalently, by the infrared value of an effective ghost-gluon coupling |103l4l05| . It has 



been speculated that this should also be possible at a non-zero momentum instead [37|, l81 



but this possibility has not been followed yet. Therefore, here only the zero momentum case 
will be discussed. Of cou rse, in a finite volume, and due to the trivial zero-modes of the 
Faddeev-Popov operator 
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the ghost propagator cannot be evaluated at zero momen- 
tum. Therefore, in this context always the momentum closest to zero will be used, and it is 
understood that the infinite-volume limit has to be taken eventually. 

The idea that this quantity could be of interest is supported by the fact that in almost 
all differin g va lue of (|23j) a l so le ads to a differing value of the ghost propagator at 



zero momentunQ B, IH, [lOdllS - lllO 



However, without knowledge of all Gribov copies, 
this cannot be made an exact statement at the current time. 

This observation can now be used to construct a set of gauges, the Landau-i? gauges 



The following construction has been found to work on small finite volumes [24 . l62 



and also at strong coupling [ill 



The situation in the thermodynamic limit will be 
discussed below. The basic ingredient of the Landau-i? gauges is a second, strictly positive, 
non-perturbative gauge parameter B, which can be selected within a certain range to be 
determined below. 



To implement the gauge, an auxiliary parameter b is needed, which is defined as 

hm 



(25) 



By definition, this quantity is renormalization-group invariant. In two and three dimensions, 
fi can actually be chosen to be infinite [81], since the correlation functions do not diverge, 
and then G{oo, oo) = 1, see section 13.2.31 The parameter b is taken to be defined Gribov 
copy-wise, i. e., on each individual Gribov cop}0. Thus, strictly speaking, b is, as F{A), a 



^^For simplicity, in the following always zero monientuni is used, with the understanding to replace it by 
the lowest momentum accessible. 

^^In a finite volume, the actual momentum is not zero. Therefore, the fact that a Gribov copy is not 
rotationally and translationally symmetric has to be taken into account in practical calculations, see (23 . 



6i,|8H 
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functional of the configuration, and should also be written as b{A) and the appearing ghost 
dressing functions as G(A, p^, u^). This will be suppressed throughout. 



The gauge is now defined [Sjj by the condition that on each residual gauge orbit the 
Gribov copy is selected as a representative which has a value of b, which is closest to the 
gauge parameter B. If there exist more than one Gribov copy satisfying this condition, once 
more a degeneracy is encountered, and in this before, a random choice should be 

made. The element of randomness has to be such that it selects in equal measure copies 
which have a positive or negative distance to B. This gauge construction implies that 
< b >= B, and thus determines the averaged ghost dressing function at zero momentum. 
If B should be selected such that this cannot be realized, e. g. because B is larger than b 
for all copies and configurations, the resulting ghost propagator is as close as possible to B. 



In this sense, it is possible to define a minB and a maxB gauge [81[, by always selecting the 
Gribov copy with the smallest and largest value of b for each orbit, respectively. 

It should be possible to construct a corresponding weight for the path integral, which 
implements this condition. But neither the existence nor the form of such a construction 



has yet been proven, there exists only proposals for it [62, ISO . 

Given the path integral in Landau gauge restricted to the first Gribov horizon fl2Tl) . one 
proposal to do so is by modifying the path integral to 

<0> = lim lim j VA^VcVcOiA^, c, c)e {-d^Df) e" Z'^*^^'' 

X exp l^-i 1 d'xd'ydlc\x)dlc'^{y)^ - ZbI^ \ . (26) 

Here, the representation of the ghost dressing function at zero momentum as 
G"^"(0) = \\myD'^a\p) = \\mv'{c\-p)c\p)) 

= lim-*- [ d'^xdVe'^^^~yHc''(y)c''(x)) 
p^o V J 

= lim^ I d'xd%dpye'P^'-y\(^{y)c%x)) 



- li^:^y" dVye'P^''-y^dy{c''{y)d^^c''{x)) = (^^^ j d'^xdSdie{y)dy[ 



has been used, and = B guarantees the correct renormalization. The volume factor 
V remains due to the translational invariance of the zero mode. The Gaussian weight thus 
inserted in f l26|) then enforces that in the limit the ghost dressing functions has to reach the 
pre-determined value, satisfying fl25|) . for each Gribov copy. Note that this additional term 
explicitly breaks perturbative BRST symmetry, and is non-local. Furthermore, the order of 
the limits is important. The statement of this inclusion is, however, merely that the ghost 
dressing function should become B at zero momentum, which can therefore be considered 
a boundary condition. This permits to implement this constraint also without explicitly 
taking into account this non-local term in the path integral in self-consistent calculations, 
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as will be done in section [3^ Of course, in lattice gauge theory just the appropriate Gribov 
copy will be selected. It should be noted that it is likely possible, analogously to conventional 
covariant gauges 0], to remove the limit in rj in (126|) . This yields then a gauge averaging 
over non-perturbative Gribov copies, with a Gaussian weight centered at Zb'J- 

However, the expression is only appropriate if, up to a measure zero set, all con- 
figurations have at least one Gribov copy which can satisfy the conditiorj^. Otherwise 
gauge-invariant quantities would be affected. If a i? value is chosen for which this is not 
the case, the maxB or minB gauges are explicit examples of such a case, the expression f l26|) 
cannot be correct. Instead, at best it can only be imposed that the condition < b >= B 
can be satisfied, but not b = B. This can be implemented in the path-integral by using a 
Lagrange parameteiEl A(i?), like temperature, as a gauge parameter, which gives the 



form [62, |80|, [87 



<0> = lim y VA^VcVcOiA^, c, c)e {-d^Df) e-/'^*^^^ 

X exp (^N{B) + A(S)i j dVydlc\x)dlc\y)^ , (27) 

where is a, possibly orbit-dependent, normalization, such that gauge-invariant quantities 
are not affected by the averaging. If the normalization is orbit-independent, it can be 
absorbed in the path integral measure. Otherwise the same statements about fl27|l hold true 
as for fj26|l . In general, the Lagrange parameter A may depend on global properties of the 
theory, in particular the value of the zero-momentum ghost propagator for A = 0. Indeed, 
for at least this one value A = this gauge condition is definitely valid, as it is just implies 
to average over the residual gauge orbit. In addition, it is trivially possible to define 

A = Z^^B, (28) 

and give the Lagrange parameter the meaning of an unrenormalized gauge parameter. Note 
that there are still some values which cannot be imposed, like negative values of S, much 
like negative temperatures (usually) cannot be employed. It thus still remains to identify 
the permitted set of B values. The minimal Landau gauge appears in such a setup at A = 0, 
since it is equivalent to averaging over all Gribov copies with a fiat weight. This implies 
that due to (125]) the minimal Landau gauge is a fixed-point of the renormalization group 
with respect to the both the perturbative gauge parameter of the covariant gauge and this 
non-perturbative gauge parameter. 

At the present time there is neither proof, nor rock-solid evidence that this eliminates the 
Gribov-Singer ambiguity in the same way as the absolute Landau gauge, i. e., by uniquely 
identifying a Gribov copy for at least some value of 5, which would imply for this/these 
value(s) of B the admissibility of the construction (12^ . However, when averaging over 
the residual gauge orbit using the prescription f l27|) the Gribov-Singer ambiguity is lifted 



^^I am grateful to D. Zwanziger for a discussion on this topic. 

^^Which on a finite lattice will additionally depend on the volume and the discretization. 
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by construction: At least in the same form as the minimal Landau gauge, it provides a 
prescription how to treat the residual gauge orbit. This follows, since if indeed there exists a 
set of residual gauge orbits with non-zero measure where there are Gribov copies degenerate 
in this condition, then, as in minimal Landau gauge, just random choices will be made. 
The gauge-fixing conditions then returns to the same operative level of gauge definition 
as for the minimal Landau gauge. Of course, in such a case the non-perturbative gauge- 
dependence between different Landau- -B gauges will be trivially just a lower limit to the 
variability which would be obtained when resolving any further degeneracies, like in the 
absolute Landau gauge. Thus, all three gauges, minimal Landau gauge, absolute Landau 
gauge, and Landau- -B gauges, are complete resolutions of the Gribov-Singer ambiguity: In 
all cases it is fully specified how to treat Gribov copies. 

After specifying this gauge, it is now an important question what the actual distribution 
of b is. A particular interesting question is, whether there are holes in the distribution of b, 
whether it has a finite width in the thermodynamic limit, and so on. Results from functional 
continuum calculation [37| suggest that the possible range is [Bq, oo), with some lower limit 
Bq not precisely determined yet, and including the limit B oo. The latter would be 
required to be implemented by another limit in fl26|l to be taken after the limit r/ — )■ 0, and 
also in fl27|) . However, even in the opposite extreme case that the interval should shrink to 
a single value, the construction is not wrong, since the ghost dressing function has to take 
some value at zero momentum, being it finite or infinite. In this case, all what happens, is 
that the Landau- gauges are reduced to the minimal Landau gauge. 

Thus, the permitted range of B is primarily of practical importance. The distribution of 
6 as a function of volume is shown in figure El First of all, the distribution is very strongly 
asymmetric, with a long tail towards large values, and thus quite different from the Gaussian 
one found for F{A). It is furthermore strongly volume-dependent. However, as can be seen 
from the result for two dimensions where on small volumes almost no Gribov copies exist, 
this is to quite some extent a 'trivial' finite volume effect, which is also seen in the minimal 
Landau gauge. These just stem from the fact that in a finite volume the lowest momentum 
is non-zero, but decreasing with volume. 

On the other hand, the renormalization properties shown in figure [7] are much better than 
for F{A). If the discretization errors are not too large, i. e., the number of lattice points 
is sufficiently large, then there are essentially no cut-off effects visible for the distribution. 
This is to be expected from a renormalization-group invariant quantity, which b is. From 
available results, it appears that a cut-off larger than 1 GeV is at least necessary for reaching 
the appropriate behavior, and at least a lattice size of 20^^. 

Hence, it remains to get a better view on the development of the permitted interval for B. 
A useful possibility is to determine the width of the gauge corridor [(min6), (max 6)], given 
by (max b) / (min b) . Taking the ratio removes any over-all rescaling effects due to trivial 
renormalization and volume effects. This ratio is shown in figure [81 The results shown are 
only approximations in two respects. Because of the inherent problem to find all Gribov 
copies with any given method, the results shown give only a lower limit to the size of the 
interval. That is an issue which has to be taken into account at any rate for any results 
displayed in this work: In almost all cases only lower limits of the possible gauge variations 
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Figure 6: The distribution of b for a fixed discretiza- 
tion of a = 0.22 fm [g^. Results are shown for two 
(top- left panel), three (top-right panel), and four di- 
mensions (bottom panel). The absolute value is nec- 
essary due to some peculiarities of the method in de- 
termining the b value in the algorithms employed here 



1121 Ill3l | , which are negligible when determining lat- 



tice averages, and in particular for the thermodynamic 
limit (ilfe^l. 



are available, and though an extrapolation in the number of Gribov copies is possible in 



principle yS, Ill4j . there is no guarantee. The effect of taking different search spaces into 



62, 87 



account is illustrated in figure M A detailed discussion can be found in 

Irrespective of this, the corridor in all dimensions, even in the case of two dimensions with 
its small number of Gribov copies, quickly opens with increasing volume. The increase of the 
ratio of upper to lower bound of the corridor is actually not affected by trivial finite volume 
effects, but can be affected by non-trivial finite volume effects. Still, this is the strongest 
impact of the choice of Gribov copies for the investigated volumes which has been found so 
far. Though it is tempting to conclude from these results that the interval of permitted B 
values is indeed larger than a point, and has an infinite upper limit, experience has shown 
that this cannot s o ea s ily b e decided without investigating the volume dependence over a 



much wider range [115| - |118| . In particular, convincing evidence has been foun d that at least 



in three and four dimensions the lower limit is necessarily finite |115l Ill6l . Ill8| , and thus the 



value of (minfe) has to flatten out. This effect is not seen yet. Furthermore, the dependence 
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Figure 7: The change of (1^^ when the cut-ofF is in- 
creased for /i = oo in two and three dimensions, and 
jjL — 2 GeV in four dimensions. The top-left panel 
shows a volume of (5.7 fm)^ in two dimensions, the 
top-right panel a volume of (3.1 fm)'^ in three dimen- 
sions, and the bottom panel a volume of (1.3 fm)"* in 
four dimensions 62 1. 
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is on the discretization is non-trivial. While in two and three dimensions the corridor width 
increases with finer discretization, this is not the case in two dimensions. This reminds of 
the case with the number of Gribov copies, shown in figure [H where two dimensions also 
behaves differently than three and four dimensions. 

Note that the construction so far has been entirely based on selection of Gribov copies, 
without returning to the findings in functional calculations. Whether indeed, as in functional 
calculations, the upper limit of the B corridor is infinite in the thermodynamic limit cannot 
be decided yet. The evidence in favor or disfavor of this possibility will be discussed at 
length in section 14.1.11 Some rather general comments on whether this is possible at all 
will also be given in section 14.3.21 However, the discussion here shows that the expectation 
value of the ghost dressing function at zero momentum can be brought into contact with 
Gribov copies. If either fl26p or (I27|) can furthermore be shown to be correct in general, this 
would imply that selecting the zero-momentum value of the ghost propagator in functional 
calculations is indeed equivalent to treating Gribov copies. 
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Gauge corridor width in two dimensions | Gauge corridor width in three dimensions 




Gauge corridor width in four dimensions 




The fact that giving a definite prescription how to deal with Gribov copies can be based 
on a desired behavior of the ghost propagator at zero momentum is already showing that 
the Gribov-Singer ambiguity can be resolved in this way in a well-defined manner, in the 
same sense as the minimal Landau gauge resolves it. It thus remains to formally understand 
how this is resembled in the approach starting with the continuum theory. 

2.6. Gauge relations 

Before turning to the more practical matter of calculating the correlation functions, it 
provides further insight to investigate the relation between the two quantities F{A) and b. 
A first insight is given by the projection of the first Gribov region to these two coordinates. 
This is shown in figure (TUl It is visible that, for the largest volumes studied here, there is 
essentially no correlation between the two coordinates. Thus, it can be expected that when 
using one of the parameters to fix the gauge, the results for the other parameter will be 
essentially its average value. Thus, it will have a similar value as in the minimal Landau 
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Figure 9: Various observables [H^] as a function of the size of the search space [sil llOC| for a 36'^ lattice 
with lattice regulator a = (1.6 GeV)~^. Note that also the number of configurations (here 1632) taken into 
account can modify this result. Lower left panel: Double gauge orbit and configuration averaged value of 
F{A). Lower middle panel: Double gauge orbit and configuration averaged value of h. Lower right panel: 
Number of Gribov copies found. Upper left panel: Width of the F{A) corridor. Upper right panel: Width 
of the B corridor. 



gauge. 

The next observation comes from classifying the different Gribov copies. One possibihty 
is to select the fundamental modular region copies, at least as good as the algorithm can 
identify them. The second possibility is to select the largest value of h. In minimal Landau 
gauge, for sufficiently large volume s, th i s qu antity is dominated by the lowest eigenvalue 
of the Faddeev- Popov operator 2% Illy, Ill9l |. and thus this would be in one-to-one corre- 
spondence with the Gribov horizon. If this were to hold in general, h can also be used to 
characterize the Gribov horizon. 

It is then seen in figure ITO] that the set of copies belonging to both categories shrinks with 
increasing volume and cut-off. Thus, less and less gauge orbits have a Gribov copy belon ging 



to both categories. Hence, it appears that though such a common boundary exists |65l. I97l|. 
it is only a small part of the total boundary of the Gribov region. But this may also be 
an artifact: The number of copies being found satisfying both criteria is found to increase 
when the search space increases \Q2\ , so there is the chance that this particular Gribov copy 
of each orbit is often just not found among the large number of copies. However, this would 
also imply that the standard algorithms used for minimal Landau gauge will not sample this 
type of copies very well. 
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I Two dimensions, V=(2.2 imf, a-'=4.47 GeV | | Two dimensions, V=(5.7 fmf, 3-^=1.73 GeV | | Two dimensions, V=(11 imf, 3-^=0.895 Gev" 




ln|b| ln|b| ln|b| 

Figure 10: The first Gribov region projected to the coordinates F{A) (after subtraction of one and multi- 
phcation by -1000) and b. Top, middle, and bottom panel show results for two, three, and four dimensions, 
respectively. Always three different volumes are shown, with increasing volume from left to right. Red points 
denote copies on the would-be horizon (see text), green points are in the fundamental modular domain, blue 
points have both these properties, and black points neither. More details can be found in j6^ . 



Because both coordinates appear to be quite unrelated, it is a natural possibility to select 
combined gauge choices. E. g., the Gribov copy for a configuration could be selected which 
maximizes 

xF{A) + yb, (29) 

for some value of x and y, or any other function of the two coordinates. However, it does 
not appear to be possible to require the satisfaction of two independent exact conditions for 
F{A) and b. Another possible choice would be just minimizing instead of maximizing F{A), 



which could be called an inverse Landau gauge [8jJ, Ill4| . Here, these possibilities will not 
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be pursued further, and only the set of Landau- i? gauges together with the absolute and 
minimal Landau gauge will be investigated. 



3. Methods 

It is actually a quite non-trivial problem to determine the correlation functions beyond 
perturbation theory. Various methods have been employed for this. However, much more 
successful than any individual method was the combination of two or more of them to 
compensate the respective disadvantages. For the task of describing gluons, the combination 
of lattice gauge theory and functional methods has been most useful, and will be employed 
here. Alternative approache s employed are in particular methods based on effect ive a ctions 

|121| . but also stochastic quantization approaches |65l. l94j . |122[ |. both 



of which will not be detailed further. 

The functional methods can be split further, in the quantum equations of motion on 
the one hand, and functional renormalization group methods on the other. Though the 
quantum equations of motion will be mainly used here, the related approach using functional 
renormalization group equations will be shortly introduced in section 13.2.21 

An enormous advantage of all of these approaches, in contrast to models, is that pertur- 
bation theory is always included. Because of asymptotic freedom, it becomes manifest for 
Yang-Mills theory at large momenta. Conversely, this implies that use can be made of the 
wealth of perturbative results in all of these methods. 

3.1. Lattice gauge theory 

One of the most successful methods to investigate gauge theories in the non-perturbative 
domain is lattice gauge theory ji], liif. In this approach, space-time is discretized and only a 
finite volume is investigated, usually with the geometry of a hypercube having N sites in each 
direction space in intervals of length a, and having thus a physical volume a'^N'^ = L'^ = V. 
Effectively, thus, quantum field theory is approximated by a finite system and thus quantum 
mechanics. Only by taking the appropriate thermodynamic limit in the end, the original 
quantum field theory is recovered. However, the existence of this limit is still an unproven 
assumption. The results shown in figures [T]|l] and [UlfTOl in the previous chapter were all 
obtained with this method. 

This method is particularly useful to determine gauge-invariant quantities 0,1431. Fur- 
thermore, in this case it is known how to establish the connection to Minkowski space-time 



3l| . The most dramatic successes of lattice gauge theory have been obtained by making 



use of the fact that it is very well approachable by numerical methods [43|, permitting to 
evaluate path integrals directly. Though this process just delivers final numbers, it permits 
insight into the hadronic spectrum and many other results, and is essentially only limited 
by the possibilities of the numerical algorithms. Though these made enormous progress over 
the years, there are still substantial challenges. These include very light dynamical fermions, 
isospin breaking, finite density, disparate scales, supersymmetry, critical slowing down, and 
other problems. 
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Here, however, the focus is different and is on the possibihty to determine correlation 
functions directly by numerical evaluation. With the numerical approach to lattice gauge 
theory, it is only possible to evaluate full correlation functions, defined by ([8]). To obtain 
vertex functions, the connected contributions have to be determined, and these have to be 
amputated. The latter has quite some implications in Landau gauge, as will be discussed in 
more detail in section 13.31 

3.1.1. The action 

The detailed approach to obtain the correlation functions is a multi-step process 43 



Usually, the formulation of Yang-Mills theory using lattice gauge theory is in terms of the 
group elements rather than of the algebra-valued gauge fields themselves. In principle, any 
group representation of the theory is equally valid, however, some turn out to be more 
amendable to numerical calculations than others. A dramatic example is the case for the 
gauge algebra su(A^). In this case a formulation in terms of the group SU(iV) is rnuch mor e 



efficient with the available algorithms [43|| than for the gauge group SU(A^)/ZAr [I23l4l25 
which would be more appropriate for the standard model pj]. Since the gluonic correlation 
functions are invariant under group transformations which are mapped to the identity in 
the algebra, e. g. center transformations for SU(A^), the notion of algebra and group will be 
used synonymously in the following. 

In general, the gauge fields v4^(x) are thus represented by the corresponding group- valued 
link variables U^{x) = exp(iA^rtt) for each lattice site x and direction fi. Of course, from 
these the original algebra-valued gauge fields can always be reconstructed, though this in 
general produces errors of the size (9(a^), where a is the lattice spacing. E. g., for SU(2) the 
relation is given b}@ 



A; = ^tTT'^{U,-U^) + 0{a') (30) 



where /3 = ANp/g'^ is related to the bare coupling g, which has dimension a~^~, and Np 
is the dimension of the fundamental representation of the gauge algebra. Similar relations 
hold for other gauge algebras. 

The discretized lattice action Sl is then given in terms of these link variables. The form 
used here is the so-called Wilson action M 



X,U>fl \ ^ / 



Uf,, = U^{x)U,{x + ef,)Uf,{x + e,)+U,{x)+, 

with the unit vector in the hypercubic lattice in the direction fi and the (gauge-invariant) 
plaquette f/^^. The size of the hypercubic lattice is A^"*^. This lattice action agrees with the 
continuum action S up to corrections of order a^, and thus coincides with it in the limit 



^°Care has to be taken if the center symmetry of the gauge group is broken [126[. In this case, (pO| can 
at best be used for configurations with a positive, real Polyakov loop. 
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a — )■ 0. This behavior can be improved l43|, but this will not be ne cessa ry here. Studies of 
correlation functions using improved actions can be found, e. g., in |l27l . 

Using a Monte Carlo simulation, a set of (independent) gauge orbits {G}, being a set 
Q = {U^} of dN"^ link variables each, is generated. For this purpose, various methods 
are available, which will not be detailed here, as they are not specific for the purpose of 
calculating correlation functions. See 43| for a comprehensive introduction, and for the 



particular algorithm used for most of the results presented here [112 



It should be noted that there are quite a number of technical problems still associated 
with the generation of the configurations, not the least to determine the value of the cut-off, 
and thus of a as a function of /3. For the present purpose, standard methods are sufficiently 
well known j^, |4^, and will not be detailed here. The important technical point is that in 
all cases here the scale is set by giving the intermediate distance string tension the value 
(440 MeV)2, see for details e. g. 0, H [lij. 



3.1.2. Fixing the gauge 

A typical Monte-Carlo algorithm generates a set of Nc independent field configurations 
[Ufj] using the weight exp(— S'l), i. e., each configuration contains link variables for all 
directions and lattice sites. Each of these link configurations are in a random gauge, and have 
to be transformed into the desired gauge first before the gauge-fixed correlation functions can 
be evaluated. Otherwise, the gauge-dependent correlation functions av erag e to correlation 
functions of Gaussian random variables, according to Elitzur's theorem 3J . 

The gauge-fixing to satisfy the the Landau gauge condition ([6]) and to restrict to copies 
inside the first Gribov region can be performed by minimizing the expression 



Y,^tU,{x) ~ 1 + const. ^ A^Al + 0{a^), 



(31) 



x,ii,a 



which agrees with (123|) in the continuum limit up to prefactors [96]. For this purpose, a 
gauge transformation g{x) on the lattice is searched for, such that the gauge fields obtained 



from the gauge transformed links 



g{x)U^{x)g{x + e^)"*" satisfy the Landau gauge 



condition. To check whether this condition is indeed fulfilled, it is possible to calculate, e. 
g., the gauge- fields by (15(1 and check whether they are transverse [96|. In fact, despite that 
these fields can be transverse up to the numerical precision, due to the discretization errors 
the corresponding continuum fields will only be transverse up to order (9(a^). 

Actually, this direct local evaluation is not very sensitive to long-wavelength fiuctuations. 



Therefore, a more useful quantity to evaluate is [96 



(32) 



^ Yl l9ix)Uf,{x)g{x + e^y - g{x + e^)Uf,{xyg{ 

Xv,l^=ii^J, 



X 
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where is the size of the hypercubic lattice in direction a. A t ypical good numerical 



condition on this quantity is to require it to be below 10 [96|, Ill2 



To determine the gauge transformation g{x), several algorithms have been developed, 



see 



96 



for a detailed discussion. A particularly useful example of such an algorithm for 
contemporary lattice sizes is the stochastic overrelaxation algorithm. Essentially, it updates 
g{x) at each lattice site in turij^ by either choosing the new g{x) such that it minimizes (131 p 
locally or by only rotating it such that fl3T|) is left unchanged. Which of both possibilities is 
performed is chosen randomly. The probabilities assigned for both possibilities significantly 
influence the performance of the algorithm, and depend on the volume, discretization, and 
dimensionality of the lattice investigated 90|, |96|. It is m ost convenient to permit the 



algorithm to adapt these probabilities during run-time Ill2l| . However, the best value for 
the probability is also very much configuration-dependent [63|. Therefore, no optimal choice 
exists, and the best choice is one which optimizes the algorithm on the average. 

The gauge fields obtained so minimize (13T|) and therefore fulfill the perturbative Landau 
gauge condition (jS]). Furthermore, it can be shown that they are in the first Gribov region 
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by noting that at the minima of (IHT]) the Hessian of its second derivatives has to be 
positive (semi-)definite, and that this Hessian coincides with the Faddeev-Popov operator 
( 120|1 . Actually, getting out of the first Gribov region using lattice calc ulations is not entirely 
trivial. This can be achieved using, e. g., stochastic quantization 129| or gauge- fixed Monte- 



60, 130 



Carlo simulations 

However, without further specification no minimum is preferred and thus a random 
Gribov copy is obtained within the first Gribov region. Choosing this copy for proceeding 
is therefore equivalent to implementing the minimal Landau gauge. To obtain other gauges, 
it is necessary to specify this copy further. Since for this purpose knowledge of more copies 
is necessary, this requires access to further Gribov copies. Up to now, no systematic way 



of genera ting new Gribov copies is known. Therefore, a stochastic approach is used [81 



lOOl Ill4| . For this it is important to note that any kind of local algorithm which attempts 
to minimize fl3T|) will find the closest minimum to its starting point, or even with some 
amount of randomness in its search, a minimum very close to its starting point. Therefore, 
restarting the gauge-fixing algorithm with a different initial guess for g{x) will in general 
produce another Gribov copy, if there exists more than one copy, though with an unlucky 
guess the previous one will be founds. By repeating this sufficiently often, a set of Gribov 
copies is generated. In this set of copies then the copy which matches the non-perturbative 
supplemental gauge condition, e. g. absolute Landau gauge or a Landau- -B gauge, as good as 
possible is selected to represent the gauge orbit. Since this is a high- dimensional optimization 



problem, it is in general not possible to guarantee finding all copies [63[. Therefore, gauges 
like the absolute Landau gauge or the maxB gauge can never be reached with certainty, 
though sophisticated algorithms can be devised to make the search at least somewhat more 



^^This can be accelerated by using a checker-board update instead of a lexicographical update. 

^^This assumes perfect numerical accuracy. On finite machines, numerical artifact copies will be found, 
but can be identified with knowledge on th e precision of the algorithm 63| . Further artificial Gribov copies 
exist due to the lattice regularization jl31 *|, but these become irrelevant when taking the continuum limit. 
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efficient |63l . Il09l . Ill4l . |132H134| . For gauges which do not strive to find an extreme value, hke 
a Landau- -B gauge with some non-extremal value of B, this problem only manifests itself in 
the generation of more stochastic noise. 

Any of these methods finally selects one particular Gribov copy with particular, gauge- 
fixed links for each configuration. This yields the set of gauge- fixed configurations {{m^}}. 
These represent the gauge orbit@ Note that no ghost fields need to be introduced. Gauge- 
fixing on the lattice, in the original spirit of gauge-fixing, is due to selecting a representative 
for each gauge orbit. 

3.1.3. Correlation functions 

Once such a set of Gribov copies representing the gauge orbits has been specified, cor- 
relation functions can be determined. This is done by determining the correlation function 
on each gauge orbit's representative and then average over all gauge orbits. Thus, for some 
correlation function - being the expectation value of the product of n fields A^, the 
lattice result for the correlation function from Nc gauge orbits is given by 

fi^-^" = = ^E^i-^- (33) 



c 



where the index i labels the representatives of the Nc gauge orbits. In the limit of Nq oo 
the function obtained is indeed the correlation function |9|. At finite Nc^ statistical fiuctu- 
ations will be present, and a stochastic error can be assigned to the correlation functions 



at each point [9|, |136| . Furthermore, on any finite lattice, there will be systematic errors 
due to the finite volume and discretization and the presence of the non-trivial geometrical 
structure introduced by choosing a hypercube. Only by taking the thermodynamic limit (or 
extrapolating to this limit) the continuum correlation functions in an infinite volume can be 
obtained. These lattice artifacts will be discussed in more detail in section 13.1.41 

Here, the correlation functions in momentum space are in most cases more interesting 
than those in position space. Therefore, instead of evaluating the correlation functions first 
in position space and then Fouri er-tra nsforming it is more convenient to transform the gluon 
fields directly, which is done by jl37 



A;:(p) = e— J^e'^^^''— A;:(x), (34) 



X 



where on a finite lattice the components of P have the integer values — iV^/2-|-l , . . . , 
and the components of X are the (integer) coordinates in the lattice ranging from to 
A*"^ — 1. Due to the periodic boundary conditions in position space, the fields are also periodic 
in momentum space, and thus only the momenta O...A^^/2 are independent. 



^■^If the aim is to implement gauges averaging with a specific prescription over gauge copies, like general 
covariant gauges or Hirschfeld-like gauges, this requires to average over the set of Gribov copies found with 
the required weight functions and over all Gribov regions. Such gauges will not be considered further here. 
See, e. g., (ifil. M. IHl. IstI. flM fl35t. 
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The physical momenta p to which the integer lattice ones P correspond are 



2 P.TT 

p, = -sin^. (35) 

The prefactor in fl34|) comes from the implicit mid-point definition of the gauge-fields when 
using the links, and has to be taken into account. The gluon propagator, e. g., on a single 
configuration is then given by 

K(p) = ^^WKi-p)^ (36) 

where use has been made of the reality of the gluon field in position space and V is the total 
volume of the lattice, stemming from the normalization of the Fourier sum. The three-gluon 
correlation function is then given by 

1 ^ 



V 



A;ip)Atiq)A;ik), (37) 



with p + q + k = 0, and so on. 

A bit more complicated is to determine quantities involving the ghost fields. In the 
lattice calculations a gauge copy is selected directly, and the Jacobian is never required, and 
thus no ghost fields appear. Therefore, the only possibility to obtain correlation functions 
involving ghost fields is by returning to the path integral ([7]). Integrating out the ghost 
fields gives expressions for correlation functions involving the ghost fields in terms of the 
Faddeev-Popov operator (120|) . 

To make use of these relations requires to determine the Faddeev-Popov operator on the 
lattice. It is defined in terms of its action on an arbitrary function a;" as a function of the 



gauge-fixed link variables as [97 

M{y, xr'uj,{x) = cIyI iG''\x)uj,{x) + Af{x)ujb{x + e^) + B;\x)uj,{x - e^)) J 
G^\x) = J2^Ti{T^,T'}{u,{x) + u,ix-e,))) 

A^J'ix) = -2tr(rV^M^(x)) 
Bf{x) = -2tr(rVn+(x-e^)), 

where c is a normalization parameter, depending on the normalization of the generators t"' 
of the gauge algebra. 

E. g., integrating out the ghost fields for the correlation functions of a ghost and an anti- 
ghost, the ghost propagator, then yields its representation in terms of the Faddeev-Popov 
operator as 

< c^ix)c\y) >=< M"^-^ (x, y) > . (38) 

Therefore, to determine the ghost propagator, and actually all correlation functions involving 
a ghost and an anti-ghost field, requires to invert the Faddeev-Popov operator. This needs 
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to be done numerically in lattice simulations. Similar to the case of fermion fields 43 



there are two nume rical possibihties widely used, both based on a (bi-) conjugate gradient 



algorithm |l38l . ll39| . exploiting the positive-semi-definiteness of the Faddeev- Popov operator 



in the first Gribov region^_Tliey differ by the source vector, point source [24, Ill2l Ill3| or 
plain- wave source lOOl llOTI . Ill2j . The advantage of a point source is that the inversion 
has to be done only once to get the full Fourier spectrum of the ghost propagator after 
Fourier transformation, while this has to be done separately for each mode with a plain- 
wave source. The drawback of the point source is that it suffers from significantly larger 
statistical fluctuations, limiting its use essentially to only calculating the propagator. A 
quantitative compariso n of both methods for the ghost propagator and the ghost-gluon 
vertex can be found in 112| . 

The ghost number is a conserved quantum number in Landau gaugj^. Hence, knowledge 
of the inverse Faddeev-Popov operator is then already sufficient to construct all higher order 



ghost correlation functions. E. g., the ghost-gluon vertex is given by [143 

< A^^{x)c^{yy{z) >=< A';,{x)M''^-\y,z 
Once more, this can be proven by integrating out the ghost fields 



> 



(39) 



Thus 

calculation of the propagators 



in general, the constr uctio n of vertices is a straight-forward generalization of the 

the relation fl55]) already shows that the 
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However, 

correlation functions obtained are the full correlation functions. To compare it to the vertex 
functions, lattice results always have to be reduced to their connected and amputated part. 
The connected part is not a problem for correlation functions up to order three, as any 
one-point correlatio n fun ction vanishes in Landau gauge, but can become very cumbersome 



at higher orders |90l . Ill2| . The fact that the correlation functions are non-amputated implies 
that only the tensor structures transverse in all gluon momenta can be obtained in lattice 
calculations in Landau gauge, but not the longitudinal ones. To isolate then the individual 
transverse tensor structures requires projection of the connected correlation functions with 

and then amputation by division with appropriate products of 
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appropriate tensors 

pr opag ators. Details of this procedure can be found in 143| for the ghost-gluon vertex and 



m 
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for the three-gluon vertex. Higher vertices have not yet been determined in lattice 
calculations, mainly due to the statistical noise which increases quickly with the number of 
fields involved. 

It should be noted that the correlation functions of an individual configuration will not 
exhibit translational nor rotational symmetry, and thus a propagator will depend on both x 
and y separately, instead of only the difference x — y. Only after performing the average f l33|) 
these symmetries will become manifest in the limit Nc — oo, within discretization errors 
and hypercubic artifacts. Thus, the averaged propagators will indeed depend only on x — y. 

3.1.4- Lattice artifacts 

The previous description produces the correlation functions for a finite set of momenta. 
This set ranges from a lower cut-off, which is essentially given by 1/L, where L is the 



There is a second symmetry of the ghost fiel ds in Lan dau gauge which relates ghosts and ant i-gho sts 
[l3j . In principle, this symmetry could be broken 140l 141 1, but this does not seem to be the case |l42j . 
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extension of the lattice in the direction of the momentum, up to a largest momentum of 
size 1/a, where a is the lattice spacing. At the infrared cutoff, the system becomes aware 
of its extension, and the fact that it is a cube, while at large momenta it becomes aware of 
its hypercubic symm etry, and that there is a lattice. Thus, in principle, only for momenta 



fulfilling the relation |145 



Y < P < -, (40) 
L a 

the lattice can be expected to give reliable results. Thus the strength of the lattice formu- 
lation lies in an intermediate momentum window. An example for such artifacts in the case 
of finite volume, i. e. finite L, will be given in section 14.1.11 in figure [Ml below. 

Though not directly comparable, an example of how misleading such artifacts could be 
comes from the treatment of full QCD. There, for several hadronic observables the results 
are still far from the experimentally known values because the quar ks are heavy, despite that 



the pion mass is already quite close to the physical one (see, e. g., |146l | for a more detailed 
discussion and further references). Also, expected states still have not been seen, like pion 
scattering states. In this case also physical effects, like isospin breaking, may play a role. 
Thus, this should serve only as a warning, but should not be treated as a general permission 
to distrust the lattice results. 

A less obvious source of artifacts appears when ?7,-point vertices, instead of propagators, 
are considered: The n — 1 independent momenta cannot be positioned arbitrarily with 
respect to each other, but only in a way compatible with the geometry of the lattice. E. 
g., on a two-dimensional lattice, there are very few possibilities to construct two momenta 
such that they have the same size as their sum, as the 60° angle between them can only be 
realized for very particular sizes in lattice units [90| . Thus, it is not always possible to easily 
reconstruct momentum configurations from the continuum, and this becomes worse the 
larger the number of momentum vectors involved is compared to the number of dimensions. 
Of course, the finer the lattice is, the weaker this problem becomes, but at lattice sizes 
currently accessible, this can be a serious restriction^^ 

There are further problems. One is that, since the lattice is discrete, any function is 
actually only given by a set of numbers. Thus, it is a-priori impossible to predict how well 
the results will approximate the continuum correlation functions in the complex plane, as a 
general statement of function theory [32]. This is of particular seriousness when considering 
the continuation to Minkowski space. A second problem is that there is still no proof that 
lattice gauge theory has as limit the desired continuum theory jol , though it appears unlikely 
that this should not be the cascj. 



^^This is connected to the problem of assigning spin and orbital angular momentum to states, since the 
rotation group is the discrete hypercubic group H^, rather than the continuous group SO((i) [43| . However, 
this is of little relevance for the topic here. 

^^In fact, at times the position is taken that the lattice is the only well-defined version of Yang-Mills 
theory. However, given that the universe is finite, and gravity plays a role at some point, this is more a 
philosophical question. Even if it were the case, the continuum field theory could then still be seen as a 
possible approximation of the lattice theory. 
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A particular problem for the gauge-fixed setting followed here is that in addition to 
the conv entio nal Gribov copies there are also additional Gribov copies which are lattice 



artifacts |131| . These will disappear when taking the limit to the continuum field theory in 
an infinite volume, but may contribute otherwise. However, it appears that this does not 
provide a new qualitative dimension to the problem of Gribov copies and in the present 
treatment lattice Gribov copies are therefore not distinguished from other finite volume and 
discretization artifacts. 

On top of all these principal problems, there is of course a practical problem. Most non- 
perturbative insights, including those presented here, come from numerical calculations. 
Thus, computational resources are a very practical li mita tion, especially when considering 



accessible volumes and discretizations |90|, Ill6l - lll8l . Il45l |. Thus, the investigation of dis- 
parate energy scales is strongly restricted by computation time. A possibihty is to use 
extrapolations, but since neither finite volume nor finite discretization are approximations 
which are under analytical control, extrapola tions cannot be considered to be ultirn ately 
reliable. It is known from solid-state physics 147l | and lower-dimensional QED 148| that 
these can be easily misleading. 

Hence, the lattice approach is alone not sufficient, at least at the current time, to treat 
the problems of interest, e. g., in the standard model and beyond. To complement it, in the 
following continuum functional methods will be paired with this approach. Their particular 
strengths are exactly at the extremities of scales: in the far infrared and ultraviolet and for 
light and heavy states, as well as combinations of all this. Thus, they are best where the 
lattice artifacts are strongest, and both approaches naturally compensate their respective 
weaknesses and combine their strengths. 

3.2. Quantum equations of motion 
3.2.1. Formulating the equations 

The quantum equations of motion, or Dyson-Schwinger equations (DSEs), can be ob- 
tained for full, connected, or vertex functions. The starting point are the respective equations 
for the generating functionals f lT2|) . f lT^ . and f|T5|) . Here, the most useful one will be the one 
for vertex functions f[T^ . though f[T^ will be valuable in section 1X51 

The equation f lT^ yields an equation for the one-point vertex function for the selected 
field. Repeated derivations with respect to other fields then generate a hierarchy of equations 
for higher- n-point correlation functions by virtue of ( 1T6|) . This der ivation is essentially an 



algorithmic problem, and can be formalized and thus automatized |149l . Il50l |. It shou ld b e 



noted that the equations may differ in form depending on the order of derivatives |149 
though of course this only corresponds to rearrangements due to identities relating different 
vertex functions. 

Before being able to obtain a set of equations it is first necessary to perform the gauge- 
fixing. To enforce the Landau-gauge condition the usual Faddeev-Popov formulation ([7]) can 
be used [ij. Thus, there will be not only gluon fields but also ghost and anti-ghost fields, 
and both types of fields have to be included in the construction of correlation functions. 
This fixes the perturbative part of the gauge. 
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The further restrictions are less simple to implement. First of all, it has been shown that 
the form of the equations is invariant whether they are derived on the whole of field- space, 



or within a fixed Gribov region [87|,l9J]. The reason is that the Faddeev-Popov determinant, 
which appears after integrating out the ghost fields, is by definition zero on the boundaries 
of these regions, effectively cutting off the integration. 

However, this is not sufficient to fix which of the Gribov regions is used. One possibility 
is to assume that to replace the 9 function in (|2T]) by a 5 function [67] is a valid procedure. 



This possibility will not be pursued here, see |6J, I67H79| for investigations based on this 
assumption. 

An alternative is to enforce the selection of the Gribov region by selection conditions on 
the set of vertex functions {F^, ...,F", ...}, which satisfy the various non-local gauge condi- 
tions. As addressed in section 12. 5[ there is not yet a sufficient condition known to restrict 
the set of solutions to those belonging to the first Gribov region. Here, the assumption will 
be made that it is sufficient to require that the ghost two-point function F'^^ is negative 
semi-definite in both position and momentum space. The remaining solution manifold can 
then be reduced further. E. g., implementing the absolute Landau gauge ( 12^ corresponds 
to selecting the solution of the DSEs which has the smallest integrated diagonal gluon prop- 
agator. Implementing a Landau- i? gauge can be done by implementing fl25|l as a restriction 
for the ghost propagator, i. e., requiring it to coincide with the desired value of -B, 

B.u„,4££(44). ,41) 

Thus, also B appears here explicitly. Note that this prescription has the same effect as the 
Landau- -B gauge condition, in particular in contrast to the condition (125|) . this is a condition 
on the averaged ghost propagator and not on an individual configuration. Hence, at this 
level this cannot be a precise statement about non-perturbative gauge-fixing in a functional 
approach, despite that the same results are obtained. Of course, if either of the conjectured 
restrictions fl26l) or fl271) were correct, and DSEs could be derived despite the non-locality 
involved, this would be an alternative. This will not be done here, but is should be noted 
that a modification of type fl7r|) . because it is bilinear in the ghost fields, will only appear 
in the DSE of the ghost propagator, and will only contribute to the position-space averaged 



ghost propagator [62 



It remains to formulate the case of the minimal Landau gauge. The only proposal of 
constructing an equivalent gauge choice with continuum methods is so far based on (127|) 
with the Lagrange parameter set to zero. Assuming this to be to correct and furthermore 
to be equivalent to the above ansatz in equation (I4ip . there is then an indirect way: In 
this case, in fact a Landau-i? gauge is chosen, but the value of B is the one obtained in the 
minimal Landau gauge in lattice calculations. This, of course, requires an outside source like 
lattice gauge theory for this information. Since the minimal Landau gauge just produces the 
average value of the ghost propagator, it will yield the average value of B, thus giving an 
indirect handle on fixing the minimal Landau gauge in DSEs. This prescription is correct, 
even if the ghost propagator coincided for all Landau- i? gauges in the thermodynamic limit. 
If, however, the Landau-S gauges did not yet completely exhaust all permitted constraints 
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to select a Landau gauge, this would imply that further constraints on further correlation 
functions can be implemented as well, and would have to specified if they were averaged 
over as well. This procedure leads to results in agreement between lattice and functional 
calculations, within the approximations made, as shown in section I4.1.1[ 

A practical implementation of the Landau- i? gauge, and thus the minimal Landau gauge, 
can then be obtained using the DSE for the ghost propagator. Since it is also the simplest 



DSE, it will now serve as an example for deriving DSEs [13|, 133 . 

Neglect for a moment the problem of Gribov copies. The ghost sector of the action is 
then given by 

S,j, = I d''z-c\z)d,{5'% + gf^'Al{z))c\z). (42) 
Inserting this into f|T5|) and differentiating^ with respect to (f{x) yields 

(^Z,d^^c\x) - gr''^dlAl{x)c\x) + -^^^ = 0, (43) 

where the x-index on a (9 indicates the variable with respect to which to derive, and the 
wave-function renormalization of the ghost with the wave-function renormalization constant 
has been taken explicitly into account. The renormalization constants of the interactions 
are absorbed into the derivatives of F, and will only reappear if these are written in terms of 
vertex functions, and the coupling constant is already the renormalized one. Replacing the 
fields by their respective derivatives and dividing, after taking the derivatives, by exp(Vr), 
yields 

Z d^^c^(x) - a r'^d^ (J]^JK_ + ] +JI- = (U) 

' ^ ' ^ \Sf^{x) 6r]<^{x) ^ 6jf,{x)6r]^{x) J ^ 5c-(x) ' ^ ' 

where j and 77(77) are the sources for the gluon and the (anti-)ghost, respectively. As a general 
feature of such derivations, products of single derivatives of W appear. When deriving such 
terms again only once with respect to the fields, always at least one single derivative remains, 
which can be replaced by a classical field. When setting the classical sources to zero at 
the end, the classical fields are also set to zero and therefore these terms always vanish. 
Hence they can already be neglected at this stage of the calculation, and will not appear 
furthermore. Note that this is not necessarily true if more than one further derivative would 
be applied. 

So the remaining expression is 

z3a-c"W_,/-a;^^^||^ + ^ = o. (45) 

To obtain the equation for the ghost propagator, this equation is derived once more with 
respect to c^(|/) which leads to 

^ ^' 5c^{y)5jl{x)5f]'^{x) 5c^(?/)5c"(x) ^ ' 



^^AU derivatives used are left-derivatives. 
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Figure 11: The DSE (|^5)) for the ghost. Dotted hnes are ghosts, curly lines are gluons. Lines with a dot are 
fully dressed propagators, white circles are full vertices, and all other vertices are bare. 

The last term will finally lead to the appearance of the ghost propagator, since W is the 
Legendre transform of F. The second term of fHB]) yields the interaction part. Using 



where the minus arises due to the anti-commuting derivatives. Since all mixed two-point 
functions vanish in Landau gauge, it is possible to write down the result in position space 

+gr''dl I d'zd'wDliix - z)D'j{x - w)rf^-'''^f{y,w,z). 

Herein the gluon propagator and the ghost-gluon vertex appear. 

Replacing all expressions with their Fourier-transformeco and afterwards dropping 
/ d'^p/ {27iY exp{—ip{x — y)) produces the result in momentum space as 

where the tree-level ghost-gluon vertex 

Tt'''^P.q.k)=^gr'^q, (48) 

has been made explicit. The resulting equation can be depicted graphically, similarly to 
Feynman graphs. For the case of the ghost equation, this is shown in figure [TTl 



^^All momenta are always defined incoming and momentum conservation at the vertices is taken into 
account. Hence in principle one of the arguments of the vertices could be dropped, but since this depends 
on conventions, all are kept. 
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Figure 12: The DSE for the gluon propagator. 



The DSE for the ghost propagator (148|) already shows the generic features of a DSE 
for a primitively divergent vertex function. These differ from the ones for non-primitively 



divergent vertex functions only by the appearance of the tree- level term [33l . ll51| . in this case 



the bare ghost propagator. The further terms are self-energy contributions, which have the 
form of loop integrals like in perturbation theory, but differ from them by the appearance 
of the full vertex functions. Thus, these equations are self-consistency equations and take 



the form of non-linear integral equations [33 



Another property is the generic structure of the equations, which is 



pn pn+l -p 



n+2\ 



(49) 



where it has been used that = in Yang-Mills theory and Fq denotes the tree-level value 
of F". The important property is that F" does not only depend on all vertex functions with 
lower n, but in general also on those witlj^ n + 1 and n + 2. Thus, this set of equations, 
being the analogue of the classical equations of motion, is coupled. In fact, it forms an 
infinite hierarchy of equations for the correlation functions. 

Furthermore, the number of self-energy diagrams proliferates quickly with the number 
but at most two-loop diagrams will be obtained. In addition, in contrast to 
the self-consistency of the equations guarantees that the equa- 



of legs j33l 1149 



perturbation theory [33|, 1152 



tions of motion are still well-defined, though their practical treatment becomes increasingly 
complicated. Already the gluon equation, shown in figure [TH gives an impression of the 
complexity. An important aside is here that the sign of the ghost-loop for tree-level quanti- 
ties is opposite from the one of the gluon loop at one-loop order. Thus, the gluons provide 



^^In a theory with higher than quartic terms in the fields even higher n-point correlation functions can 
appear. 
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Figure 13: The DSE for the ghost-gluon vertex. 



an ant i- screening behavior at large momenta, while the ghosts provide a screening behavior. 
Adding both up, this will finally yield the asymptotic freedom of Yang-Mills theory, as the 
gluons dominate at large momenta 
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The simplest vertex equation is the one for the ghost-gluon vertex, shown in figure [T2i 
The particular structure of this vertex will play a great r ole in the following sections. Various 



further examples for DSEs for vertices can be found in ITsol. [l53l - [l55j . 



Returning now to the problem of Gribov copies, the form of the ghost propagator DSE 

Multiplying the equation ( HHj) 



suggests a way to implement the Landau- -B gauges [37|, |81 

by p'^ and neglecting for now the color structurj^. it becomes an equation for the ghost 
dressing function. By subtracting the ghost equation from itself at zero momentum, it takes 
the form 

^ ^ n^"'=(p) -n^""(o), (50) 



where n"^ 
yields 



Gip) 

is the ghost self-energy. 

G(/i) 



G(0) 

Multiplying then by the ghost dressing function at fi 



ff=(/x)) . 



(51) 



Since the value of G{fi) is fixed by the renormalization condition, this completely specifies the 
equation, and introduces the B parameter directly. As noted above, this method can also be 
used to solve the equations in the minimal Landau gauge. Only the absolute Landau gauge 
requires a different approach. In this case, the full solution manifold of the equations has 
to be constructed, and then the one solution has to be selected which absolutely minimizes 
the (renormalized) expression fl23|) . 

The remaining problem is then to find a solution to the DSEs with a positive semi-definite 
ghost propagator. For such non-linear integral equations this is quite a formidable task, and 
rarely possible exactly. In particular, it is in general necessary to truncate the DSEs, which 
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This win be justified by the resuhs presented in section 2) 
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usually implies to neglect most of the equations for the higher correlation functions. This will 
be detailed below in section 13.2.51 Before turning to this task, it is useful to first introduce 
a related formalism, the (functional) renormalization group equations. 



3.2.2. Renormalization group methods 

The functional renormalization group equations (FRG(E)s) [39, 42, 156 | are very similar 
in nature to the quantum equations of motion. Without going into further details, they are, 
essentially, functional derivatives of the quantum equations of motion. 

The most important element in their construction is to introduce regulator functions. 
Since physical observables in a renormalized theory are independent of such regulators, this 
independence can then be used to derive equations for the correlation functions. A solution 
is then obtained in the limit of removing the regulator at the end. If th e system of equations 
is truncated, the choice of the regulator can have a non-trivial impact 157 . Il58| . though for 
an exact solution its choice is irrelevant. 

Thus, Instead of deriving the functional renormalization group equations from the quan- 
tum equations of motion, a better starting point is the effective action F where a regulator 
has been inserted. This then yields a starting equation, similar to f|T5l) for DSEs, as a 
starting point. The FRGEs can then be derived directly from this equation. 



8A 



6S_ 



d'^x 



6 



Tl{x,y) + Rk{x,y) 6A 



+ A 



(52) 



where R is the regulator function and k is the regulating scale, which has to be sent to 
zero to obtain the full result for the correlation functions. From this generating functional, 
similar as for the quantum equations of motion, a hierarchy of coupled non-linear integral 
equations for the vertex functions can be obtained. 

However, the structure of these equations is significantly different from the ones of the 
equations of motion. One obvious difference is the appearance of the regulator function. 
More importantly, the appearing loop diagrams are always only one loop. The resulting 
drawback is that only full vertex functions appear, instead of also bare ones as in the 



equations of motion [39 



The two systems of equations are thus suffi cient l y di fferent that finding a solution to 
both of them is a non-trivial consistency check 1591 . Il60| . It has been shown that indeed 



the one-parameter family of solutions, which is obtained from the DSEs by selecting the B 



arameter, can also be obtained with FRGEs |37l . Il61| , in particular the one for B = oo 



1571]. On the other hand, the fact that indeed both systems are solved by the same 



solution is sufficiently restrictive to show that th e infrared structure of the vertex functions 



has to be qualitatively unique in the B = oo case |159l . ll60| . under certain weak assumptions. 



Furthermore, the freedom in choosing the regulator can be used to compensate for the 



effects of truncations [39|, |158[. As a consequence, the results from functional studies for 
propagators most closel y resenab ling the ones from lattice calculations up to now have been 
obtained using FRGEs 371 . Il62[ |. Nonetheless, for the solution of functional renormalization 



group equations otherwise the same applies concerning the necessity of a truncation as in 
the case of the quantum equations of motion. 
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In total, both sets of equations have, from a practical point of view, their individual 
strengt hs and weakn esses. Exploiting the combination of both of them has so far been most 



useful |15lL Il59l . llGOj . and promises to be so in the future as well. Given that both systems 
of equations are very similar nevertheless, for simplicity here only the case of DSEs will be 
treated, with only given some additional results from FRGEs where they provide particular 
useful additional insights. 

3.2.3. Ultraviolet 

At large momenta, the leading behavior of the correlation functions can be obtained 
using perturbation theory in the coupling constant, due to asymptotic freedom [l, [3o|. The 
corresponding equations are obtained by expanding all correlation functions in the DSEs in 



the coupling constant [33[. As a consequence, the contributions from correlation functions 
of order n + m in equations for correlation functions of order n are then subleading by m 
powers of the coupling constant, and can thus be neglected if the expansion order is less than 
m. Therefore, at any given order of the perturbative expansion the number of equations is 
finite, and identical to the ones constructed with conventional Feynman rules. 

Since the results in two, three, and four dimensions are interesting and have turned out to 
be important to the understanding of Yang-Mills theory, it is worthwhile to also analyze the 
perturbative expansion in all cases. However, in two dimensions the perturbative integrals 
are ill-defined, and require an infrared cut-off, not leading to a viable perturbative expansion 
(t^P^ . Thus, only the case of three and four dimensions remain, which exhibit a qualitatively 
different behavior. Therefore, it is worthwhile to consider the situation in four and three 
dimensions separately. 

The characteristic properties of perturbation theory are already obtained from the leading 
corrections to tree-level. At this order, all vertices become their tree-level counter-parts 
in the propagator equations. Corrections to the vertices are then of higher order in the 
coupling constant. The much more complex case of higher-order corrections including higher- 
or der c orrelation functions is in principle analogous, though much more challenging in detail 



The generic properties of correlation functions at large momenta can already be obtained 
from a dimensional analysis. This can be illustrated most easily for the case of the prop- 
agators. For this analysis it is important that Yang-Mills theory is renormalizabl e in four 



dimensions, and in Landau gauge the highest appearing divergence is logarithmic |164 
As a consequence, in four dimensions, the generic propagator behaves as j2| 



p^D{p)-^C(l + u\n^] (53) 



where C is a normalization factor and Avm is a dynamically generated scale. The appearance 



of it, called dimensional transmutation [13|,|30||, is a pure quantum effect. The classical theory 



^^Strictly speaking, this also applies to three dimensions (l63i|, but there the breakdown of perturbation 
theory only occurs at higher order, and is thus not relevant for the purpose here. In particular, the Linde 
problem at finite temperature _23], to be discussed in chapter [5j is also connected to this. 
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is scale invariant, and thus the propagator would have to scale like the canonical dimension 
of the fields, and thus p^D{p) would be constant. The quantum effects break this scale (or 
dilatation) symmetry of the classical theory, and give the theory a scale, Aym- Since this is 
the only scale, p^D{p) can only depend on p^/AyM- The logarithmic behavior is then typical 
for such renormalizable field theories like Yang-Mills theory. 

In contrast, in three dimensions the propagators asymptotically take the form 

p'Dip) ^ 1 + ^ + O f 41 (54) 

with some constant a. This is the simplest, asymptotically free possibihty just by power- 
counting, since g"^ has dimension of mass. The classical theory is therefore not scale- invariant, 
in contrast to the four- dimensional case. 

Though, as stated, in the two-dimensional case perturbation theory is not applicable, 
the same arguments, and the fact that g now has the dimension of mass, suggest that the 
asymptotic behavior will be 

P'D(P) ^ 1 + ^ + O (3 . (55) 



which turns out to be approximately correct |90l . Il65j . It is thus similar to the behavior in 



three dimensions. It should be noted that the possibility for the expansions and fl55]) 
are solely based on a dimensional analysis, and therefore remain valid even if perturbation 
theory is not applicable. 

Higher order correlation functions show, of course, a more complicated structure, as the 
number of external momenta increases [2|. Still, the characteristic perturbative behavior is 
logarithmic in four dimensions, and polynomial corrections to the tree-level value in l ower 



dimensions. Note that in lower dimensions, operator-product expansion corrections [30|, |164 
could mix with the perturbative contributions, in contrast to four dimensions, where they 
are sub-leading. 

Three dimensions 

At leading order in the perturbative expansion the two-loop diagrams in the gluon equa- 
tion drop out, since they are of order g^, compared to the order g"^ one- loop diagrams. A 
second important point to this order is resummation, i. e. the inclusion of all one-particle 
reducible contributions for a self-consistent solution of the DSEs to a given order |2|, Il3, 33 



In four dimensions, resummed and ordinary perturbation theory does not coincide. How- 
ever, in three and two dimensions resummation only produces effects from order g'^ on. This 
follows from dimensional arguments alone, since only tree-level expressions can contribute 
at order g"^ in the loop diagrams, such that a polynomial correction of type g'^/p appears. 
The next higher orders are of g^, and thus require a factor However, there resum- 

mation sets in, as factors {g'^/pY appear in next-to- leading-order from the leading order. 
Therefore, resummation effects do not occur at leading order, but only from next-to-leading 
order onwards. 
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The next problem is that each of the loop integrals in the propagator equations is possibly 
linearly divergent from superficial power-counting, and only their combined result is finite. 
The latter can be deduced from the fact that Yang -Mills theory has only a logarithmic 
divergence in Landau gauge in four dimensions 164| . To cope with these divergences in 
perturbation theory it is convenient to use dimensional regularizatioijfl 0, 30 1 to perform 
the integrals. However, dimensional regularization requires an analytic continuation to non- 
integer dimensions. It is a-priori not clear whether such an analytic continuation exists for 
the non-perturbative correlation functions ^6j. Even if it exists, it is generally not known 
beforehand, and thus dime nsional regularization can only be applied if an assumption on 
this structure is made 167|. The alternative is to use a different regularization scheme. 

It is thus more instructive to use instead of dimensional regularization the same regular- 
ization scheme used later for the non-perturbative calculations. This is even more important 
as, besides the lattice, no non-perturbative regularization scheme is known, which preserves 
manifest gauge covariance 86|]. This is not in itself a p roblem, as this merely implies that 



the counter-terms will also not be gauge-covarianio 1641 ] . Only the renormalized correlation 
functions will then be again gauge-covariant. Thus, studying this effect already in pertur- 
bation theory, where it is possible to compare to a gauge-covariant regularization scheme, 
provides insights into how to treat these problems in the non-perturbative calculations. 
This does not matter for the g host equation, which due to the non-renormalization of 



the ghost-gluon verte^o |170l - ll72l | is actually fully finite in three d imen sions. Performing 



the loop integral yields the result for the ghost dressing function as |l66 | 



6' 



ab 



G{p) 



6 



ab 



16p 



(56) 



where Ca is the adjoint Casimir of the gauge group. Thus, the perturbative correction is 
polynomial in the momentum. Hence, at very large momenta just the tree-level behavior 
remains, in a manifest display of asymptotic freedom. On the other hand, at a small, finite 
momentum the dressing function diverges, a manifestation of the Landau pole, and the 
explicit breakdown of perturbation theory. 

The equation for the gluon propagator is different. The cutoff-regularization produces 
spurious divergences in this ca se. F irst, it is useful to note the result obtained using dimen- 



sional regularization, yielding |166 
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Which in odd dimensions imphcitly already fixes the renormalization scheme (3. Il64 flGGf . 
^•^Note that such a procedure modifies the corresponding Slavnov- Taylor identities in a well-defined way 



33,111,068,] 

^'^An explicit leading-order confirmation can be found in [16] 
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where in the first hne the first term is the tree-level contribution, the second term comes 
from the ghost loop, the third term from the gluon loop, and the tadpole vanishes upon di- 
mensional regularization. The individual diagrams are not individually transverse, but only 
their sum reproduces this property [30j. This cancellation of the longitudinal components is 



actually guaranteed by the gauge condition [40| and by the Slavnov- Taylor identities (STIs), 
as discussed in section 1X51 

This should be compared to the result obtained using a cut-off regulator. 



.abd'^CA [ . , P^Pu\ g^CA^.ab, 



64p \^ '^'^ P^ ) Qtx'^p^ 



32p \^ '^'^ P^ J GiT^p^ 
Stt'^P^ 

where the first line is the tree-level part, the second line comes from the ghost loop, the 
third from the gluon loop, the last from the tadpole, and A is the cut-off. The finite part 
agrees with the previous result. However, there remains a longitudinal part proportional 
to the cut-off, which is spurious, and behaves like a glu on rn ass. It emerges since a cut-off 



regulator does not respect the gauge symmetry [30|, |86|, Il64j |. As a consequence, the gluon 
propagator is no longer transverse, as it ought to be in Landau gauge. 

However, this is not a problem. Such spurious divergences or other artifacts of a non- 
gauge- inva riant regulator can always be compensated for by non-gauge-invariant counter- 
terms 



164j |. In this case, the counter-term Lagrangian to be added takes the form 



= Alix)Alix)6m' 



6m 



which has the dimension of a mass squared since the coupling constant squared in three 
dimensions has the dimension of a mass. This counter-term makes the gluon self-energy once 
more finite and transverse, and in agreement with the result using dimensional regularization. 
It should be noted that these counter-terms are fully fixed by the requirements of gauge 
invariance, as in the present case, and thus no additional arbitrary renormalization constants 
appear. 

Similarly, any other kind of such artifacts of a non-gauge invariant regularization pro- 
cedure can be absorbed in corresponding counter-terms. Thus, the one-loop perturbative 
result for the gluon dressing function becomes 
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This already shows that in the non-perturbative treatment, where dimensional regulariza- 
tion is not possible without knowledge of the analytic structure, such counter terms will 
naturally appear, and have to be taken into account. 



Four dimensions 



The situation in four dimensions is different for two reasons. One is that the dressing 
functions are not finite, but actually divergent, and therefore have to be renormalized even 
when using dimensional regularization. The second is that resummation already sets in at 
one-loop order. This is due to the fact that corrections are logarithmic instead of polynomial, 
and two logarithms at different order can be combined into one, as the momentum depen- 
dence is the same. Thus, to obtain a self-consistent solution of the DSEs at a given order 
requires to perform resummed perturbation theory in four dimensions. Since the DSEs are 
inherently self-consistent, only such a solution makes sense as a perturbative approximation 
to the full solution. The resummed solutions can be obtained with a variety of methods from 
ordinary perturbation theory [2I, l^oj, but here simply an appropriate ansatz will be made. 

The form of the perturbative dressing functions in four dimensions is 0] 
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(58) 
(59) 



with the anomalous dimensions S and 7, and the parameter u, which all three will be 
determined from self-consistency, the subtraction point s, and the values Z{s) and G{s) of 
the dressing functions at the subtraction point. 

The perturbative domain is reached for momenta large compared to the scale Aym- It is 
furthermore sufficient to only inspect the large momentum tail o f the inte grals, cutting off 



the loops in the infrared at p, the scale of the external momenta |173l - ll75 



Since the external momentum is small, it is justified to approximate the internal dressing 
functions hj G{p + q) ~ G{q) and then performing the angular integrals explicitly. The 
explicit form of the ghost equation is in this limit given by 
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with y = and the wave-function renormalization constant Z3, and the upper cut-off just 
indicates the use of a regulator. As in three dimensions, there are no spurious divergences, 
even when using a cut-off regulator. The (logarithmically) divergent upper boundary of the 
integral can then be absorbed by the wave-function renormalization constant. Putting in 
the ansatze (1551) and (|5^ yields that the solutions must fulfill 



UJ 
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to be self-consistent solutions, reproducing the known relations ^ between these perturba- 
tive parameters. 

For the gluon equation, a new problem appears. In fact, there is no self-consistent 
solution to it and t he ghost e quation simultaneously, as long as the three-gluon vertex retains 



its tree-level value |173l |175[. The reason is that, due to resummation, high-order effects of 
the vertices also contribute to generate a self-consistent solution. This does not affect the 
ghost equation, due to the non-renormalization of the ghost-gluon vertex. Therefore, it would 
be necessary to solve the corresponding equations of propagators and vertices simultaneously. 
This can be addressed better with the renormalizat ion-group improved perturbation theory, 
but this is a rather lengthy calculation ^ [s^]. Here, it is sufficient to note that this yields 
the missing quantity S as 

* = (61) 
In fact, the knowledge of (58|) and (l59i) can then be used in turn to limit the structure of 



the three-gluon vertex [176| . This is important for truncations to be discussed below. 

Again, asymptotic freedom yields that these are indeed the leading ultraviolet contribu- 
tions to the propagators, and at high-energies the non-perturbative contributions will only 



be sub-leading corrections to it, falling off as inverse powers of the momenta |30l . 

As a final remark it should be noted that the gauge algebra does not enter the ultravi- 
olet behavior in a qualitative manner, but only alters some of the coefficients, and in four 
dimensions the value of the anomalous exponent at higher orders, but without changing its 
sign. 

3.2.4- Infrared 

In analogy to finding the (resummed) self-consistent solutions at large momenta, the aim 
of a non-perturbative infrared analysis is to find self-consistent solutions to the equations 
at asymptotically small momenta. Two qualitatively different types of solutions have been 
found over the time, the so-called scaling and decoupling solution^. Here, the latter name 
is chosen for the characteristic momentum-behavior of its running coupling in a certain 



scheme, which is infrared vanishing. It is also referred to as a massive solution [15l] because 
the gluon exhibits a finite screening mass. However, as discussed below this occurs also in 
the scaling case with a finite or infinite screening mass, the later being an overscreening 
behavior. Thus, this distinction is not entirely satisfactory. The most elementary difference 
is found in the ghost dressing function, which is for the scaling case infrared divergent, while 
it is infrared finite, like for a photon, in the decoupling case. Thus it is possible to refer to 
the decoupling case as the finite-ghost case, which will be done here. 

For a discussion of the infrared behavior, it is quite useful to discuss finite-ghost-like and 
scaling-like behavior separately, which turn out to be quite distinct. Though the scaling 
solution will be formally obtained in the limit of S — > oo, it is not settled which is the 



Ori ginally, as a th ird possibility a singular solution has been proposed, embodying the concept of infrared 
slavery |l3 , 177 , 178| . This possibility is at the current time not a consistent solution of functional studies 
anymore, and there are no indications using lattice gauge theory that it could exist, see chapter 21 
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relation of the scaling case to this gauge-fixing condition, in particular if it turns out that 
only finite B values or even only a single value of B is admissible in the first Gribov region. 
This will be discussed further in section I4.3.2[ Until then, the results for the scaling solution 
will be just displayed alongside with those of the finite-ghost type, which are obtained for 
finite values of B. 

In any case, a momentum is infrared if it is very small compared to all other scales. At 
zero temperature the only scale of Ya ng-Mills theory is Aym, and thus an infrared momentum 



is a momentum satisfying [l45l . Il75 



p < Aym, 

while at finite temperature infrared momenta satisfy 

p < Aym, T, 



(62) 



(63) 



irrespective of the relative size of T and Aym- If Yang- Mills theory is coupled to matter, then 
the momenta also have to be much smaller than any mass or other dimensionful parameters 
introduced by this coupling. 

The quantification of much smaller than all other scales, and thus the value of the en- 
ergy scale beyond which an asymptotic behavior is expected, is a highly non-trivial question. 
Investigations of string-breaking in four dimensions indicate that the typical length scales 
where dynamical effects can still be enco unte red in Yang-Mills theory coupled to matter is 
of the order of 50 MeV, or about 10 fm 
is known as the A^-ality domain, i. e.. 



179 



relevant 180 



In the context of Wilson lines, this domain 
the domain where only asymptotic properties are 
Another indication is that the typical scale obtained in lattice calculation 
for Aym is 25-50 MeV. Both indicate a lower distance limit from where on pure asymptotic 
behavior could be expected. However, it is well possible that the behavior is already almost 
asymptotic at much smalle r distanc e s, bu t in terms of correlation functions rather the con- 



verse seems to be the case 115- 118. 145 



Scaling behavior 



One hallmark of a scaling-type behavior is that diagrams can be classified as leading or 



non- leading when all momenta are asymptotically small |l8l| . It is then actually possible 



to provide an analyti c sol ution to the whole tower of equations on the level of the sca ling 



In 



dimensions (l51 . 181, 182 1, which under certain weak assumptions is unique 159 , 160 
particular, for the propagator equations it can be shown that only those one- loop diagrams 



are relev ant wh ich c ontain a ghost loop, and that in general ghost dominance [65|, l9J, Il75 



183l4l85| holds |l8ll |. To illustrate the method of this so-called infrared analysis again the 



propagator equations will be useful. 



The relevant part of both equations then read [9J, Il75l . |186 
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and all contributions which turn out later to be self-consistently sub-leading have already 
been neglected. Comparing to figures ITT] and fT2| these are just all diagrams containing at 
least one ghost-line on the right-hand side. 

The only further input needed is the then the full ghost-gluon vertex Y^^^'°'^'^, or actually 
only its infrared-lea ding part. The consistency of the scaling solution actually demands it 



to be infrared finite 159 



160l | , though variations as a function of the angle between the two 



independent momenta are still permissible, but will be ignored for the sake of the argument 



Such possibilities have been explored and found to be only a quantitative effect |l86l . Il87 . 
Furthermore, the results shown in section 14.1.21 indicate that any such variations are rather 
small. Note that an additional color tensor proportional to the symmetric color tensor d"'^'^ 
would get lost due to the contraction with the antisymmetric tree-level vertex. Hence, at 
least the leading infrared part will not have color off-diagonal elements for either propagator. 
Thus, the ansatz is to replace the f ull ghost - gluon vertex b y its bare form. 



Making then the scaling ansatz [QJ, Il75l . Il83l . Il86l . II88 



2 



(65) 
(66) 



it is possible to absorb the remaining tree-level term in (IMllJgv implementing the boundary 
condition of an infrared divergent ghost dressing function [175|, i. e., setting 1/5 = 0. 

The integrals (J65l) and f l66|) can now be evaluated exactly, using the dimensional regu- 
larization formula 13011 
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which entails a definite minimal subtraction scheme in odd dimensions |l66l |. The essential 
point for the validity of this formula is that, because of the absence of other scales, the 
integral is dominated by the contributions around the external moment |175i] . Since this 
in tu rn depends on the singularity structure of the propagators, and in general of the vertices 



186| . this yields that the integration turns, up to a constant pre-factor, into amap from 
the topology of the diagram to the exponent of the external momentum scale 189 . 
The integrals can then be performed analytically to yield 
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g'CAAlAzIzif^^c, d)p-^'~'^'^^--. 



(68) 
(69) 



Attempt s to determine corr esponding general formulas for the case in the presence of further scales can 
be found in tl89l. see also II60I. 
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The expressions Iq and Iz [94, II86 
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are functions depending solely on the exponents Kgc and kaa, and the structure of the 
underlying space-time manifold, symbolized by the dependence on d. 

Counting powers of momenta in fl68tl69p yields the fundamental scaling relation 9J, Il83 . 



186 



= 2k,cc + I^AA + 



d-A 



(70) 



Thi s basic relation is the single most important property of the scaling-type solution jl59i 
160l . If it is not fulfilled, the solution will not be of scaling-type, irrespective of the behavior 
of the individual propagators or vertices, as long as the ghost-gluon vertex is infrared finite. 
The remaining consistency condition 



/z(Kgc, d) = Ig{hcc, i^aa{i^cc, d), d) 



(71) 



implies that the exponents depend only on the space-time manifold. The solution of this 
con sisten cy condition for a bare ghost-gluon vertex yields k^c ~ 0.595 in four dimensions 
94 Il86l |. K rr = 1/2 or Hcc ~ 0.39 in three dimensions, and Kcc = 1/5 or Kcc = in two 
dimensions [9^]. As a consequence of ( 170|) . the gluon exponent satisfies kaa < — 1, and the 
gluon propagator is infrared vanishing. A formal second solution Kcc = 1 of (17T|) in four 



dimensions is not admissible, as it cannot be brought into agreement with renormalizability 



IO4J . Il73| . If the ghost-gluon vertex has a n on-tr ivial angular dependence at zero momentum 



the values of k can change quantitatively |187 



actors A^gA 
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Furthermore only the product of the pre- 
is fixed by the equations, and their individual values cannot be determined 
Also this product will be affected by angular variations of the ghost-gluon vertex 



A direct consequence of the setup of this infrared analysis is that the exponent k^c must 
satisfy the inequality 

d-2 

n-cc > ^— , (72) 



otherwise the basic assumptions made, 
equation is the leading te rm, is invalid. 



i. e., that the ghost loop in the gluon propagator 
If the equality applies, indeed the ghost loop is as 
leading as all other loops [188[, and the gluon propagator becomes infrared const ant, and its 



therefore appearing screening mass receives contributions from all other loops |l9lj . This 
furthermore implies that the relation f lTT]) receives corrections on top of any corrections from 
the ghost-gluon vertex. This option has not yet been widely explored, though it may be 
relevant. Note that furthe rmor e arguments exist, which suggest a lower upper bound in four 
dimensions of less than 1 192 . 
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A much more stable consequence of the scahng solution t han the numer ical value of k. 
is that the renormalization-group-invariant running coupling 175l . Il83l . Il93 
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where A^i is the size of the adjoint representation, is constant at zero momentum, as this 



follows directly from fjTOj) . Of course, in dimensions smaller than four there is no dependency 
on the renormalization scale fi. It should be noted that, except for an overall scale factor of 
the adjoint Casimir in the ru nning coupling, the underlying gauge algebra has not entered 
into this result 1651 . Il83l . Il94j . Thus, the infrared in the scaling case should be qualitatively 



Furthermore, the relation 
the ghost-gluon vertex 



175 



independent of the gauge algebra. The only possible exception occurs, if again equality holds 
in f[7 2l). in which case contributions from the two-loop diagrams could, in principle, enter 

[i9ir 

fl73p also holds at finite momenta, and can be deduced from 
It is thus sufficient to know both propagators to calculate 
the full running coupling in Landau gauge. However, to ensure the correct renormalization 
behavior of the running coupling, in particular the relation 

Zi, (74) 

which follows from multiplicative renormaliza tion alone, requires that the renormalization 
conditions for both propagators are locked as |l75l | 
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Herein, Zg denotes the renormalization constant of the gauge coupling, Z^ the wave-function 
renormalization of the gluon, and it has entered that in the renormalization scheme em- 
ployed, see section 13. 4^ the renormalization constant of the ghost-gluon vertex Zi has been 
set to one. Of course, an alternative is to introduce a scheme with an explicit (finite) renor- 
malization constant Zi of the gh ost-gluon vertex, which will then appear as a factor on the 



right-hand side of (JZH) [103|, UM 



In dimensions smaller than four, the simplest possibility to ensure (175|) is to formally set 
II to infinity, in which case Z(oo) = G( oc) = 1, just because of asymptotic freedom, as has 



been discussed in section 13.2.31 l8ll . Il66 



Based on the insight gained w i th t he propagator equations, it is possible to turn the 



results into a general rule |l59l . II6OI . Il89| . akin to the power-counting ru les u sed to determine 



the superficial degree of divergence in the renormalization process |164| . In particular, 
it is possible to calculate the superficial exponent 5 of a diagram, which determines the 
dependency on the single mom entum scale a t the symmetric point, i. e. the momentum 
configuration pf = . 



151, 159, 181 



To obtain this result, the first step is to calculate the scaling behavior of a diagram in the 
DSEs (or FRGs) at this symmetric point in the infrared. Similar to corresponding formulas 
in the renormalization program, a formula can be constructed which maps such a diagram 



onto this exponent. For a diagram having n^. external legs of field typq_] (pi it is given by 



For the case of mixed propagators, see [7J|. 
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dressed vertices r 



where rirp^ is the number of external propagators with their corresponding exponents k^-,^. 
signifying their power-law behavior, and the n^.^ jf,.^ count the number of bare and full 
vertices, respectively, appearing in the diagrams having infrared exponents t^^^^...^^^ and 
canonical exponents c,^^^...^^^, respectively. Finally, c is the overall canonical dimension of 
the diagram, and / is the number of loops. Similar to the case of the superficial degree 
of divergence 164| it is also possible to recast this expression using graph theory into one 
involving only the critical exponents k and the type of internal lines, but no longer the 



canonical dimensions 151 



For the purpose here, it is sufficient that such a mapping exists. 
With this mapping at hand, and the assumption that the superficial infrared exponent is 
indeed the correct one, i t is p ossible to determine a recurrence relation between the infrared 



exponents |l5ll . ll59l . ll60l . ll8l[ |. It follows from the fact that the exponent of the left-hand-side 
of a DSE must be the same as the dominating, i. e. most singular, one on the right-hand-side. 
The simplest example for this is the gluon equation keeping only the tree-level term and the 
ghost-loop. The corresponding statement is then for a bare ghost-gluon vertex 



Kaa = min ( 



(77) 



Since the ghost equation implies that Kgc > |l88j |. the second option is the dominating 
one, and the tree-level term is subleading. In analogy, such a hierarchy of relations can be 
obtained for all equations. To actually resolve them requires a further piece. This comes 



from comparing relations of type fl77p for the same quantities in DSEs and FRGs 159|. Since 



in FRGs all vertices are always dressed, while in the DSEs always one vertex is undressed, 
both hierarchies are not identical. This is sufficient to find an explicit resolution of the 



infinite hierarchy of inequalities |160 



This result yields the infrared critical behavior of all correlation functions at the sym- 



metric point for Yang- Mills theory in the scaling case as |181l . 1182 



p{n,m)^^^2^| ^ ^^p2^J(m-n)K^;e+(n-l)(f-2)^ 



(7^ 



where n is the number of external ghost legs (identical to the number of external anti-ghost 
legs) and m is the number of external gluon legs. An important consequen ce of this result 
is that any coupling constant defined from the higher vertices is also finite 18lL Il82| , by a 
construction similar to f l75]) . 

Note that the construction presented here is given for Yang-Mills th eory only. Vari- 



ous investigations have been performed when matter fields are included |l95l4l98| . but it 
is possible that some assumptions necessary for the construction then break down in the 

57 



matter sector |l60l ]. This has not yet been finally clarified. There are also investigations 
how to extend this analysis to nor i-sym metric momentum configurations 189| , though these 
encounter in general ambiguities |l60j |. 



Finite-ghost behavior 



The finite-ghost case was first di scussed in 
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16, 37, 64, 69, 70, 103-106, 120 



199 



200| and has a fterwards been investigated 



1891 . I20ll - l21l| . The finite-ghost case is 



m, e. g. 

characterized by an infrared finite gluon propagator and an infrared massless ghost, 
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In this case, it is more complicated to make a general infrared analysis as in the scaling 
case. Taking the ghost equation, and choosing a finite-ghost behavior by setting 5 to a 
finite value, the tree- level term is no longer canceled in flMl) . Inserting then the finite- 
ghost ansatze in a truncation at the same level as in the scaling case automatical ly y ields 



that the ghost is infrared massless and the gluon is infrared screened [15|, ll6|, l37|, Il04l | . In 



particular, in the infrared limit the appearing mass scale acts as a bare mass, and only finite 
momentum corrections make it suf ficien tly fast vanishing in the ultraviolet to ensure the 
correct ultraviolet behavior 15|, |37|, 1202 . 



The reason for this is rather simple. Investigating the appeariri g integrals in detail show 

even if not introduced 
Beyond 



that the diagrams in general harbor a mass-like term |166l . |190 



by assumption [15j. In perturbation theory, these terms cancel exactly [30|, l212 



perturbation theory, thi s is in general no longer the case. However, in the scaling case, they 
are infrared sub-leading 1661 . Il74j . ll90[ |. since they behave like m? /p^, where m is some scale, 
generically proportional to Aym- The contribution due to the ghost-loop for a scaling ghost 
and tvcc > {d — 2)/4:, e. g. in the gluon propagator DSE (jSSD, provides a stronger infrared 
divergence than just j9~^, therefore dominating the screening behavior and leading to the 
scaling behavior. This option is, by construction, switched off in the finite-ghost case, and 
thus the screening behavior becomes leading. This does not occur for the ghost due to the 
factorization of the external ghost momentum of the ghost-gluon vertex. This makes terms 
automatically more infrared divergent, promoting the screening mass to a correction of the 



wave- function renormalization [145[. Note th at th erefore at least all one-loop contributions 
in the gluon equation contribute equally [37|, 1213 . 

Thus, a finite-ghost behavior emerges naturally, as soon as any other infrared-dominating 
behavior is switched off by choosing B finitj^. Similar effects are expected to also occur in 



the vertex equations, closing the system self-consistently |214| . However, this has not yet 
been investigated in detail on the l evel of a self-consistent solution of the vertex equations, 
but the investigations are ongoing 215 . 



^^This can also be achieved, e. g., by an explicit infrared regulator like a finite volume, see [14 
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On the other hand, this mechanism imphes the possible existen ce of a scahng window in 



the correlation functions' behavior even for the finite-ghost case [37|, Il45| : Define the momen 



turn pb to be the one at which the normalized ghost dressing function p\Dg{j)%) / {j^^ Dcin'^)) 
becomes similar of size as 1/B. Then screening becomes manifest at momenta at or below 
Pb- If -B is chosen such that p_B is much smaller than the scale Ay^, then there exists a win- 
dow where the correlation functions show a scaling behavior, up to sub-leading corrections, 
in the momentum range 

< p < Aym- (79) 

Since the power-laws exhibited in this energy regime are the same as expected for the scaling 
case, this will be called the scaling window. The existence of this window follows from the 
fact that for momenta within the range (jTHD no sensitivity to the deep infrared yet exists. 

The one technical challenge associated with the finite-ghost solution is that any such 
screening mass inevitably mixes with the spurious divergences discussed in the perturbative 
case when a non-gauge-invariant cut-off regularization is used. Therefore, disentangling the 
different contributions, in princi ple, re quire s in this case a resolution of the corresponding 
broken Slavnov- Taylor identities |37l. l39l . ll68j . to obtain the screening mass, e. g., of the gluon 
propagator. In general, this is complicated, and instead usually the necessary counterterms 



are mod e led s uch that corresponding lattice and perturbative results are reproduced [37 
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or an appropriate truncation is chosen to eliminate these artifacts 
Finally, as is clear from the fact that all free parameters of Yang-Mills 
theory are fixed once Aym is fixed, the gluon screening mass is not an independent parameter 
of the theory. Furthermore, the perturbative cancellations enforce that the screening terms 
have to drop faster than l/p^ at large momenta, thus recovering the masslessness of the 
gluon in the high-energy limit. 

Note that the effective wave-function renormalization of the ghost is less arbitrary. Since 
its value is uniquely fixed to the B parameter, no problem would emerge even if spurious 
divergences were encountered in the ghost equations. 

3.2.5. Intermediate momenta 

At the current time, in order to obtain solutions at all momenta, it is necessary to 
truncate the system of equationfl In general, the main ingredient is to drop all equations 
for higher order correlation functions from a certain, usually low, order n. This requires to 
make assumptions for the various higher n-point correlation functions still appearing inside 
the equations, and which are no longer determined self-consistently. Various schemes have 
been employed for this purp ose. Such schemes include mo deling the remaining vertices along 
some guiding principles |37|, Il57l Il66l . 1 1731 . ll75L Il95l . 1215 , truncations based on the desired 



concepts |13l . l214j |. lattice input [187|, and others [13|, llJ, |39| . The possibilities are essentially 



properties of the solutions inferred fr om other sourc es isl . 104 . 105 , 218 |. methods based 
on a bac kground- field formulati on 15l . Il22l . |202| . |219|. l220 l. effective actions 7^, analyticity 



only limited by the necessities imposed to answer a particular question. 



'Very few systems, like the l+l-dimensional Schwinger model, are known where this is not necessary 
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As an example, the truncation, with shght variations, to be used below in the chapters H] 
and [5] will be discussed. Alongside the guiding principles will be presented which motivate 
this truncation. 

The aim in the present case will be results for the propagators. This requires their 
respective equations, which are then truncated at one-loop level. Thereby models for the 
three-point functions are required. 

The first cent r al require ment is to establish contact with resummed perturbation the- 



ory 



37, 166, 173, 175, 216 



The main motivation to make this the central guide-line for 
the truncation is that this is the fastest way to make contact to experimentally accessible 
quantities, instead of using non-perturbative quantities like bound-states as a guiding prin- 
ciple. The second advantage is that this makes the non-perturbative approach here a direct 
extension of the perturbative case, guaranteeing a smooth transition from one to the other. 

The goal is to have agreement to perturbation theory at one-loop order. A consistent 
truncation scheme therefore requires to keep the equations for the propagators and include 
all one-loop diagrams, but to drop all two-loop diagrams and all higher-order equations. 
This level of truncation is the currently most commonly usecQ. Thus, it is necessary to 
specify the ghost-gluon vertex and the three-gluon vertex. 

One possibility to achieve this is to leave the ghost-gluon vertex bare, which even turns 
out to be a rather good approximation of its full behavior, see section 14.1.21 This is not 
the case for the three-gluon vertex. As discussed in section 13.2.31 reproducing leading- 
order perturbation theory already requires a modification such as to provide the correct 
self-consistent behavior for the gluon propagator in four dimensions to order . A po ssible 



asymptotic dressing realizing such a behavior for the three-gluon vertex is, e. g., [216 



az 
a 
^(ti) 



^0) 



2KAAil + 35) 



-2 - 65 + aa^ 

1 + 25 
-9 - 195 



255^1 



with k"^ = {p + q)'^, r}(t,p the tree- level vertex, and Zi the renormalization constant of the 
three-gluon vertex, which is introduced to ensure the correct renormalization of the vertex. 
This ansatz explicitly preserves the Bose symmetry of the vertex. This is the construction 
which will be used in section O for finite temperatures. The choice of the constant az is 
uniquely fix ed by requiring the correct renormalization properties and resummed perturba- 



tion theory [216[. The parameter ac is then free and chosen such as to make the vertex 



See for first progress in going beyond this level 195l l217j . 
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freeze in the infrared. This estabhshes a truncation scheme which is self-consistent in the 
perturbative domain. Note that in three dimensions the choice F"^ = 1 is sufficient to 
ensure the correct perturbative behavior. The only remaining problem is the presence of 
spurious quadratic divergences due to the cuto f f reg ularization. For this truncation, these 
can be de alt with by explicit counter-terms 2161 . l22l| , or by explicit and implicit subtraction 



166, 221 



A bit more elegant is the possibility to provide an additional dressing of both the ghost- 
gluon and the three-gluon vertex [37|. In case of the three-gluon vertex, it is given b} 



Vf^t^iP^q^k) = r'"'rfj(p,q,k)r-\q,p) 
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(81) 

fuviP,q] k) 



fuv{p,q;k) 



q^ 



g2 + A2 



uv 



where p is the external momentum and q is the loop momentum. This form completely 
eliminates the spurious quadratic divergences. The formfactors fi with their corresponding 
parameters ensure that at large momenta only the ultraviolet completion F"^ contributes to 
ensure the correct perturbative behavior. This completion can be chosen, e. g., as f lHOj) or 



as 



871. 1178 



F-(g,p) 



1 [G{p' + Aj^)G{{p + q)' + Aj^)] 



{l-a/5-2a) 



(82) 



Z, [^(p^ + A?«)^((p + g)2 + A?J](i+-) ' 
with a free, but conventionally chosen as 35, and Zi the three-gluon vertex renormalization 
constant. This choice breaks Bose symmetry, and also does not reproduce the results for 
the three-gluon vertex obtained using lattice gauge theory, shown in section 14.1.21 below. 
However, the difference at the level of the propagators is very small, in particular in the 
scaling case [37 . 

This choice also requires a modified ghost-gluon vertex to fully cancel the spurious 
quadratic divergences. Furthermore, in the scaling case the only contribution in the far 
infrared are those containing ghosts, and thus the transversality of the gluon propagator 
must be guaranteed by the ghost loop alone. This is not necessary in the finite-ghost case. 
This is ensured by the construction |37| 



r"" (p, q,P + q) = ip,^A{p, q,k)-i{p + q)f^B{p, q, k) 



^3) 



A{p, q, k) 
B{p,q, k) 

fmip, q; k) 
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1 -^fuv{p,q] k) 
^ fmip, q; k) 



(p2 + A|^)(g2 + A2^)(P + A 
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'^^Note that fuv is not the same as the one in [37[, and it is more generally applicable. I am grateful to 
Markus Huber for pointing out a potential problem with the original version. 
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Ghost-gluon vertex ansatz at the symmetric point || Three-gluon vertex ansatz at the symmetric point | 
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Figure 14: The two truncation ansatze for the ghost-ghion vertex (left panel) and the three-gluon vertex 
(right panel) in four dimensions. The functions A and B for the ghost-gluon vertex are given in (I83p . the 
function F"^ in ([5^ , the improved vertex in (ISTl) , and the Bose-symmetric one in (150)) . The pa rame ter Aju 
is 300 MeV, and Ajjy 20 GeV, the dressing functions are for the scaling case, using the fits in [l73l |. 



with the same form factor fuv as for the three-gluon vertex. Note that, in principle, the lon- 
gitudinal part should be determined by the STIs, but generically d epen ds on the four-point 
ghost-gluon scattering kernel, unknown at this level of truncation 175l |. Making an ansatz 
to obtain transversality in the scaling case is therefore the only possibility here. However, 
as emphasized, this longitudinal contribution is completely irrelevant when calculating the 
non-amputated correlation functions, and in particular physical observables. It should be 
noted that the transversality of the gluon propagator is guaranteed by this ansatz only at 
asymptotically large and small momenta. However, the residual longitudinal component at 
intermediate momenta, being a truncation artifact, is negligible 37|. Truncations which en- 
sure transversality at all momenta have been proposed in 15, 215l , though at the expense of 
disagreement with perturbation theory at large momenta. A truncation having both virtues 
is in principle possible to construct, though great care is necessary to avoid the introduction 



of infrared- ultraviolet mixing [174| . Il90| . Such a truncation would make use of the not yet 
included longitudinal tensor structures of the three-gluon vertex, which could be chosen to 
ensure transversality at all momenta. Since for the present case the effects appear negligible, 
this will not be done here. How well transversality can be achieved once the three-point 
vertices are included self-consistently has not yet been explored. However, for a final con- 
struction the four-point vertices and the two-loop terms in the gluon equation could play a 
role. 

These various ansatze for the vertices at the symmetric point are shown in figure [HI to 
provide a comparison to the results from lattice gauge theory presented in section l4.1.2[ It is 
visible that the vertices which have the ultraviolet cancellations of the quadratic divergences 
built in drop off sharply at large momenta, since there the divergent contribution from the 
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integral kernels are needed to generate the correct ultraviolet running. The Bose-symmetric 
ansatz, which requires explicit counter-terms, shows directly the logarithmic running, which 
is also encoded in the ultraviolet dressing part of the improved three-gluon vertex fl82i) . Note 
that this logarithmic rise may be altered if the ansatz were enlarged to include the other 
three p ossib le transverse (and ten more longitudinal) tens or structures o f the three-gluon 



14J], for which first studies have been performed in |154j . Il95l . |215 



vertex 

On the other hand, compared to the results shown in section 14.1.21 the bare ghost-gluon 
vertex, depicted also in figure [TH shows at the symmetric point a rather decent resemblance 
of the re sults obtained below in l attice calculations, and is in agreement with self-consistency 
studies |l53l . 1 1591 . II6OI . Il86l . Il89| . but the three-gluon vertex fails this test. Nonetheless, at 



the level of propagators this turns out to be of minor importance, in both the scaling and 



finite-ghost case [37||. Therefore, this departure from the presumably better results of lattice 
calculations seems to be of minor importance at the level of propagators. For other quantities 
this may not be the case, and for a systematic improvement of the truncation three-point 
vertices would have to be the first to be treated, see e. g. |195| for such an investigation. 
These are only two examples of em ployed vertex constructions. V arious alternatives and 



modifications have been explored |l5l . II66I . Il75l . Il86l . llSTl . l22lL |222| . The results turn out 
to depend only quantitatively on the employed version of the vertices, after fixing B. The 
only major exception is finite temperature, and there only for one polarization of the gluon 



216| . This will be discussed in more detail in chapter |3 Reiterating, the guiding prin ciple 



here has been to obtain the correct resummed perturbation theory, renormalization [37| , and 



for the scaling case infrared self-consistency. Conceptually different guiding principles have 
also beenemployed, in particular those guided by reproducing a particular infrare d behavio r 
il,lio3,lio5|. Going beyond such modeling of the vertices is an ongoing endeavor |l95l . 12231 



It should be noted that using FRGs instead of DSEs, self-optimizing flows and regu lators 



can be used to at least partly capture the effects of four-point functions j37|, Il62l l223j . This 



yields the currently best quantitative agreement between functional and lattice methods on 



the level of propagators [37 



3.2.6. Numerical methods 

This fully fixes the equations to be solved here. However, this leaves a coupled set of 
non-linear integral equations. Solving them is still a non-trivial numerical problem. Various 
methods exist to s olve these e quati ons, mainly based on iterative methods |200l . l224| . as well 
as local 174j . |225| and global 226| Newton algorithms. 



In solving the DSEs, there are two main challenges. One is the numerical integration of 
the loop integrals. The second is finding a solution for the self-consistency equations. 

Concerning the loop integrals, it turns out that these are somewhat non-trivial numeri- 
cally, especially for the scaling solution. The main problem is that the relevant integration 
range features three distinct problems. One is that a very large range of momenta is con- 
tributing to the integrals, when it comes to a quantitative precise determination of the 
pre-factors of the power-laws in the ansatze f l68ll69|) . Thus, an exponentially large momen- 
tum range has to be covered. The second is that the integral kernels feature singularities in 
the external momenta. In fact, this is the original reason why the infrared analysis discussed 
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before is working. The third is the necessary cancellation of divergences, which are of two 
kinds. The first is the subtraction of the spurious and physical logarithmic divergences. The 
second kind are phantom divergences, i. e., divergences which appear when the integral in 
the absolute value of the loop momenta are made, but are exactly canceled when the integral 
over the relative angle between external and internal momenta are performed. 

To deal with all of these problems, the following procedures have been found useful 



225| . Take the case where divergences are subtracted by a counter-term. The integral 
kernels can be separated into a part / depending only on the loop momenta's size g, one 
part h depending only the internal momenta and the external momenta p, and one part g 
depending also on the relative angle 6. A suitable rewriting of the integral is then given by 

n(p)-n(s) = jdqj defip,q)h{q)g{p,q,9)- j dq j dO f{s,q)h{q)g{s,q,e) 
= J dqh{q) (^f{p, q) j dOgip, q, 6) - f{s, q) j dOgis, q, 9)^ . 

In the case where / should be one, which actually sometimes occurs, the subtraction of 
the two g functions should furthermore be done before the integration. Note that this 
combination inside the integral is formally p art o f the regularization process [l64l ]. and can 



be thought of as related to a BPHZ scheme |221 



Furthermore, the integral in q can never be extended to real infinity in a numerical 
treatment, and it is often also not possible to evaluate it at exactly zero momentum. Thus, 
in practice a numerical upper cutoff A and lower cutoff e are included. To actual perform 
the numerical integration, the q integral is then split into 

^ I + I (84) 
42I 



to deal with the singularities of the integral kernel I. It is then only required to choose 



the integrator. A useful choice |225[ | is a Gauss-Legendre integrator, which samples the 
boundaries of the integration domain well, and thus the singularities appearing at the ex- 
ternal momenta. To deal with the exponentially distributed strength of the integral, it is 
furthermore useful to spread the points instead of on a linear scale o n a lo garithmic scale. 



This type of integrator, with probably adapted domains of integration |227l | , has so far been 
sufficient to deal with the numerical integrals encountered in such studies. 

It remains to solve the self-consistency problem for the sought for function D, e. g. 
symbolically for a single equation given by 

^ \ + U{p,D), (85) 



Dip) 



neglecting counter terms. To deal with the fact that D are functions, two approaches have 
been found useful. One is approximating the logarithm of the function by an expansion up 



^^For vertices |153l | or in lower dimensions [187| additional singularities may occur also in the angular 
integral, requiring a similar splitting there. 
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to a fixed order in an Chebychev basis 174 . l225l . |226| . The second is to to calculate the 



functio n only at a fix ed number of points and interpolate between them, e. g. with cubic 



splines |l87l . l200l . |224| . It has been found useful |226j to split the sought for function into 
a function which embodies some known properties, like the perturbative behavior or the 
qualitative infrared behavior, and another function multiplying this fit function. By this, 
only a function of order one has to be calculated, which stabilizes the numerics. 

To find the self-consistent solution with such an approximation an iterative p rocedure is 
used. In the case of a function basis, a Newton method can be used 
self-consistent solutions. In the second case a fixe d-point iteration 



20C 



224 



225 to find the 



is appropriate. 

In both cases, adaptive global versions 1871 . 12261 ] are useful to stabilize the iterations if no 
good starting guess is available, and the radius of contraction of the iteration procedure is 
not yet known. 

However, there are two drawbacks with this approach. One is that such an approximation 
is never exact, and in particular the results are usually only continuous, but not continuous 
differentiable over the whole momentum range. Thus, the final accuracy which can be 
obtained in solving the equatio ns is necessarily limited, though given sufficient computing 
power can still be rather good 226(1 . Secondly, in both cases the solutions are only known 
in a certain interval, and extrapolations outside this interval are necessary when performing 
the integrals, which due to the sum of internal and external momenta will necessarily also 
probe outside this domain. At momenta where perturbation theory is applicable, this is not 
a problem as the unfitted part can then be substituted by perturbation theory. However, 
at other momenta this is a pro blem, as extrapolation is usually not very stable. Here it 
has been found useful 224j . |225| to employ analytical knowledge, like the scaling behavior 
fl68tl69p . with free coefficients to extrapolate. The free coefficients are then chosen such as to 
guarantee continuity at the point where the numerical so lutio n ends, and are then adapted 



during the iteration process, possibly including damping [226 



If a system of equations should be solved, it is either possible to iterate all equations 
simultaneously or one after the other. However, the numerical stability may differ, depending 
on the details of the system. 

With these methods, it was so far always possible to solve this type of DSEs. However, 
one should be wary as these integral equations are rather forgiving concerning erroneous 
coding or inconsistent truncations. Secondly, it has rather often be encountered that a 
truncation lead to inconsistencies which surfaced by a non- convergence of the iteration pro- 



cess. A particular example [166| is e. g. that if the gluon loop in the equation for the gluon 
propagator is not sufficiently suppressed at mid-momenta due to a not appropriately chosen 
three-gluon vertex, it dominates and tries to lead to a sign change in the gluon propagator. 
When using an expansion of the logarithm of the propagators, this is impossible, leading to 
a failure of convergence. 

Last, but not least, most of these steps are rather similar, and can thus be automatized 



228l | , though still great care is required when dealing with both the physical and numerical 



particulars of a given system, as discussed here. 
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3.2.7. Dyson-Schwinger equations on a lattice 

One of the main advantages of DSEs and FRGs is that the can treat all momentum 
scales equally well. Thus, a natural possibility is to use them to extrapolate lattice results 
beyond the limited range of momenta accessible in lattice calculations fl40p . It is in particular 
interesting to determine the results of the functional equations at larger and larger volume, 
as well as finer discretizations, to ensure that for the range accessible to lattice calculations 
the correct behavior is reproduced. 

In practical calculation, the finite volume is usually intro duced as a sphere enclosing a dis- 
cretized hyper-cube, to simplify calculations 1451 . \224l |229| . Furthermore, the discretization 
is performed in momentum space, rather than in position space, to take the formulation 
of functional equations into account. Thus, the volume-dependence and discretization- 
dependence cannot be expected to be quantitatively the same as in lattice calculations. 
Nonetheless, good agreement at finite volume and discretization, as well as an understand- 
ing of the qualitative behavior, make the functional equations then a useful tool to perform 
these extrapolations. This type of extrapolations becomes more important when investigat- 



ing matter, in particular chiral fermions |230 



So far, this has only been performed to investigate the volume-dependence of the scaling 
case to predict at which volumes a scaling behavior can be obse rved instead of the finite-ghost 
case which will appear necessarily on any finite volume 145| . Such results are illustrated 
in figure [T3 This gave also rise to the notion of effective infrared exponents, which are 
det ermin ed under the assumption of a scaling window f l7^ . and taken as function of volumes 



90L Il45| . Investigations in two dimensions indicate that especially for the gluon propagator 



the volume dependence is rather well reproduced [9 

3.3. Slavnov- Taylor identities 

Though fixing the gauge is a difficult problem non-perturbatively, it also has an advan- 
tage. In particular, it restricts the structure of the vertex functions. A particular example 
of such a restriction is that the longitudinal part of the gluon propagator has to vanish in 
Landau gauge. This follows immediately from 40 



Pf,PuD^^ = p^p^ < A^A^ >= FT < {d^A^){d^A^) >= 0, 



^6) 



where FT denotes Fourier transform, and the Landau gauge condition in Fourier space, 
p/i^a _ Yisis been used. Another example is that the part of the non-amputated ghost- 
gluon vertex longitudinal in the gluon momentum has to vanish upon contraction with the 
gluon momentum 



= p^,< A^{p)c{q)c{p - q) >=< Pt,A^{p)c{q)c{p - q) > . 



17) 



Similarly, it can be shown that in the limit of g — )■ 0, t he vertex self-energy is not divergent, 
yielding a finit^ ghost-gluon vertex for all momenta 170 . 171 . 



^■^Finite in the sense that the renormalization constant is finite. 
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Gluon propagator in a finite volume from DSEs ] Ghost dressing function in a finite volume from DSEs ^ 




Figure 15: The results for the gluon propagator (left panels) and ghost dressing function (right panels) as 
a fun ction of volume (top panels) and of discretization at fixed volume of (13.5 fm)^ (bottom panels), from 
[l45| . This is only to illustrate the possibilities of DSEs. A detailed discussion of the solution properties is 
given below in section 2] 



These results follow from the constraints imposed by the gauge fixing on the gauge 
freedom. They can be generalized, yielding the Slavnov- Taylor identities (STIs) jH. In fact, 
these can be formulated just based on the presence of the gauge-symmetry of the original 
Lagrangian. After fixing a gauge, they yield a hierarchy of algebraic condi tions on the vertex 



functions, which can be obtained from the Zinn- Justin equations [39|, l86|, [16 



5T ST ST ST d.A^^ST 

+ 1 T-^T— = cut-on terms 



SK'^SAf, SL'^Sc'' 3 SU" 
/5r _ _ST_ 
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Herein K and L are sources for the gluon and ghost fields in the effective action F, re- 
spectively. The cut-off terms on the right-hand side represent contributions which are non- 
generic, and depend on the regulator used. In case of gauge-symmetry preserving regulators, 
they will be zero, but non-zero for a cut-off regulator. However, similar to the functional 
equations, these are a coupled infinite hierarchy of conditional equations. Only in perturba- 
tion theory the hierarchy, as for the functional equations, decouples ji], [ssj . Therefore, it is 
generally not possible to decouple these additional conditions. This can be illustrated [37 



by the STI for the three-gluon vertex, which is given by |144 



< Alc'd^c^ > - < A^^c'd^c' > 



-g {r'^' < AIAYc' > +f'^' < AlAld^c' >) . (88) 

Here, terms from the cut-off regulator have been dropped. The last term is evidently always 
of higher order in the coupling constant than the other terms, and thus does not contribute in 
perturbation theory. However, in a non-perturbative setting it cannot be ignored in general. 
So far, however, all no n-pe rturbative inv estigations of this STI have kept this term either 



at best approximately |231| or not at all |103L llOSl . Il44l . Il75| , since as a scattering kernel its 



structure is enormously comphcated. Furthe rmo re, it appears to be consistent to truncate 



the STIs at the same level as the DSEs [37|, Il75| . Thus, for the here employed truncation 



scheme, it would be self-consistent to drop it. 

The question naturally arises whether these conditions further restrict the solution man- 
ifold for the functional equations. In general, the answer is affirmative Q . However, Landau 
gauge is special in this respect0 (3^]- 

As already indicated by fl86|) . f l87|) . and flHHj) . the hierarchy of STIs involve the con- 
tributions longitudinal with respect to the gluon momenta of vertex functions. On the 
other hand, the non-amputated correlation functions have all gluon legs contracted with a 
transverse projector from the gluon propagators. Hence, no longitudinal tensor structures 
contribute to these. Since the DSE hierarchy of equations for the non-amputated correlation 
functions is closed, the longitudinal parts of the correlation functions cannot contribut^. 
The STIs, on the other hand, involve both the transverse and the longitudinal contributions. 
Since the transverse contributions are completely determined, the STIs are then constraint 
equations for the longitudinal contributions as a function of the transverse contributions 



only [371, Il62j . In particular, the STIs can be solved for the longitudinal tensor structures 
as functionals of the transverse tensor structures. 

This is to be expected: The STIs are relations manifesting the constraints due to the 
gauge symmetry. As such, they are conditions which are automatically fulfilled by the correct 
solutions, and thus do not contain any independent information compared to the hierarchy 
of DSEs. In particular, when determining physical observables, all correlation function are 



It is possible that very severe infrared divergences, proportional to d^A^, could appear in the vertices 
when taking the Landau-gauge limit, which spoil this argument. That this is not the case is a rather weak 
regularity assumption, and there are no indication that this occurs [37l.[l6^. 

■^^On the lattice, this is manifest, since only full non-amputated correlation functions are accessible. 
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fully contracted with transverse projector^, and hence the longitudinal tensor structures do 
not carry any physical information. Thus, the STIs cannot constrain physical observations. 
However, when attempting to solve the functional equations for the amputated correlation 
functions this will only be possible in general when, at least at the same level of truncation, 
the longitudinal tensor structures will fulfill the corresponding STIs. Otherwise additional 
truncation artifacts may be introduced. 

Two remarks must be made here. First, the hierarchy of equations for the amputated 
correlation functions requires the determination of both transverse and longitudinal tensor 
structures, and for the latter the STIs provide additional information. Second, any kind 
of truncation potentially violates the hierarchy of DSEs at one point, and therefore also 
the STIs in general are not fulfilled anymore, even at the same order. From these, the 
only available information which can be gained is on properties of truncations, but not on 
the correct result. An example is perturbation theory, in which case the violations of the 
STIs are well present, though always of one order higher in the coupling constant than the 
employed order of perturbation theory. Thus, it is possible to show that the violations occur 
because the expansion is truncated, but it is not possible to infer, e. g., the size of higher 
order contributions just from these violations. 



3.4- Renormalization 

An important technical issue to be returned to is that besides the spurious divergences 
discussed earlier there are the divergences which are physical, i. e., they result from the fact 
that Yang- Mills theory (at least in four dimension s) can only be a low-energy effective theory, 
but not a consistent theory of its own This can be directly seen in pert urba tion 

theory [2|], but lattice calculations also confirm this beyond perturbation theory [ol. Il37 . 

This fact can be hidden in the renormalization process, which permits to encode this 
unknown physics in the renormalization constants [2|. In the perturbative expansion of a 
renormalizable theory, like Yang-Mills theory, it can be shown that this is possible with 
a finite number of independent renormalization constants. For covariant gauges, like the 
Landau gauge, this is equivalent to multiplying correlation functions by appropriate chosen 
renormalization factoro, at least as long as no matter fields are involved. If this is also 
possible beyond perturbation theory has not yet been proven, though no evidence to the 
contrary exists. Since due to asymptotic freedom the process of renormalization, once mul- 
tiplicative renormalizability is assumed, can be performed essentially in the same way as in 
perturbation theory, it appears likely that this is correct. 

As a consequence, for t he purp ose s her e the renormalization can essentially be performed 



as in perturbation theory [13|, [IJ, l37|, Il75| , and will not be discussed in detail other than to 



state the renormaliz ation conditions. In particular, the renormalization can be perf ormed 



using counter-terms 
larization prescription 



174j. 1 1751 . Il87l | or even more conveniently using a BPHZ |164| | regu- 



190l |. The only important constraint is that the truncated system 



®Which in a perturbative construction is known as the LSZ formalism, if the asymptotic state space is 
available [30l |. 

■^^This not true for all gauges, in particular not for non-covariant gauges [s^. 
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of DSEs shows the correct quahtative transformation behavior under the renormahzation 
group. Enforcing the coincidence with leading-order perturbation theory in the ultraviolet 



in the present truncations guarantees this actually even quantitatively [37|. Results from 



lattice calc ulatio ns automatically show the correct renormahzation behavior, up to lattice 
artifacts 

It should be noted that Landau gauge is one of the gauges having optimal renormahzation 
properties in the sense that the least number of independent divergent renormahzation 
constants occur, and al l of them are at most logarithmically divergent with the cut-off in 
four dimensions 
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^ ^ Indeed, there are at most two independent renormahzation 

factors [21, Il75 |. which can be chosen, e. g-jto be the one of the running coupling and the 
wave-function renormahzation of the ghostcj. Furthermore, in dimensions lower than four 
all renormahzation constants, and thus correlation func tions, are finite, though still some 



loop-graphs are divergent and require regularization |l64l | . In perturbative diction the theory 



becomes super -reno rmalizable, though from the point of view of correlation functions it is 



actually finite |166 



A last consideration concerns the scheme dependence 0] of the results. The stan dard 
scheme used in the calculations here is the miniMOM schemj^. given in [l75l . Il93j . It 
essentially implements the condition flTSjl to lock the gluon and ghost wave-function renor- 
malizations, and a condition on either propagator at an arbitrary, non-zero momentum 
can be used to yield the remaining wave-function renormahzation. Furthermore, the finite 
ghost-gluon vertex renormahzation constant is set to one. The finite three-gluon vertex 
renormahzation constant, by virtue of the ansatze ( 180|) and (182|) . cancels out, and the finite 
four-gluon vertex renormahzation condition is dropped anyway. This leaves only the renor- 
mahzation condition for the running coupling, which is implicitly defined by the relation 

There are two additional remarks to this scheme. First, if the different non-perturbative 
gauges indeed exist as described in section [531 the propagators will differ for different choices. 



though the difference at very large momenta will be power-like suppressed |63[. Still, this 
implies that the renormahzation conditions will differ between different non-perturbative 
versions of the Landau gauge, introducing a gauge- dependence. The second is that the 
definition fl7^ for the running coupling is not unique. In particular, it is possible to de- 
fine a running coupling which will have an infrared fixed- point irrespective of whether the 
propagators show a finite-ghost or a scaling behavior 37|, |204| . This is in contrast to the 
definition (1731) . for which the running coupling is only infrared finite for the scaling case, 
but not for the finite-ghost case. Since this implies a difference between the corresponding 
/3-functions, this yields that both correspond to different renormahzation schemes. In fact, 
since the former always produces an infrared fixed-point while the latter can imply either 



^^Since possibly the resolution of the Gribov-Singer ambiguity may be related to the specification of at 
least one of the two correlators to some extent, it may be that one of the two renormahzation constants can 
be related to the non-perturbative gauge-fixing procedure. This would be satisfying, as it would leave only 
one independent condition for the one independent parameter of the theory. However, at the present time, 
this is pure speculation. 



Which is in spirit very close to the conformal scheme proposed in 234| . 
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a fixed-point or a zero crossing of tlie /3 function, as demonstrated in the next chapter, 
it follows that the scheme transformation between both schemes is singular. Such sin gula r 
scheme transformations are possible, and can affect as here the presence of fixed points 235 . 



4. Describing gluons at zero temperature 

With the chapters [2] and [3], the stage is set to actually determine correlation functions. 
In the following the results of such calculations will be presented. Also, the interpretation 
of these results in terms of physical properties will be briefly discussed. 

4.I. Correlation functions 
4-1.1. Propagators 

The simplest correlation functions which are non-zero in Landau gauge are the propa- 
gators. Before investigating their momentum dependency, a first important question is the 
justification of whether it is permitted to factorize the color structure in the form of a unit 



matrix, as it is usually done. The so-called primitiveness assumption of group theory |236 
suggests this, as does perturbation theory |2|. Furthermore, the DSE for the ghost propa- 
gator shows that the non-tree-level color structure of the ghost-gluon vertex is not relevant 
for the ghost propagator. Still, a general proof is lacking. 

Explicit lattice c alcu l ation s have also provided substantial support for color-diagonal 



propagators [89|, |90|, Ill2l . Il65| . As an example in figure [12] some of the color-off- diagonal 
components of the propagators in two, three, and four dimensions for various gauge algebras 
are shown. In all cases, within statistical errors, the off-diagonal propagators are compatible 
with zero. Thus, assuming the propagators to be color-diagonal seems to be justified, and 
therefore will be done henceforth. 

The first propagator to be studied is the gluon propagator, also as function of the di- 
mensionality, and for the various gauges and cases. The results are shown for both the 
propagator and the corresponding dressing function in figure [17] for two dimensions, in fig- 
ure [18] for three dimensions, and in figure [19] for four dimensions, also comparing different 
gauge algebras. Most of the results are from lattice calculations, though there are also some 
results provided from the available functional calculations using DSEs and FRGs in two, 
three and four dimensions. For the lattice calculations four different choices of the non- 
perturbative gauges are shown, the minimal Landau gauge, the absolute Landau gauge, and 
the maxB and minB gauge, to explore both a typical choice of 5, and extreme possibilities. 
Since the implementation of gauges beyond the minimal Landau gauge are numerically very 
expensive until now only much smaller volumes are available. 

The first result visible is that there is very little dependency on the gauge algebra, at 
least up to finite volume effec t s, and qu antitative dependencieJ^. The former turn out to be 



very similar for all cases |ll5l . Il65l . l247| |. It thus suffices to stay with the su(2) case for now, 
for which most data is available. In case of DSEs, the gauge algebra dependency is entirely 
through the combination q'^Ca at this level of truncation, and therefore the dependency 



^°Note, however, [2481 l249l | for su(2) and su(3) in four dimensfons. 
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Figure 16: Some components of the o ff-diagonal gluon propa gator (left panels) and ghost propagator ( right 
in two (top panels) [ool. Il65j |. three (middle panels) |ll2 . 165 [. and four (bottom panels) [sol . l237| 



panels) 

dimensions, for various gauge groups. The results in four dimensions are renormalized at 2 GeV. Here and 
later in the legends the groups used for the lattice implementations are given. 
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Figure 17: The su(2) gluon propagator (top left panel) and its dressing function (top right panel) in two 
dimensions. Shown are results for minimal Landau gauge (MLG) [90], absolute Landau g auge (ALG) [63l |. 
and the maxB and minB gauge [g^ . [sij . compared to DS E results in the scaling case |l87| . Note that 
in two dimensions the finite-ghost case seems n ot to exist |l87l |238| . The bottom panels show the same, 
but for gauge algebras su(3-6) and g2 in MLG 165[. Further results in two dimensions can be found in 



can be abso rbed in a dependency on the renormalization scale, and perfect 't Hooft scaling 



is manifest [l75l . Il94 |. except for perturbative corrections in four dimensions. Of course, 
by construction, leading-order (resummed) perturbation theory is reproduced at sufficiently 
large momenta for all cases. 

At first glance, there is very little difference between the propagator for the different 
gauges. This is in line with the expectations from figure [2|, which already indicates that 
the total integral of the gluon propagator is rather constrained. However, since the measure 
in (1251) strongly suppresses the infrared, this is not a necessary result, and large differences 
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Figure 18: The su(2) gluon propagator (top left panel) and its dressin g fun ction (top right panel) in three 
dimensions. Shown are lattice results for minimal Landau gauge (MLG) [240|, abso lute Landau gauge (ALG) 
[g^l, the maxB and minB gau ge 62, 81|, and DSE results for the scaling case [l66| . A finite-ghost- type DSE 



soluti on ca n be found in 201]. The bottom panels show th e same, but for gauge a l gebras su (3-6) and g2 in 
MLG [lei]. Further lattice resuhs can be found in [H, HI, [ml [ul [113, [12I Im [24l[j243j 



in the infrared could be possible [63|. This is indicated by the DSE results, which show 
a qualitative difference between the scaling and finite-ghost case, with a vanishing gluon 
propagator at zero momentum in the former case. Concerning this statement, there are two 
points to be cautious about. One is that in every finite volu me, even in the scaling case, 
the gluon propagator has to be non-zero at zero momentum 145|. As shown in figure! 



however, only in the two-dimensiona l cas e the gluon propagator at zero momentum appears 



to extrapolate to zero [62, [63|, [8l|, [9Q|, Ill7l | and in this case for all gauges. On the other hand. 



the gluon propagator appears to be non-zero in three and four dimensions in all gauges 



63l . l81l . III5I . [ml, III8I . I239L l243l | . The fact that the result of the minB gauge undershoots the 
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Figure 19: The su(2) gluon propagator (top left panel) and its dressing function (top right panel) in four 
dime nsions, renormalized to Z(1.5 GeV^)=l. Shown are lattice results for minimal Landau gauge (MLG) 
[llSf , absolute Landau gauge (ALG) , a nd th e maxB an d mi nB gauge [HI] . The bottom panels show the 
same, but for gauge algebra su(3) in MLG 237 [ and ALG llSl , c om p ared to scahng and finite-ghost DSE 
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one of the absolute Landau gauge demonstrates again the irrelevance of the zero-momentum 
behavior for the absolute-Landau-gauge condition. A comparison at finite momenta shows 
that the integrated weight of the propagator in absolute Landau gauge at the same lattice 
setting is smaller than for the minB gauge j62|. 

This seems to be ruling out the scaling case in four and three dimensions in the first 
Gribov region, but, as discussed in section I3.2.4[ the actual value of the exponent kaa is 
dependent on the truncation, and its value changes if the ghost-gluon vertex is different 
from a bare vertex in the infrared. For both ansatze employed here, the ghost-gluon vertex 
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the value of be determined 



becomes constant in the infrared. In this case [OJ, Il86 
to be approximately 0.595 in four dimensions. Thus, the value of haa is then about 1.19. In 
three and two dimensions, two possible solutions are found. In three dimensions the values 
are 1/2 and about 0.39, yielding kaa as 3/ 2 ari d 1.28, respectively. The functional equations 



yield consistent solutions for both values 166j . In two dimensions, the possibilities are 1/5 
and 0. The latter value, yielding kaa = 1 with an infrared fini te gluon propagator, appears 



to be ruled out by the lattice data, which prefer kaa ~ 1-4 [90j, |239| . However, in three 
and four dimensions it is possible to obtain a scaling solution with an infrared finite gluon 
propaga tor w ith a ghost-gluon vertex, which differs quantitatively from a bare one in the 



infrared (l87j . Hence, at most, the conclusion which can be drawn from the gluon propagator 
alone is that the truncation in three and four dimensions is possibly not viable, if scaling 
is realized in the first Gribov region. Since the ghost-gluon vertex indeed shows infrared 
variations, as will be shown in the next section, this is a possibility to reconcile scaling 
with the gluon propagator results, but a full understanding requires further investigations 
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for both lattice and continuum methods. Furthermore, the ghost propagator also has to be 
taken into account, which will be done below. In addition, there are further subtle issues 
concerning this question, which will be detailed in section 14.3.21 Thus, up to this point, no 
final conclusion can be drawn. 

The results in three dimensions in minimal Landau gauge illustrate the concept of the 
scaling window, f l79p . rather nicely. After reaching the maximum, the gluon propagator first 
decreases, over a momentum range, showing an approximate scaling behavior with a local 
K.AA larger than one, before it finally reaches an infrared finite value when the screening sets 
in. In four dimensions, this window seems to be closed, i. e. pb ~ ^ym- If in two dimensions 
at very small momenta still screening sets in, for which there is no sign ^239|], then the scaling 
window in this case would be very large. 

The second elementary propagator is the ghost propagator. It, and again its dressing 
function, is shown in figure |2T] for two dimensions, in figure |22] for three dimensions, and in 
figure [231 for four dimensions. Again, there is little dep endency on th e gauge algebra, and 
thus the following will be restricted to the su(2) case |115l . Il65l . |247| |. The ghost exhibits 
a rather stron g dependency on the gauge. The propagator is least divergent in the minB 
Landau gauge |63l. ISlL llOOl llOOl Il34| . at least for the same lattice settings. Different volumes 
and discretization artifacts in the renormalization in four dimensions can lead to an apparent 
different ordering. In the case of mi nim a l Landau gauge the dressing function is infrared 
finite in three and four dimensions 



115, 116, 118 



dimensions 90, 116 



on the volume 90, 112, 145 



but does not appear to be so in two 
I t is a lso quite visible that the infrared enhancement strongly depends 
The features of the finite-ghost case is rather well reproduced by 



the finite-ghost solution of the DSEs |37| . The behavior is thus that of a massless particlqfll. 



The situation for the maxB gauge is drastically different [8l|. In this case, the ghost 



propagator is much more divergent than even the scaling case, and more divergent than the 
results in any of the other gauges. Such an over-scaling can actually also be obs erved in 



functional calculations in a finite volume at finite cut-off for the scaling case [229|. In the 



latter case the over-scaling originates from a mismatch between continu um a nd finite- volume 



regular izat ion, and can thus be regarded as a kind of lattice artifact [145|. Given that in 
the other gauges the ghost propagator also becomes less divergent with increasing volume, 
it is to be expected that this also applies to the maxB gauge. If then in the infinite-volume 
case the ghost propagator becomes the scaling one or a finite-ghost one is currently an open 
question SjJ. Since the functional calculations naturally provide such a scaling case, there 
exists a motivation that this should be the case. However, as will be discussed in more 
detail in section I4.3.2[ it may also be that it is not possible to achieve scaling within a single 
Gribov region. This will require further investigations. 

It should be noted that the various gauges only start to differ significantly below about 1 



GeV. As to be expected [63[, in the essentially perturbatively dominated energy range above 



^^Actually, the lattice data can not exclu de the possibility of a logarithmic divergence, which has been 

lid 203| . It will be assumed here henceforth that this is not the case, 



proposed for the finite-ghost case 
as the results from most calculations are much better in agreement with an infrared finite ghost dressing 
function in the finite-ghost case. 
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Figure 21: The su(2) ghost propagator (top left panel) and its dressing function (top right panel) in two 
dimensions. Shown are results for minimal Landau gauge (MLG) absolute Landa u ga uge (ALG) [63l |. 
and the maxB and minB gauges 81 1, compared to DSE resultsin the scaling case |187| . Note that in 
two dimensions the finite-ghost case seems not to exist |l87i |238| . The bottom panels show the same, 
but for gauge alg ebras su(3-6) and g2 in MLG 165[. Further results in two dimensions can be found in 
73. [llll . 1ll6[ [l28j. 



1 GeV the non-perturbative gauge dependency does not alter the results. This indicates that 
the contributions decay faster than the leading-order perturbative contribu tions. However, 
when performing an operator-product expansion of the propagators 255l . l256| . the alterations 
would contribute. 

As repeatedly stated, these results have a significant dependency on lattice artifacts. 
Studies of such arti f acts due t o vo lume effects and discretization effects can be found, e. g., 
in 62, 0, 112 . 116 . 117 . 193 , 246 1. The major contribution of the artifacts is due to finite- 
volume effects. In the figures [T7] to dn] and [21] to [23] always very large volumes have been 
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Figure 22: The su(2) ghost propagator (top left panel) and its dressing function (top right panel) in three 
dimensions. Shown are lattice results for minimal Landau gauge (MLG) ^1 , abso lute Landau gauge (ALG) 
[g^l, the maxB and minB gau ges 81 1, and DSE results for the scaling case |l66l |. A finite-ghost- type DSE 



soluti on ca n be found in 201'. The bottom panels show th e same, but for gauge a lgebras su(3-6) and g2 in 
MLG [lei]. Further lattice results can be found in (62l. Isj. [i H [Ti l [Tlil . [iSi |24^ 



used. To illustrate the dependency on the volume, results for different volumes are shown in 
figure [211 It is is visible that the volume-dependency for the gluon propagator is apparently 
stronger than for the ghost dressin g function , making the infinite-volume extrapolation of 



the latter a rather subtle problem [116l . |239| . It is also possible to perform studies using 
functional methods in a finite volume and at finite discretization, which also show sizable 
finite volume artifacts 



145 



as discussed in section 13.2.71 
Even if the gluon propagator and the ghost propagator would be infrared vanishing and 
infrared divergent, respectively, this is not sufficient for the existence of the scaling case of 
the functional results. The characteristic property of scaling is the sum rule (ITOl) . and thus 
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Figure 23: The su(2) ghost propagator (top left panel) and its dressing function (top right panel) in four 
dime nsions, renormalized to G(1.5 GeV^)=l. Shown are lattice results for minimal Landau gauge (MLG) 
[llSf . absolute Landau gauge (ALG) [62], and t he m axB and rninB gauges 62 1. The bottom panels show 
the same, but for gauge algebra su(3) in MLG [237] and ALG [1185, compared to scali ng and finite-g host 



DSE results and scaling FRG results ,3 7 J. Further results can be found in |iiMl63.l6i. [io3l - [ioi 
12 ifl57l[l85l . fl93L[20a 



20i . [21 I [24in243 - l250[ [25i [2& 



ute- gnost 



the infrared constancy of the running couphng (173|) . This couphng is also of generic interest, 
as it is the characteristic strength parameter of the theory and the most direct contact 
to perturbation theory [193j, and thus experiment j2j. For the su(2) case, this couphng 
is shown for the four gauges and three dimensionahties in figure [251 Shown alongside is 
the volume-dependent behavior of the left- hand -side of the sum rule fl73|) . using volume- 
dependent effective infrared exponents 90l. Il45[ |. The results show a finite-ghost behavior 



in three and four dimensions for minimal and absolute Landau gauge as well as for minB 
gauge, at least for the given volumes. This apparently rules out the original idea that the 
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Figure 24: Dependency of the gluon propagator (left panels) and the ghost dressing function (right panels) 
on the volume in two (top panels), three (middle panels), and four dimensions (bottom panels) [62|. All 
results are from lattice calculations at the same lattice spacing a = 0.22 fm and in minimal Landau gauge. 
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Figure 25: The runnim 
(middle-left panel) |63[ 



coupling ([75)1 for the gauge algebra su(2) in two (top- left panel) [63, [fill |90| , three 



dimensions for the minimal, absolute, and 



l8l| . and four (bottom- left panel) [6S 
maxB gauge. For comparison, results from DSEs for scaling in three dimensions are also shown jl6t|. In the 
right-hand plots the corresponding behavior of the sum-rule (j70l) is shown. Note that the lowest momentum 
points are not used to calculate the exponents 9Q] . Solid lines give the value of the sum rule in the scaling 
case and dashed lines for the finite-ghost case. In two dimensions both lines coincide. 
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absolute Landau gauge could be connected to the scaling case [37|, l63|, |65|, 1243] , at the very 



least for the volumes investigated so far. Howe ver, in any finite volume a finite-ghost-like 



behavior is expected for the running coupling [145[. Hence, only a detailed analysis and 
extrapolation of the sum rule (ITOll . shown in the right-hand side of figure [25| can distinguish 
between the different cases. Indeed, the sum rule appears to confirm the statement that 
scaling is at best seen only in two dimensions, but neither in three nor in four dimensions, 
as has also been established indirectly on large lattice volumes including only relatively few 



points for the ghost propagator at low momenta |115Hll8l |. In two dimensions the situation 
is rather special, since the sum rule coincides for both the finite-ghost case and scaling. Thus 
here the propagators, showing in all cases a scaling-like behavior, have to be included in the 
discussion as well, as detailed above. 

The situation looks somewhat different when going to the maxB gauge. In both two and 
three dimensions the coupling shows over-scaling, i. e., it appears to diverge. As noted above, 
this may be a lattice artifact. Still, it implies that because the distribution of b is continuous 
over the gauge orbit, see figure El in all volumes studied the value of B can be adjusted such 



that a scaling-like behavior is obtained |8ll . It remains to be seen whether this is also 



possible on significantly larger volumes. In four dimensions, the coupling seems to saturate, 
though this requires also a more systematic investigation of the systematic errors: Since in 
four dimensions the number of Gribov copies at fixed lattice extension is largest, see figure 
[H it cannot yet be excluded that over-scaling occurs also in four dimensions. Given the 
experience with including more and more Gribov copies in the maxB gauge, see figure [HI 
this appears to be at least a viable option to be studied further. Nonetheless, over-scaling 
in the infinite- volume limit would be an unexpected result. 

Since the running coupling is both gauge-dependent and scheme-dependent, it is inter- 
esting to also investigate the underlying /3-function. The latter has the advantage that 



it, at least in principle, can be connected to scattering amphtudes j42|], though a scheme 
dependency still remains. The /3-function is defined as 

= (89) 

where g is connected to the running coupling (175]) by = Air'^a. The resulting /3-functions 
are shown in figure [26l It is well visible that in next-to-leading order (NLO) perturbation 
theory describes the running of the coupling rather well up to g around five, correspondingly 
a around 0.6, for the lattice data. The functional results are only up to leading order (LO) 
consistent with perturbation theory, and therefore start to deviate quicker from the NLO 
result. It is also well visible how the /3 function approaches the Gaussian fix-point in the 
ultraviolet. In the infrared for the minimal and absolute Landau gauge lattice data and the 
finite-ghost functional results a zero of the /3 function emerges at a finite g, corresponding 
to the maximum in the running coupling. At this point the coupling, however, does not 
stop running, but decreases once more, and the /3-function becomes positive. This indicates 
that the derivative of the /3 function at this point is non-zero, and becomes dominant. In 
the very far infrared it appears that the /3-function reaches another, likely trivial Gaussian, 
fix-point when the running coupling becomes zero. It would be highly interesting whether 
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Figure 26: The /3-function in four dimensions for su(2) in the minimal Landau gauge (top-left panel) an d th e 
maxB gauge (top- right panel) (g^I, and in absolute Landau gauge for su(3) in the bottom- left panel 118 1, 
compared to NLO perturbation theory [ij, and the functional continuum results for su(3) for the finite-ghost 
and scaling case in the bottom- right panel [STj . Both axis have been rescaled by an arbitrary constant for 
visibility, and therefore the absolute values should not be taken from the figures. 



this occurs in an analytic or non-analytic way, but this has not yet been investigated. In the 
maxB gauge a zero-crossing is not yet seen, but much larger volumes are required to check 
whether this remains true. Of course, for the scaling solution, the /3-function has indeed a 
non-trivial infrared fixed-point, but the running coupling also exhibits a shallow maximum, 
and thus the /3-function arrives at its infrared fixed point from above. This can esse ntially 
not be resolved on the scale of the figure [261 cind could be a truncation artifact |l73 . 



Thus, from these investigations it should be concluded that in the infinite-volume limit 
both the finite-ghost case and scaling case can be found using functional methods. In the 
minimal Landau gauge, within systematic error, in three and four dimensions a finite-ghost 
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behavior is obtained, and in two dimensions a behavior which is at least very close to scaling. 
The same appears to be the case for the absolute Landau gauge and the minB gauge, though 
in these cases the systematic uncertainties are larger due to the higher computational costs. 
Finally, the results for the maxB case are not yet conclusive, and require further systematic 
studies. Of course, the lattice statements apply to calculations in the first Gribov region. 
Whether this marks the existence of the scaling case in three and four dimensions an artifact 
of the functional equations is a subtle question, and will be discussed in more details in 
section I4.3.2[ and requires that the non-perturbative gauge-fixing procedure for functional 
equations described in section 12.5.41 and 13.21 indeed lead to an equivalent gauge in both 
functional and lattice calculations. 

Before continuing, it should be noted that the functional results show by construction the 
correct perturbative behavior. By comparing it to the lattice results, which automatically 
contain up to lattice artifacts all orders of perturbation theory, it shows that already at 
momenta around 1.5 GeV resummed leading-order perturbation theory becomes an adequate 
description of the propagators within a few percent. Given that there are proces ses fo r which 



a perturbative description is not valid even for substantially higher energies [257|, this is 
rather early. On the other hand, the onset of Gribov effects at momenta starting from 
around 1 GeV shows that significant non-perturbative effects come into play already very 
close to the perturbative domain, leaving only a rather small window where higher-order 
perturbation theory can substantially improve the description of the propagators. 

In total, the results for the propagators exhibit clearly different behaviors, depending on 
how the non-perturbative gauge-fixing is treated, being it either in the form of explicit Gribov 
copies in lattice calculations or of undetermined parameters in the functional equations. The 
presence of the finite-ghost case is definitely established for the minimal Landau gauge, and 
rather well reproducible by functional calculations for the same choice of B in both three 
and four dimensions. 

4.1-2. Vertices 

The next-simple correlation functions beyond propagators are the three-point vertices. 
Two of them exist in Landau-gauge Yang-Mills theory. One is the ghost-gluon vertex, the 
other one is the triple-gluon vertex. 

Of both, the ghost-gluon vertex r^'^™'"^ i s th e significantly simpler one. It is related to 
the corresponding correlation function b\j^ |l43 | 



< A';ip)c\qyik) >= D^,^\p)D'^'^{q)D"^'^fTf'^''^f{p, q, k)5{p + q + k). (90) 

Because of Lorentz-symmetry, the ghost-gluon vertex can be expressed in terms of two tensor 
structures 



170 



Ff -^^(p, q, k) = ig {q,A^'%p, q, k) + p,B'^'%p, q, k)) , (91) 

which are chosen here conveniently to be one along the anti-ghost momentum, and one 
along the gluon momentum. At tree-level, A"'^'^ = /"^^ and B"'^'^ = 0. Due to the explicit 



Connected and full correlation functions coincide in Landau gauge for three-point functions 112[. 
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(transverse) gluon propagator, of the two possible Lorentz structures only the one not longi- 
tudinal in the gluon momentum appears in the unamputated correlation function ( 190|) . and 
in particular in functional equations projected transversely. So there is only one relevant 
tensor-structure to be kept. At tree-level, this structure has a color-structure of f""^^. In 
principle, others can appear beyond tree-level. This has so far only been investigated to some 



extent, indicating the appearance of a d"'"^ stru ctur e for su(A^ > 2) gauge algebras [198 



236[. However, since only the component 



i. e., for the only ones with non- vanishing d"'^'^ 
proportional to /'^'"^ is relevant for the functional equations at propagator level, here only the 
component proportional to the tree-level color structure will be discussed. This leaves one 
single scalar dressing function. In all gauges implementing the Landau gaug e condition 



it can be shown that the associated renormalization constant is finite 170 . Furthermore 



beca use of a ghost-anti-ghost symmetry in Landau gauge [l^, which is presumably intact 
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the ghost and anti-ghost legs can be exchanged without changing the value of the 
dressing function, i. e., A{p, q, k) = A{p, k, q). 

Lattice results in two, three, and four dimensions for su(2) for this dressing function in 
various kinematic configurations are shown in figure [27] for the minimal Landau gaugj^. 
In all cases, the vertex dressing only slightly deviates from tree-level at mid-momentum, 
and becomes comparable to tree- level in both the infrared and ultraviolet, though might 
be slightly smaller than tree-level in the infrared. The enhancement at mid-momentum in- 
creases with dimensionality, but never exceeds about 50%. Thus, the ghost-gluon vertex is 
indeed very close to tree-level, as commonly assumed in functional calculations. However, 
a comparison with figure [13] immediately shows that the ansatz to remove quadratic diver- 
gences from the ghost loop is different in the ultraviolet. From this it can immediately be 
deduced that the removal of these divergences in the gluon equation has, in principle, to be 
done with the necessary counter-terms for the non-gauge-invariant regularization, instead of 
some modification of the vertices. 

For the scaling case, functional results are availabl e, w hich are obtained in a semi- 
perturbative manner, i. e., by using a skeleton expansion 181 [but replacing the propagators 



with the non-perturbative ones from the previous section |l53l |. Thus, these results are not 
self-consistent, and will only give a first estimate of the non-perturbative corrections. The 
result is also included in figure [27] These are found to not devi ate s ignificantly from tree- 
level as w ell. Even when modifying the input vertices drastically 186| , this does not change 



153, 169 



In four dimensions, r esult s for su(3) are also available, from both lattice 2581 . l259| and 
continuum calculations |l53j. They show the same qualitative behavior, and almost the 
same quantitative behavior. 

In summary thus the ghost-gluon vertex is essentially constant. Since this is also the case 
in two dimensions, it can be expected to be the case in both the scaling and the finite-ghost 
case. As noted, this is also anticipated in t he scalin g case from very general consistency 



arguments between the DSEs and the FRGs |l59l . Il60| . The slight angular variations at low 



momenta may however be of quantitative importance, as they could alter the value of the 



There are not yet any results in other gauges available. 
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Figure 27: The tensor component A of the ghost-gluon vertex for two (top panels) [o^ , three (middle panels) 
[sot - and four (bottom panels) [s^ dimensions, in minimal Landau gauge. In three and four dime nsions, 
results from a semi-perturbative evaluation of the DSEs in the scaling case are also shown |l53| . The 
left panel shows the result in the symmetric momentum configuration ~ = fc^ for three and four 
dimensions. In two dimensions, this configuration is not possible on a quadratic lattice except at special 
points (ooj . and the result for zero gluon momentum is shown instead. The right panels show the results for 
the gluon momentum orthogonal to the ghost momentum for the largest volume shown in the left panel. 
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exponent Kgc |l87l |. In particular, such changes can decrease this exponent, and therefore 



move it closer to one generating an infrared finite gluon propagator even in the scaling case 
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The three-gluon ver tex i s a much more complicated object, due to its much more involved 



Lorentz structure 0, [ll^, EM, H-[263 . It is related to its full correlation function by 



112, 260 



< A^^ip)Aliq)A;ik) >- 



D^y{p)Dt;'\q)Dj;^fvff;f{p, g, mp + q + k). 



(92) 



The correlation functi on h arbors four independent Lorentz structures, which are transverse 
in all gluon momenta 144| . and on top one or more color structures. The latter have been 
investigated, and it has been found that at least one other non-trivial structure can exist in 
principle [73|, Il98| . The various structures have been investigated using functional studies 



in the far infrared in the scaling case [15J], showing differences not too large for the various 
Lorentz structures, though this remains a subject of further investigation. 

In lattice calculations, only t he proiection on the tree-level color and Lo rentz structure 
has been investigated so far 89, 112 . 260 - 263l |. which is given in general by jll2| 



A^^^ip,q,k) 



^X'^''''''{p,q,k) < Alip)Aliq)A<^^ik) > 



■^t\,L,A3,abc 



ip)m{q)D%{k)Tl'^^'''^^{p,q,k) 



(93) 



The in dex L at the tree-level vertex denotes the usage of the lattice ver sion of the tree-level 
vertex 1371 ]. to reduce artifacts from violation of rotational symmetry 112| . By definition, 
the function would be one if the full and the tree-level vertex coincided. It should be 
noted that this is the only contribution of the three-gluon vertex which contributes in the 



gluon loop of the gluon propagator's DSE |112| . Results in two, three, and four dimensions 



are shown in figure [28] for su(2) in the minimal Landau gauge. The results are v ery noisy 



and thus very high statistics are necessary for an accurate result [89|, |90|, Ill2| . making 
the determination of this vertex a computational challenge. The results finally show that 
the qualitative behavior is rather similar for all dimensions. In particular, in all cases a 
suppression at momenta around the typical scale Aym of a few hundred MeV is found. In 
two and three dimensions a sign change occurs at about 100-300 MeV. In four dimensions, 
data at these momenta are not available yet, but the tendency is of the same type. 

In the very far infrared, in two and three dimensions a very strong enhancement (with 



negative sign) is found 89|, |90|, Ill2| . In fact, the behavior is very much reminisc ent of th e 



infrared divergence found in the scaling case for the vertex in functional studies |l8ll . Il82 



Since again in the two-dimensional case minimal Landau gauge is at the very least very close 



to the scaling case, it is even possible to quantify the infrared exponent |90j , which turns out 



to be in agreem ent i r i wit h the formula fl78p . confirming the analysis of all vertex functions 



in the infrared 18ll . Il82| . This is shown explicitljo in figure Remarkably, the fitted 



^''Note that the momentum configuration is strictly speaking not the same as used in the corresponding 
calculations using functional methods, which would be at the symmetric point. However, for the Bose- 
symmetric three-gluon vertex, this could be of minor importance. 
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Figure 28: The tensor component ([Ml) of the three-gluon vertex for two (top panels) foo'l, three (middle 
panels) (soj . and four (bottom panels) [ 89,] dimensions, in minimal Landau gauge. Momentum configurations 
are as in figure [?fl Points with an error l arger tha n 1 have been suppressed. Some results in four dimensions 
for gauge algebra su(3) can be found in |26Cll4263j |. for one gluon momentum vanishing. 
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Figure 29: The tensor component of the three-gluon vertex in two dimensions for the gauge algebra 
su(2) in the minimal Landau gauge at low momenta (90| . The points come from different lattice volumes 
and discretizations. The line corresponds to a function 0.17p~^'^, the exponent coming from (|78p . 



)re-factor is about an order of magnitude smaller than the one of the propagator power-law 
90i . Thi s has also been found in semi-perturbative functional calculations to be the case 
Moreover, when determining the four-gluon verte x wit h the same method, the 




pre-factor is found to be a further order of m agni tude smaller |155| . while at the same time 



showing the self-consistent exponent of fl75]) 181 



Though in higher dimensions it cannot be expected that the enhancement is a genuine 
divergence due to the screening occurring in the propagators at very low momenta, at the 
very least, the triple-gluon interaction strength in this channel is strongly enhanced at low 
momenta. In the scaling case, it has been found to be infrared divergent in all dimensions 



18ll . Il82j , as w ould be all other tensor structures, as an infrared analysis of the functional 
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equations |181| . semi-perturbative studies |l54j |. and approximate self-consistent treatments 
have shown 



Calculations of the four-point functions using lattice gauge theory have not yet been per- 



formed due to the strong statistical noise in such higher-order correlation functions [90|, |112 



An even more restricting problem is that disconnected contributions have to be removed ex- 
plicitly, beginning with the four-point functions. Functional calc ulati ons have started to 



explore higher n-point functions using a semi-perturbative ansatz |155| . Beyond this point 



in both lattice and continuum calculations usually only bound-state properties have been 
calculated, either directly from the correlation functions in lattice calculations from their 
exponential decay in Euclidean time j^, Hsj , or using Bethe-Salpeter or Faddeev equations 



90 



in functional methods 13|, 1264 |265| . Rather recently, dynamic hadronization has also 
received attention for an effective possibility t o in clude higher-order correlation functions 



in way of a self-consistent ef fective th eory |39l . |266| . For the only physical bound-states in 



Yang-Mills theory, glueballs |267l |268| . funct ional met hods have met quite severe technical 

while lattice studies have been quite 
Further objects of interest so far are 



problems, which are still a matter of resear ch |269l-l272 



successful in determining their spectrum [273 
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scattering processes and bound-state coupling constants, though their direct determina tion 
in terms of correlation functions becomes increasingly complicated 



Even with the yet limited range of higher order n-point functions results, conclusions can 
be drawn both for physics and for technical aspects. From the technical point of view, the 
fact that the ghost-gluon vertex should be enhanced at mid-momentum and this particular 
tensor structure of the three-gluon vertex should be suppressed, has been inferred from 



functional studies under the assumption that the two-loop terms are negligible [166|, |222 
The reason is that otherwise the opposite signs of the ghost and the gluon loop destabilize 
the system, and in particular tends to drive the gluon dressing function into the negative, 
which is by construction not permitted. This is now confirmed a-posterior to be correct by 
the results shown in this section. 

The sign change of this particular tensor structure of the three-gluon vertex also has quite 
far-reaching implications. Since it is the only contribution relevant for the gluon loop in the 
gluon propagator equation at one-loop order, it can provide a sign change of the whole loop 
in the infrared. Thus, the gluon physics changes from being ant i- screening at large momenta 
to screening at small momenta. In accordance with the observations in section I3.2.4[ the 
screening mass from the ghost-loop in the critical limit of Hgc = (c? — 2)/2 is then enlarged 
by the gluon-loop, as both now have the same sign. This is also of profound importance for 
the finite-ghost case, where the gluon loop is one of the leading contributions, and the sign 
change helps to provide a positive screening mass. 

The second consequence is even more profound. The three-gluon vertex behavior shows 
that at the hadronic energy scale the emission of a gluon from a gluon is strongly suppressed, 
at least for the tree-level tensor. However, it becomes strongly enhanced at momenta close 
to zero, i. e., close to the the light-cone. At these small momenta the gluons are screened 
by either a finite or an infinite mass, thus not contributing anymore. This could be an 
important hint into how the infrared enhancement threatening to be unit arity- violating in 



parton structure functions could be remedied by non-perturbative effects [276 



4-2. Schwinger functions, mass, and confinement 

Having the propagators now available, it is an interesting question what kind of insight 
can be obtained from them on the properties of the elementary particles. In particular, their 
masses and widths are of great interest, as is their analytic structure. This information is 
encoded in the propagators, and can be obtained by determining their space-dependence. 

The most direct access to properties like masses and widths is granted by the Schwinger 
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function 32 . This function is defined ai 



1 f°° 

A{t) = - / dpocos{tpo)D{pl) (94) 

Jo 

for a particle with propagator D[p^). The notation po indicates that the dressing function is 
eva luate d at zero spatial momentum. On the lattice, the corresponding expression is given 



by [217 



^ Po=0 ^ ^ 

where both t and Pq are in lattice units, i. e., integer, and Nt is the lattice extension in time 
direction. Note that the sum extends over the whole momentum range, and includes the 
parts of the propagator reproduced by periodicity. 

In case of a stable particle with a simple pole mass M, having an Euclidean propagator 

D{v) = -T^^ (96) 
and thus a pole at p = ±im, the Schwinger function is given by 

AW = (97) 

Therefore, a simple exponential decay is expected in Euclidean time for a massive, stable 
particle. It is worthwhile to note for the investigations at finite temperature in chapter |5] 
that such an exponential decay for a bosonic, stable particle with zero width is also expected 
in the spatial directions for the soft mode, i. e., for A(z), since its propagator is l/{0^ + m'^). 
Here, the energy is zero, as it is a soft mode, as are the other spatial momentum components 
orthogonal to the direction along which the integral to determine the Schwinger function is 
performed. 

Beyond tree-level, the propagator fl96l) is modified by the appearance of a cut, starting 
at p = ±2im [30]. At lowest order in perturbation theory, a possible analytic form for such 
a propagator in four dimensions is given by 

= 2 j_ 2 ^r,f 2 2^ (98) 
n(p^m^) = -g^ [ — p + WlH — atanh ' ' 



2v^ V I y + p2 

where g has dimension of mass. This form is motivated by leading-order perturbation theory. 



and occurs, e. g., for a scalar theory with a three-point coupling [278 



'^^The Schwi nge r fun ction can also be evaluated directly in position space without the detour over mo- 
mentum space 43, 273| . 
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If the particle is not stable, the poles are moved off the first Riemann sheet onto the 
second Riemann sheet at m + iV /2. The cut then starts at ±2iM, where M is the mass 
of the particles in which the original particle can decay, assuming for the mome nt on ly this 



two particles in the theory |30|. In this case, the propagator (p8|) is modified to |278 



(u A^) - n ^ (^m 



(99) 



n(p^m^A^ 



4m2 

1 H --atanh^ 

p2 



If A2 
+ -ln( 1 + — 



pz 



{2JS? + + p^) atanh 



4A2+4m2+p2 



4m2 + p2 



(A^ + m2) +p2 



where there are now two couplings of dimension mass, g and h, describing^ a self-interaction 
and a decay channel with strength h, and A is the cutoff, of which this renormalized 
propagator is independent. 

However, since the gluon (nor the ghost) behaves necessarily like a physical particle, after 
all it is not gauge-invariant and thus not physical, its propagator may not be of either form 
fl98|) or fl99l) . There have been va rious prop osals, which form it may have instead. One, the 



so-called Gribov-Stingl type \uX l279l - l28lj . h as comple x poles on the first Riemann sheet 



but may have vanishing residues at the poles |2 791428 This behavior can be described by 



a meromorphic function |32 



D{p) 

Its Schwinger function is given by 

e2 

A(t) = — -( 

^ ' 2m3 sin (20) 



e2 + /p2 



+ 2m? cos (20) p2 -)- ^4 



~tm cos((/)) 



sm I 



ffn? 

+ tm sin (0)) H — sin 



— tmsin (0)) 



(100) 



(101) 



The Schwinger function therefore exhibits oscillation with a period determined by the angle 
and the mass parameter m. Its positivity violations are signaling therefore the instability 
of the particle. The propagator fllOOp illustrates at the same time the difference between 
the concept of screening mass and pole mass. The screening mass is defined as the inverse 
square-root of the propagator at zero momentum, and given by 



Dm 



m 

e 



(102) 



^^Note that not all parameters are independent. The reality of the Euclidean propagator fixes one of the 
parameters m, F, M, 5, and h as & function of the others, since the decay width is not an independent 
quantity. 
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Since e depends on the wave-function renormalization, this immediately shows that a screen- 
ing mass is never renormahzation-group-invariant, and can therefore not be a physical ob- 
servable quantity. The pole mass, i. e., the location of the poles of fllOOp is given by 

ime^''^, (103) 

and is thus a set of complex conjugate poles, but in general with non-zero residuum. This 
mass is renormalization-group invariant, if there is no additive mass-shift to m, and can 
thus be, at least in principle, a physical mass, though in the present case it may be gauge- 
dependent. Note that in the limit of real masses, e has to go to m and / to one, or otherwise 
a double-pole would emerge. 

There is another concept which appears when using the Euclidean correlation function. 
Assume that the propagator can be more generally written as 

D{p) = — (104) 

with some wave-function renormalization constant Z. At first sight, the condition = 
—MijP'Y could look like an indicator for the pole mass of the particle described by the 
propagator. However, this is only possible if the function M(p^) only depends on the real 
part of p^. This is certainly the case for (1^ . but is not necessarily the case for ( 11041) . 
Thus, the point —MijP'Y only gives a would-be pole mass, and the correct pole mass is 
only obtained when taking into account the full dependence of M{p'^) on complex momenta 



32, 37 



None of these forms can actually fit a scaling type behavior. The results from functional 
studies suggest a form which has a cut on the imaginary momenta axis starti ng a t zero 



momentum. Functional forms which provide such a structure are given bjcj |32|, 1166 



= (i+/;^_,rt..>. /"v(/). (105) 

where /uv encodes the perturbative tail. While the first version works better in three 



dimensions |166l |. the second is more adequate for four dimensions [3^. Both have no simple 
complex or real poles on the first Riemann sheet, but have additional poles on further 
Riemann sheets. If one assumes that the finite-ghost case is a modification of the scaling 
case, a fit form like 

^'^P^ = ail + fp^ + gp^^^^^J '^^^P ^'''^ 

is suggested. This form can retain the intermediate scaling behavior in a scaling window 
(179]) to some extent. 



Note that similar functional forms are already encountered in QED [Uliil, but also appear in (near- 



)conformal theories, as are considered for unparticles 2,83 ] , and to some extent technicolor models in beyond- 
the-standard-model scenarios Q. 
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Figure 30: Various forms of propagators in the Euclidean domain and of the corresponding Schwinger 
functions. To illustrate lattice artifacts, the Schwinger functions are also shown for the case of a 64^ lattice 
with a = 0.13 fm. The parameters for the propagators are: tree-level propagator (|96p m — 1 GeV (full line); 
double pole (fTUO)) m = 1.03 GeV, e = (1.03 GeV)^, / = 0.9, (j) = 0.245 (dashed hue); stable particle with cut 
l^m^l GcV, g = 70 MeV (short-dashed line); unstable particle with cut ([Ml) = 1 GeV, T = 250 MeV, 
M = 300 MeV, g = 470 MeV, h = -484 MeV (short-dashed-dotted line); intermediate scaling form (1X06)) 
m = 1 GeV, KAA = 1.19, a = {1 GeV)2+2«A^, / ^ (1.3 GeV)^, g = 1 GeV-2-2«^A^ = 1 (dashed-dotted 
line); scaling ([TCT5|) kaa = 1-19, / = -0.996, = (1 GcV)-2«aa^ ^ 1 (dashed-triple-dotted line). 



Most importantly, all of the forms permit Wick rotation [32|, essential to transfer the 
results back to Minkowski space. To illustrate the differences, both the propagators and 
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Figure 31: The Schwinger function for the gluon in 
two (top- left pan el) [H, |8l|, , three (top-right panel) 
(63 , 24Cll . 277| , and four (bottom panel) dimensions 
[62, ISSl, for the minimal Landau gauge, the absolute 
Landau gauge, the maxB gauge, and the minB gauge 
for su(2). Also results from DS Es ar e shown for the 
scaling case in three dimensions (l66| . Functional re- 
sults for four dimensions can be found for both the 
finite-ghost case and scal ing in 37[ and further lattice 
results also in 28^ l285l | . These are all qualitatively 
and quantitatively very similar to those shown. 



the corresponding Schwinger functions for all of these forms are shown in figure [301 It is 
nicely visible that the form of the Euclidean propagator is very similar in all cases, and only 
the Schwinger function shows significant differences. The stable form, with or without cut, 
shows an essentially exponential decay. The unstable particle has a zero at 1/(2M), and 
decays afterwards slower, due to the lighter mode in the propagator. Still, this long-time 
decay is essentially exponential. The double-pole structure shows an oscillatory behavior, 
superimposed with an exponential decay. The cases with full or intermediate scaling both 
also show a zero-crossing, but no oscillations, and a decay which is not exponential in 
time. Thus, it is the long-time behavior which ultimately is the most powerful possibility to 
distinguish between the analytic structure of the possible propagator forms. 

After this rather general introduction, the Schwinger function can be obtained for the 
different propagators in Yang-Mills theory. The first object to apply this formalism to is 
the gluon propagator. As figure [H] shows, the Schwinger functions for the various cases are 
qualitatively identical. Furthermore, the DSE results for both the finite-ghost and scaling are 
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also qualitatively identical [37] . However, the results are strongly affected by finite- volume 
corrections, and it requires a certain minimal volume to observe the most characteristic 
feature, the zero crossing of the Schwinger function at a range of a half to a few Fermi, 
depending on the number of dim ensio ns, the gauge algebra, and whether the behavior is of 
finite-ghost type or scalinj^ 62 . 277 1. Besides this first zero crossing, no further crossings 



have been observed in functional calculations, and thus there is no oscillatory behavior, as 
would be expected for a double pole structure. Furthermore, the decay at large times is 
not found to be of exponential time for both cases, thus making an analytic structure like 
(11061) the best description. For the lattice calculations, statistical uncertainty makes it yet 
impossible to decide which of the possible forms with a zero crossing, f l^ . fllOOl) . or (11061) . 
is most appropriate. In any case, the gluon Schwinger function is negative above a certain 
distance. 

This implies that positivity is manifestly violated for the gluorJ^. This implies the 



absence of the gluon from the asymptotically physical Hilbert space [10|, l32|. However, 
positivity violation is not equivalent to confinement. The case of an unstable particle (19 9 p 
shows that positivity violation is not implying confinement. Whether confinement is im- 
plying positivity violation is not known. However, a confined particle is necessarily absent 
from the asymptotic physical Hilbert space. Thus, the positivity violation implies already 
this necessary condition. 

Note that this does not prevent the existence of a stable asymptotic gluon state in the 
unphysical part of the Hilbert space of the theory; it is just excluded from the physical 



part of the Hilbert spac^J 12 1. Furthermore, the screening mass found in the finite-ghost 
case is not implying the existence of a pole mass of the gluoiF^. as illustrated by the form 
(I106p . compatible with the results from functional calculations. The only implication of the 
positivity violation is thus just that the gluon will not appear as an asymptotic state in the 
physical part of the Hilbert space. 

Though not linked to a particle, it is also interesting to investigate the Schwinger function 
of the ghost. A complication in its determination is that the propagator diverges at least as 
fast as the one of a massless particle. Thus, its Schwinger function has to be renormalized. It 
is shown in figure [321 and it is visible that is suffers from quite strong lattice artifacts and/or 



has a significant gauge dependence. This has to be investigated further [62[. Nonetheless, 
after taking into account the divergence at zero momentum by appropriately normalizing 
it, it is found that the Schwinger function becomes constant at long distances, instead of 
decaying to zero as for ordinary particles or the gluon. Thus the ghost implies correlations 
over arbitrarily large distances. This automatically implies the absence of a mass gap in 



the spectrum of the full Hilbert space [10|, ll3| , despite the presence of one in the physical 



^^Similar results have also been obtained using an effective action approach [6J,|70|. 
^^Note that a gluon propagator which vanishes at zero momentum is implying (maximal) positivity vio- 
lation (sij . 

^''This observation can already be indirectly infe rred from the violation of the Oehme-Zimmermann su- 
perconvergence relation in perturbation theory [l3[ |286| . 

^^If one exists, it would be be gauge-parameter-independent, as described by the Nielsen identities |287l |. 



though in general not gauge-independent 287[. 
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Figure 32: The renormalized Schwinger function for 
the ghost [sll in two (top- left panel), three (top- 
right panel), and four (bottom panel) dimensions, for 
the gauge algebra su(2), and for the minimal Landau 
gauge, the absolute Landau gauge, the maxB gauge, 
and the minB gauge. 
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subspace. In particular, this implies the failure of cluster decomposition for Yang-Mills 
theory on the level of gauge-dependent elementary particles and the absence of a 'behind-the- 



moon'-problem [12| . The presence of these infinite-range correlations are then a consequence 
of the necessity to maintain gauge invariance over long, and also space-like, distances. The 
Schwinger function itself is then also negative at large distances. Thus, the ghost trivially 
violates positivity, and also does not belong to the physical state space. 

It is another interesting question to investigate the Schwinger function of the coupling 
constant. Though it does not describe a particle, this is essentially a measure of the in- 
teraction strength as a function of distanc^. The result, removing trivial momentum 
factors before forming the Schwinger function, is shown in figure [331 At short distances, the 
Schwinger function is relatively large, while at long distances, even in the two-dimensional 
case where the coupling freezes out, the corresponding Fourier transform vanishes. However, 



'In QED, it corresponds to the potential, but this is no longer the case in a non-Abelian theory [Sj 
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Figure 33: The Schwinger function for the couphng 
constant [b^I in two (top- left panel), three (top- right 
panel), and four (bottom panel) dimensions, for the 
gauge algebra su(2), and the minimal Landau gauge, 
the absolute Landau gauge, the maxB gauge and the 
minB gauge. 
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also here lattice artifacts may play a role. The Schwinger function is always positive, despite 
that the propagators from which it is formed violate individually positivity. Furthermore, 
the Schwinger function also exhibits distinct minima, even when taking the relatively large 
errors into account. Whether these have any meaning is at the current time rather unclear. 
Finally, the Schwinger function for the maxB gauge differs from the others, which is likely 
linked to the fact that it is still infrared divergent for all lattice volumes investigated here. 

4-3. Confinement scenarios 

Explicit violation of positivity, manifest in the Schwinger functions, shows that gluons 
are not part of the physical asymptotic spectrum. This makes any kind of asymptotic gluon 
states part of the unphysical Hilbert space. Assuming that the S matrix only transforms 
physical states into physical states, this implies the absence of gluons from the asymptotic 
observable world. Given that this is a necessary assumption to prevent a gauge anomaly to 
occur [11], this assumption appears reasonable. 
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The emergence of gluons as rather well-defined, quasi-massless states in high-energy 
scattering experiments is not at odds with this result. Indeed, the gluon propagator behavior 
at large energies is dominated by the perturbative part, or by the tree- level part in lower 
dimensions. Thus the corrections to the perturbative picture of a well-defined gluon state are 
sub-leading at sufficiently high energies, and only manifest themselves at energies where the 
effective interaction permits to resolve the propagation of the gluons over sizable distances. 
This duality of a highly non-perturbative, screened state at low energies and an almost 
perturbative state at high energies is again a consequence of asymptotic freedom. However, 
even at the shortest distances the gluon has always a non-perturbative dressing, which, 
though not quantitatively relevant, ensures its confinement as a necessary consequence of 
having only gauge-invariant observable states: There is no almost local definition of a gluon 
in contrast to the case of a photon 
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However, this is not showing that gluons are confined. Thus, it is worthwhile to compare 
the results for the gluons to scenarios describing a confinement mechanism, which yield pre- 
dictions for the correlation functions. Three such scenarios have been of particular relevance 
recently, the ones of Kugo and Oi ima fli.l288i]. of Gribov and Zwanziger [17|, |65|, l97| , and 

There exist also othe r sce narios which give predictions 



of Gribov and Stingl [17|, I279l - l281 



on the analytic structure of the correlation functions [13|, l289j , though those seem not fully 



compatible with the results obtained, and thus will not be detailed further here. 

There are also many other proposed mechanisms, but most operate on different entities, 
like collective (topological) excitations 
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^ ^ Relations to them will be detailed more in 

section 14.41 

In all of the following it should be permanently kept in mind that there is no full sys- 
tematic control of any of the methods employed to determine the correlation functions, and 
thus the following is an evidence-driven discussion. 

4.3.1. Perturbative BRST 

Before starting with these scenarios, it is worthwhile to reconsider how in perturbation 
theory the unphysical degrees of freedom, like ghosts, are removed from the spectrum, and 
therefore loosely speaking confined [I, [l3]- This mechanism at the perturbative level is 
not sufficient to confine transversely polarized gluons, leaving the illusion that these are 
physical. Therefore, it cannot be sufficient beyond perturbation theory. The only exception 
is two dimensions, where the perturbative mechanism works, but this is due to the lack 
of transverse degrees of freedom of the gluon in two dimensions, rather than a qualitative 



change [90 



In Landau gauge, the perturbative removal of ghosts and longitudinally and time-like 



12 . This mechanism utilizes the 



polarized gluons is achieved by the quartet mechanism 
BRST symmetry of the perturbatively quantized theory to structure the complete space in a 
physical and an unphysical subspace, with the ghosts and longitudinal and time-like gluons 
residing in the unphysical subspace [2, lul ■ 

The corresponding BRST symmetry transformations are ^ 



?BRST 



6XDfct{x) 



(107) 
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^brstC^Ix) = -6X^g,r'%{x)c,{x) (108) 
^BRSTc'^lx) = SX^d^Al^ix), (109) 

where for the Landau gauge the appropriate hmit ^ — )■ has to be taken and SX is an 
infinitesimal constant Grassmann parameter. These transformations leave the gauge-fixed 
Lagrangian ([7]) invariant. This defines the BRST-operator s as 

SbrstF = 5XsF, (110) 

by its action on any field F. 

It is possible to linearize the transformation by introducing an auxiliary field, the Nakanishi- 
Lautrup field 5" pjj. This manifestly establishes the nil-potency of the BRST transforma- 
tion 

'^Irst = 0. (Ill) 

Without this field, the nil-potency is manifest only on-shell. 

To find the corresponding algebra, it is necessary to include another symmetry: Rescaling 
the ghost fields by a scale transformation exp(s) and its anti-field by exp(— s), leaving all 



other fields unchanged, is a symmetry of the Lagrangian [11|, Il2|. This gives rise to the 
conserved ghost number Qg, in analogy to the fermion number. As the ghosts are the only 
fields carrying them, it is necessary that all observable final states must have ghost number 
0. Then, counting ghosts shows that the BRST charge Qbrst has ghost number 1. Together, 
this establishes a closed algebra 

{Qbrst, Qbrst} = 

[iQc, Qbrst] = Qbrst 
[^QG,^QG] = 

as a global residual of the gauge symmetry in Landau gauge, besides the global color sym- 



metry [12(] and the ghost-anti-ghost symmetry [13]. The important consequence of (11 11 1) is 
that the BRST charge is also nilpotent. 

A well-defined nilpotent charge splits the state space directly into three disjoint parts 



ll[. The states which are not annihilated by the BRST-transformation form a subspace 
Qi, carrying BRST-charge. By acting on these states, daughter states in a subspace Q2 
are generated which are annihilated by the BRST-charge. The last possibility are states 
which are also annihilated by the BRST-charge but are not generated from parent states. 
These form a subspace Qq. Physical states must be gauge invariant and are therefore 
annihilated by Qbrst? which is essentially a gauge transformation with the ghost field acting 
as the transformation parameter. In addition, any states in Q2 do not contribute to matrix 
elements, since they have zero norm. Therefore the physical subspace is 



ifphys = KerQBRST/ImQBRST = Qo- (112) 



101 



It is this subspace in which the perturbatively physical transverse gauge bosons exist, while 
forward polarized gluons and anti-ghosts belong to Qi and backward polarized gluons and 
ghosts belong to Q2. This can be seen directly using the Nakanishi-Lautrup formulation of 
the gauge-fixed Lagrangian 0, 12|. Due to the relation of Q2 and Qi, the unphysical degrees 
of freedom are connected by BRST transformations and are thus metric partners. They form 
a quartet under this charge, and are thus said to be confined by the quartet mechanism. 
Hence in perturbation theory the physical subspace Qo contains only transverse gluons, and 
perturbatively unphysical degrees of freedom are confinecl^. 



4-3.2. Non-perturbative BRST and the Kugo-Ojima confinement scenario 

The basic idea of the Kugo-Ojima confinement mechanism is to extend the BRST con- 
struction such that not only tr ansve rse gluons, but also in general any colored state is moved 
to the unp hysic al subspace [l2|,|288|. E. g., pairing transverse gluons with gluon-ghost bound 



states 



12, 290 



in a quartet would make them unobservable. 
The basic ingredients needed for this construction are threefold. 

One is the existence of a mass gap in the physical spectrum. Since the kno wn lowest - lying 



glueballs all have quite a significant non-zero mass, this appears to be fulfilled |274j . l29ll . |292 
For QCD this would actually be known from experimento pLi]- At the same time no mass 
gap is permitted in the complete state space, including the unphysical states. As has been 
shown in section the absence of a mass gap can be read off the ghost Schwinger function 
for all cases investigated. This is therefore also fulfilled. 

The second ingredient is that the global color charge after gauge-fixing is well defined 
and unbroken. If all other preconditions of the construction are fulfilled, this is immediately 
implied if the ghost propagator has a stronger than massless pole with negative residue, i. 



e., its dressing function must diverge to positive infinity [12|, |37j. Evidently, this is satisfied 
in the scaling case. But to the finite-ghost case the Kugo-Ojima scenario in its original form 
cannot be applied. There have been several investigations whether it is possible to construct 



a modified Kugo-Ojima mechanism also in this case [293l - l298l |. but this has not been finally 



settled. Thus, the remainder of this section only applies to the scaling case. 

The third, and probably most complicated, requirement is that there exists, in the sense 



of differential geometry, a globally well defined and unbroken BRST charge [37| ■ Moreover 



the field transformations induced by this charge must take the same form as the perturbative 



oneS fll07lll09|) [12]. 

The perturbative definition of the BRST charge is trivially not sufficient in the non- 
perturbative setting. The reason is that perturbatively a BRST transformation mediates 



^•^In principle it is pos sible to have states in Qo with non-vanishing ghost number, which would render 
the theory ill-defined 12] ■ This does not seem to be the case for Yang-Mills theories. 

^^Note that for QCD in the chiral limit with a massless pion this would have to be reconsidered. 
Actually, when the original work of Kugo and Ojima was done [12j the existence of Gribov copies was 
unknown, and the perturbative definition of the BRST charge was used. Since only the algebraic properties 
of the charge are relevant to the construction, this supplemental condition is sufficient to perform the 
construction even in presence of Gribov copies [37,]. If this condition is not fulfilled, the original Kugo-Ojima 
scenario fails. 
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between different gauge copies separated infinitesimally in a covariant gauge. It shrinks to 
the identity transformation in Landau gauge, since perturbatively only one representative of 
the gauge orbit fulfills the Land au gaug e condition. Since there are Gribov copies, connected 
by large gauge transformations 17|, |45| , a well-defined BRST transformation cannot simply 
be just the identity. This is also clear from the fact that massless transverse gluons are put 
into the physical state space by the perturbative construc tion , thus being in conflict with 



29S 



experiment. A more formal argument can be found in |64l . 

The question is then whether there exists a non-perturbative extension of the BRST 
ch arge such that its algebra coincides with the perturbative one. In fact, there is one jSTl . 



59 



ch is based on first performing a lattice regularization of the theory, and then 
taking the continuum limit. The result is that one obtains a non-perturbative BRST charge 
in the fully quantized and renormalized theory, i. e., with all renormalization constants fixed. 
This Neuberger-von Smekal construction even leads to the same algebraic form of the non- 
perturbatively well defined BRST charge. Indeed, together with the assumption of a mass 
gap in the physical spectrum this is then sufficient to imply the Kugo-Ojima construction if 
an infrared divergent ghost dressing function, and thus a globally well defined color charge, 
exists, which appears to be the case in functional calculations of the scaling case 37 . 

The remaining problem is how to relate the argument to the present situation, in partic- 
ular concerning lattice calculations. The reason is that the Neuberger-von Smekal construc- 
tion is based on an average over all Gribov copies satisfying the Landau gauge condition 
(P). This implies an aver age o ver all Gribov regions, weighted by the signed Faddeev- Popov 



determinant [57|, l59|, l60|, ll3Cll | , and is thus equivalent to a full non-perturbative evaluation 



of the path integral ([7j) in the sense of Hirschfeld 55|]. Since the DSEs and FRGEs are 
identical whether they are applied to a single Gribov region or to the whole of all Gribov 
regions 



65, 87 



the existence of the scaling solution implies that it is associated with the 
sum over Gribov regions by virtue of the Kugo-Ojima construction as outlined above 371 . 

The question is then how to relate this scaling case to lattice calculations presented here. 
If this were possible at all, there must exist some prescription to select Gribov copies inside 
the first Gribov region such that all correlation functions take the same value as in the 
Neuberger-von Smekal construction an thus be equivalent to the average over all Gribov 
copies. This may appear hopeless at first sight. However, there are a number of arguments 
which motivate that such a relation may indeed exist. 

The first is the simple fact that the scaling solution has a positive semi-definite ghost 
dressing function, both in momentum and real space. In general, outside the first Gribov 
region this does not need to be the case. An explicit counter example is given by 1+1- 
dimensional Coulomb gauge, which can be solved exactly 82|- In this case it was also shown 



that for a finite number of Gribov regions the scaling case is only realized if the restriction is 
made to the first Gribov region, when only a finite number of Gribov regions is considered. 
The second argument is that scaling appears to be realized insi de the first Gribov region 



in two dimensions, for reasons not yet fully understood [72|, |90|, l9J, Ill6l . Ill7l . Il86l . l238| . Since 



the construction of the glob al BR ST proceeds in two dimensions in the same way as in 
higher dimensions 57, 59, 60, 13d |. in this case even the minimal Landau gauge seems to be 
equivalent to a sum over all Gribov copies. 
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The third argument is that the finite-volume version of scahng can be reah zed at lea st 



for small volumes inside the first Gribov region using the Landau- -B gauges [8ll, lllll, Il45 



All of these arguments are not sufficient to ensure the possibility of scaling in the first 
Gribov region, nor can they explain how the appropriate cancellation of Gribov copies should 
occur such that this is possible. However, it is sufficient motivation to indeed investigate 
whether such a construction exists. In particular, it requires a realization of the idea of a 
simultaneous formulation on both the lattice and in the continuum, e. g. like fl26|l or fl27|) . The 
Landau- i? gauges provide at least a possibility in principle how this could be realized, and 
an explicit construction description to verify this possibility: Trace the maximal attainable 
value of 5 as a function of volume inside the first Gribov region. If it diverges in the 
infinite- volume (and continuum) limit, scaling could be realized. Check then whether the 
correlation functions exhibit a scaling behavior. Unfortunately, current lattice calculations 
can at best give evidence that such a behavior could be correct, and, even worse, could never 
falsify it. Furthermore, this does not exclude that a different selection procedure for Gribov 
copies could still yield scaling inside the first Gribov region, even if it were not realized in 
some Landau-S gauges. An alternative would, e. g., be to select Gribov copies such that 
the infrared behavior of the running coupling is as constant as possible. 

It thus remains at the current time a speculation, based on the arguments presented, that 
such a connection could exist. This working hypothesis is the central element to make contact 
between the results obtained for the correlation functions and the Kugo-Ojima construction, 
as well as the possible scaling Landau-i? gauge on the lattice and the DSE/FRG scaling 



case. This is a matter of ongoing research [62|, |87|. However, the realization of the Kugo- 



Ojima scenario is not necessary for gluon confinement, as noted above, but it would be very 
desirable, as it would immediately extend the construction to all colored quantities, and 
would immediately permit an interpretation of Yang-Mills theory as a local quantum field 
theory beyond perturbation theory [186]. 

As pointed out above, a vital ingredient is the divergence of the ghost dressing function. 
Thus, in the finite-ghost cases the Kugo-Ojima construction is not possible, and in fact even a 



local off-shell version of the standard BRST symmetry is not implementable |64j. l69l . l300l . l301 



Still, proposals exist to provide even in the finite-ghost case a non-local symmetry which 
could possibly replace BRST symmetry an d paye the way to a construction analogous to 



the Kugo-Ojima one in the scaling case [7l|, l298l . |299| . |302| - |307 



4-3.3. Gribov-Zwanziger confinement scenario 

A scenario based on a completely different point of view, but leading to rather similar 
predictions for the correlation functions, is the one of Gribov and Zwanzig er Jl3,|63. 



The original basic idea started by restricting to the first Gribov region [17|]. In a semi- 
perturbative way it was then shown that this restriction yields an infrared divergent ghost 
dressing function with a negative residue. When it was later realized that there are also 
Gribov copies inside the first Gribov region it became clear that this required further refine- 
ment. 

The next step was the assumption that it is permissible to implement a selection of Gribov 
copies inside the first Gribov region by replacing the 6 function in fl?T]) by a (5-function 67 . 
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This reduces the region of field configuration space in the functional integral to the first 
Gribov horizon. This was motivated by the idea that most Gribov copies are located there, 
since most of the volume of the hypersurface of Landau gauge in the first Gribov region 



65, 97 . From 



of the total infinite-dimensional field configuration space is located there 
the current point of view this corresponds just to a particular choice among Gribov copies, 
though without a formal justification up to now. 

Assuming this procedure to be valid, it was possible to cast this condition with the help 



of auxiliary fields into a renormalizable and local, albeit complicated, Lagrangian [67|, |299 



Already at tree-level this Lagrangian exhibits an infrared divergent ghost dressing function 
and infrared vanishing gluon propagator, in qualitative agreement wit h th e results in the 
scaling case [67|. This is not altered by perturbative corrections [73|, 1 184 l308l . |309[ | and 



remains true beyond perturbation theory 76 



This qualitative agreement even becomes quantitative beyond tree-level in the non- 
perturbative domain [74']. Implementing the same restriction of perfect tree- level cancellation 
like in the conventional scaling case in certain DSEs in the Gribov-Zwanziger framework, 
it was found that the correlation functions exhibit the same qualitative behavior. In par- 
ticular, the ghost dressing function still diverges in the infrared. Depending on the precise 
behavior of the ghost-gluon vertex, even the same infrared critical exponents are found for 
the Faddeev-Popov ghost and gluon propagators. 

Thus, so far this appears just to reproduce the scaling case. Given the results of lat- 
tice calculations in the minimal Landau gauge in three and four dimensions, it has been 
investigated whether it is possible to impleme nt a finite-ghost behavior also in the Gribov- 
Zwanziger Lagrangian 64. l69l. l70l l72L I209L l309| . This is indeed possible, by imposing conden- 



sation of certain combinations of the auxiliary fields needed for localization. An una mbigu ous 



determination of the values of the condensates is an intricate problem j6J, |70|, l79|, |310|. In 



fact, a possibility would be that the value of these gauge- dependent condensates could be 
determined by the treatment of Gribov copies, e. g. corresponding to the B parameter of 
Landau-i? finite-ghost gauges. This possibility has not yet been explored. Note however 
that dimension-two condensates, though being employed also beyond the Gri bov-Z wanziger 
scenario 



256| . may be problematic in an operator product expansion setting |311 



Irrespective of the actual values of the condensates, the generic result of p ert urbat ive 
calculations in these cases coincide qualitatively with the finite-ghost case 6J, |70|, l209| . i. 
e., an infrared photon-like ghost and a screened gluon. In addition, in these cases also 
a Schwinger function showing the same qualitative behavior for the gluon propagator is 



found [6J, l70|. Thus, the results are in good agreement with lattice results in three and 



four dimensions. Interestingly, at least the perturbative treatment in two dimensions breaks 



down [72(1, i. e. exactly for the case in which it is not yet clear whether a finite- ghost behavior 
can be achieved by a choice of Gribov copies in lattice calculations. 

If indeed the imposed 5-function on the field configuration space is equivalent to the 6- 
function, it is also possible to provide a connection to the results in absolute Landau gauge 
65j . The basic assumption is that all correlation functions of a finite number of field variables 



take the value which is most probable, based on a maximization of the corresponding en- 
tropy. This is exactly analogous to what indeed happens in the minimal Landau gauge, thus 
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indicating a possible formal connection. In this case, the results of the Gribov-Zwanziger 
Lagrangian for these correlation functions should coincide with the ones of minimal Lan- 
dau gauge. If furthermore all copies become equivalent for such finite-order polynomials of 
field operators, like discussed in section I2.5.3[ the result would also agree with the absolute 
Landau gauge on the lattice. 

4.3.4- Gribov-Stingl scenario 



A more pragmatic scenario is given by the one of Gribov and Stingl [17|, l279l - l28ll | . In this 
scenario, the existence of a double pole structure similar to f llOOp is postulated. However, 
the generic propagator must have the form 



D{GS 

where Z{p) satisfies the condition 



Zip) 



+ 2m^ cos(20)p^ + 



Z (ime^^'^) = 0, 



(113) 



(114) 



1. e.. 



the residue vanishes at the poles. Under these conditions the appearance o f physica l 



gluons can be avoided, and such solutions have been argued for using DSEs [2791 - 1281 



However, in case of the full solutions of the functional equations, as presented in section 
14. 2[ the corresponding analytic behavior is not compatible with the Gribov-Stingl type. 
Whether it is compatible wit h lat tice calculations is yet unclear, since though it fits the 
gluon propagator rather well |239j . it is not yet possible to discern whether the Schwinger 
function on the lattice shows the necessary oscillatory behavior. 

4.4. Relations to topological structures 

Many popular and powerful confinement scenarios are no t based on corr elation functions, 
but rather on topological, i. e. collective, gluon excitations I8OI . I312h314| . The se scenarios 
in tu r n ar e usually based on one particular type of ex citations, like vortices 180| . monopoles 



312, 314 



merons 



313|, and possibly others |315| . From this non-exhaustive list it is 



evident that not even the underlying dimensionality of the relevant structure is yet decided. 
Still, there has been quite a number of investigations to find relations betw een the variou s 



excitations, yielding e. g. d eep connections between vortices and monopoles |180l . I316l - l319 



even at finite temperature [320 



Given that these scenarios a re quite s uccessful in describing properties of non-local ob- 
servables like the string tension I8OI . |312| it is of importance to understand how correlation 
functions contain the same informatiorl^4 That this information is contained follows directly 



^^Of course, there is no guarantee that any kind of truncation in functional equations will preserve any 
of this. Given that the modern truncation schemes presented here maintain not only the shape of the 
correlation functions b ut al s o ef fects which are usually attributed t o topolog ical excitations, like chiral 
symmetry breaking 3, 174L 321 1, possibly the formation of a string 195L I9l\ . or the ry' mass 322 1, it is 
likely that at least a considerable amount of the topological information is captured. 
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from the fact that the exact correlation functions are sufficient to reconstruct the complete 
generating functional using the equality ( 1T7|) jij, and in turn the correlation functions are 
fully described by the functional equations by virtue of the relation (|T2D [Ss']. In particular 
cases, this has even been demonstrated explicitly for solvable models [13|, I323| . 

Besides such formal arguments, there is an explicit possibility to investigate this relation. 
This is done by obtaining the correlation functions in a topological background field and 
comparing the results to the ones obtained in the full Yang-Mills vacuum. This has been 
done using analytical and modeling methods, which indicate that the characteristic features 
of correlation functions are infiuenced or even dominated by topological excitations. In 
particular, the infrared behavior of the ghost propagator has been linked analytically via 
the Faddeev- Popov operator to topological excitations [4^ |324| , as have been the properties 
of the gluon propagator |255 . 

Complementary to analytical methods, lattice simulations provide a possibility to access 
the topological co ntent of a fie ld configuration by smearing and projecting, at least in an 



approximate way 



The so treated field configurations become self-dual, which 



180L I325L 132 

classifies them as topological [2|, |327| . Evaluating the correlation functions in such smeared 
or projected configurations provides then a possibility to assess which generic features are 
dominated by topological contributions. 



One possibility for such a smearing prescription is the so-called APE smearing |32(: 
APE smearing on the lattice is implemented by iteratively replacing links by a weighted 
average over neighboring links 



UJx) alluix) 



X) 



projected to the group 



where "projected to the group" implies that the non-group element U'^ found after addition 
is replaced by the group element U'^ closest to the result, where the distance is given by 
trU'i^U'^, with no summation implied. A typical value for a, which is also employed here, is 
0.55. Choosing a to be zero is referred to as cooling. 



Results for such a calculation are shown in figure IMl in minimal Landau gauge 32 8| . The 



absence of the ultraviolet tail of the correlation functions after a few smearing steps results 
from the removal of the ultraviolet fiuctuations by the smearing process. Furthermore, the 
behavior at or below Aym is qualitatively almost unchanged, showing that indeed the low- 
momentum behavior of correlation functions does not only contain topological information, 
but is rather dominated by i10. Simil ar re s ults have been obtained using projecti on m ethods, 
in particular for center vortices [285i, l329l . |330|, and also in Coulomb gauge 



99, 330 



This indicates that the collective excitations are indeed refiected in the properties of the 
correlation functions at low momenta. Even though more exact analytical mappings have yet 



^^This also supports that the truncations presented here captures a significant part of the topological 
information. 
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Gluon dressing function | 



Ghost dressing function | 




Figure 34: The gluo n dr essing function (left panel) and ghost dressing function (right panel) evolving under 
APE smearing 62. l328l |. Already at about 100 smearing steps, the configurations are essential self-dual 
[325| . and thus approximately pure topological configurations. Results are for a V — (2.5 fm)^ lattice at 
a = 0.21 fm and su(2) in the minimal Landau gauge. 



to be developed fully |44, l62|, |255| . |324| | , this makes clear that topology-based and correlation- 
function-based approaches to the non-pert ur bativ e properties of Yang-Mills theory are just 



two (equivalent) facets of the same physics [289|. Therefore, also the question whether a 



description in terms of elementary particles or of collective excitations is more appropriate 
is rather a question of practicality for a given question than a conceptual difference. 

As a side-remark, it is also interesting to investigate whether the correlation functions 
depend on the net topological charge. This is exemplified in figure |35l where it is seen that 
the dependence is probably only quantitative. This is of considerable practical importance, 
given that topological properties have a much longer relaxation time in lattice simulations 



than most other quantities |331 



4-5. Beyond Yang-Mills theory 

Of course, the approach discussed here is not limited to gauge bosons in Yang-Mills 
theory. Two rather important extensions are the inclusion of matter fields and the deter- 
mination of the properties of composite objects, in particular bound states. However, these 
topics are beyond the scope of this review, and will only be briefly introduced here. 

4.5.1. Matter fields 

Adding matter fields to Yang-Mills theory is rather straight-forward |2|1. Since Landau 
gauge, in contrast to e. g. 't Hooft gauge outside the Landau gauge limit |2|, |2J], does not 
involve the matter fields in any active way, all that has been said so far on gauge-fixing 
remains valid. This does not imply that the process remains completely unchanged. E. g., 
the corridor of the permitted B values may change in the presence of matter fields. An 
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Gluon dressing function | 



Ghost dressing function | 




Figure 35: The gluon dressing function (left pane l) and ghost dressing function (right panel) as a function of 
topological charge after ten APE sweeps [H^, l328l | . The topological charge has been measured as the volume 
integral of -F^j^-F^^ [273[. Results are for a V = (1.7 fm)^ lattice at a = 0.21 fm and su(2) in minimal Landau 
gauge. 



indirect example of this has aheady been foun d in the Higgs phast 
Higgs- Yang-Mills system in both Coulomb 



99, 319 and Landau 



^1 of a fundamental- 
gauge. In this case 



the average properties of the Faddeev- Popov operator change qualitatively compared to the 
confinement phase. Though this has not yet been investigated in detail, this makes it likely 
that also the corresponding 5-corridor is different in both phases [l^]. 



The same goes for the viability of the truncations in functional calculations [l60l . Il98 



not 

to mention technical problems in lattice calculations when including fermions |43|]. However, 
these are technical rather than conceptual problems. 

Thus, there have been plenty investigations including matter fields. The emphasis is of 
course on fermions in t he fu r idamenta l represen tation, as these are necessary for an investi- 
gation of QCD [l3|, [ij, Il95l . 12171 . l32ll . I333l4340[ ^. Less o f ten the properties of adjoint and 



fundamental scalars have been investigated [2J, llGOl . Il66l . Il97l Il98l . |342| . l343l | , though they 
are much more amendable to lattice calculations j^. Investigations of fermions in different 
representations than the fundamental, however, have received significant attention lately 
[oj, as these are candidates for theories exhibiting an approximate conformal behavior in 
the infrared, and thus are possibly relevant for, e. g., walking technicolor extensions of the 
standard model. Howe ver, determinations of the corresponding correlation functions are 
just beginning 1344 . 



^^Strictly speaki ng, the H iggs phase is continuously connected to the confinement phase when a (lattice) 
cut-off is imposed 319l . 332 1. but the notation will be kept for the sake of simplicity. 

^^This only includes references which to some extent take into account the Yang-Mills correlation functions 
and are in Landau gauge and only focusing on the properties of the quarks. There is a wealth of literature on 
calculations using various ansatze for these correlation functions, see [isl l264l l34lj for introductory reviews. 
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Figure 36: Examples of the gluon propagator and ghost dressing function for theories wit h different ma tter 
content. Top- left panel: Two flavors of fundamental fermions and gauge group su(3) [liillHilMl on 
a 24^ X 48 lattice and Wilson clover fermions with k ~ 0.13575. The scale is fixed by the intermediate- 
distance string tension. Top-right panel: Two flavors of adjoint fermions and gauge group SU(2) [2^ l346ll 
on an 12'^ x 24 lattice with Wilson fermions and —amo — 0.95. Scale is fixed by a 2 TeV techniglueball [25[- 
Bottom-left panel: The confinement phase for a doublet of complex fundamental scalars for SU(2) in the 
confinement phase 2J| on a 24^ lattice with k = 0.25 and A = 0.5. Bottom-right panel: The Higgs phase for 
a doublet of complex fundamental scalars for SU(2) in the Higgs phase 2j| on a 24^ lattice with k = 0.32 
and A = 1. In both cases the scale is fixed by a 250 GeV Higgson ium l 24ll . All results are in minimal Landau 



gauge. Results using functional methods can be found, e. g., in jl60l . Il95lll97l . l321 | 



A sample of results for propagators for various cases are assembled in figure 123 Going 
into details or even listing a complete set of references just for Landau gauge is far beyond 
the scope of this review. It suffices to say that one can investigate fro m them a phletora o f 



phy s ical phen o meria like , e. g., the dynamical generation of mass |14J . l258l . I259L I266L l335 



8881 . 1844 18451 . I847L l848 



chiral symmetry breaking and restoration at finite temperature 
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2661 . l349l . I350L the de confinement transition [23 71 . I350l - l353| . center observ ables like the 
Polyakov loop |237l. l354| , conceptu al access to the string tension |l95l - ll97l . l355l |. co nfinemen t 
ofmatter fields |l95 - 197 . 354 . 356 1, color-superconducting phase s at l arge densities 227, 357 - 



m, the H-ggsSTctlii^lild W;ared couformal.ty BS 

However, already from the brief presentation in figure 



36ll |362| |. only to mention some. 



it can be deduced that the presence 
of matter fields can be both, almost irrelevant and of qualitative impact for the properties 
of the Yang-Mills sector. 

This rapidly growing field shows the versatility of correlation-functions-based methods. 
The wealth of possible physics results has yet only been tapped to a small extent, and will 
certainly be the focus of many investigations in the time to come. 

4-5.2. Composite objects 

The natural quantities to be determined from the point of view of experiments, beside 
the ultraviolet properties of elementary particles, are the properties of bound states. In par- 
ticular, quantities like masses and decay constants, but also widths, have therefore received 



a lot of attention 13|, ud, Um Ml 



In principle, these quantities can be obtained from the corres ponding c orrel ation func- 



tions, e. g. a glueball mass from the 4-point correlation function [269|, l271l . I273j . Since the 



masses of colorless bound states are gauge- invariant and renormali zatio n-group invariant, 
they can even be obtained without explicit gauge fixing ji], [Hi 43, 273 1. In particular for 
lattice calculations use has been made of this possibility. Despite the numerical costs to 
simulate quarks with physical r nasses it was even possible to obtain results essential at 
the physical parameters of QCD 363|-|365|. Using this approach, information on the gauge- 



dependent internal structure of composite object s is lost, a nd at best it is possible to expand 



3661 . 1367j . In particular, relative momenta 



the states in terms of gauge-invariant operators 
between the compounds are not accessible without fixing a gauge. 

Similarly, it is also possible to use functional methods to determine the gauge-invariant 
masses of colorless bound states. Instead of using the DSEs or FRGs for the corresponding 
correlation fun ctions dir ectly, it is possible to employ the Bethe-Salpeter 2171 . l333l l334j . l368| 



and Faddeev |265L l369l | equations to obtain the mass spectrum of mesons and baryons 



131. Il4 1264 1341 



respectively. This has been widely discussed in the literature 

The lattice results suffer primarily from statistical and systematic uncertainties, the 
latter mostly due to the extrapolation to the continuum and infinite volume, as well as 
to physical quark masses. The functional calculations suffer in contrast from the errors 
introduced by the truncations. Thus, comparing both approaches, the uncertainties have 
completely diffe rent origi ns. Thus, once more, the comparison of both methods can give 



mutual support [230|, |370|. In fact, the results of both methods are in acceptable agreement 
with each other and with the experimental values, thus giving even more confidence to have 
access to a theory like QCD even at low momenta and beyond perturbation theory than a 
single method alone could. 

It should be noted that the results using a hierarchy of correlation functions, as with 
Bethe-Salpeter equations, are rather insensi tive to the details of the infrared behavior of 
the gauge- dependent correlation functions |333l . l339l . |370| - |372| . This is as it should be. 
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given that this infrared behavior is essentially dominated by the gauge-fixing procedure. 
Explicit calculations show that, when taking the gauge-dependence of all involved correlation 
functions into account, the results for the finite- ghost case and the scaling case, and thus 
ossibly for different non-perturbative gauges, agree within the errors, and with experiment 



333| . This also provides the possibility to understand how gauge- invariance is recovered for 



176 



physical states, as has been done previously for QED 

This concludes the construction of the framework at zero temperature, which gives an 
impression how to start with the elementary degrees of freedom and obtain the correlation 
functions. A discussion more detailed than the short presentation here how to finally arrive 
with this approach at gauge- invari ant, qbseryable quantities can be found in various reviews 



and textbooks, e. g. [9|, ll3|, M, |43|, UM, UM Ml 



5. Describing gluons at finite temperature 

In the previous chapters the properties of gluons have been determined using correlation 
functions. It has long be en d iscussed whether these properties fundamentally change when 



the system is heated up [373| . In the case of Yang-Mills theory it is rather well established 
that a phase transition oc curs at a n on-zero temperature Tc of th e order of Aym, for any gaug e 
group investigated so f arjsii, H , including SU(iV < 9) Sp(2) and Sp(3) |376l.l377i. 



E7 [377[, and G2 |378N38nl . Only the order and other details of the phase transition turn 



373, 381 . Based on 



out to be dependent on the particular gauge group. 

One of the central ideas about this transition was deconfinement 
the argument that the running coupling evaluated at the temperature becomes small for 
increasing temperature it was assumed that confinement as a strong coupling phenomenon 
should cease. At first sight, this is supported by the fact that thermodynamical quantities, 
like the pres sure j loga rithmically approach a Stefan-Boltzmann behavior with increasing 
temperature 



374, 382 



It is of course clear that this can only be superficially correct. On the one hand, Haag's 



theorem forbids the absence of non-perturbative interactions [10|. On the other hand, this 
would require a gauge-invariant formulation of asymptotic, almost local colored states at 
finite temperature. Though the existence of such a formulation has not yet been excluded, 
the impossibility of such a bleaching at zero temperature for non-Abelian gauge theories 
Tol makes this rather unlikely. Furthermore, the impossibility of stable asymptotic states 



of elementary particles at finite temperature due to the Narnhofer-Thirring theorem |383 
alleviates this question anyway [384:|. 

After introducing in the following the formalism to describe gluons using correlation func- 
tions at finite temperature, it will be shown that indeed gluons are not deconfined. In fact, 
they do not develop an asymptotic physical pole, as will be read off from their Schwinger 
functions once more. However, it will be necessary to distinguish the possibilities of trans- 
versely and longitudinally polarized gluons with respect to the heat bath four-momentum. 
It will furthermore be argued that persisting gluon confinement is not in contradiction to a 
Stefan-Boltzmann-like behavior of thermodynamical bulk quantities. These will be found to 
be dominated by hard processes, i. e., at energies of the order of the temperature. Since this 
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becomes a hard scale with increasing temperature, indeed their leading behavior will begen- 
erated by leading-order perturbative contributions. However, the Linde problem [23|, l373l | of 
such a perturbative treatment will be cured by the non-perturbative effects. In addition, the 
formal infinite-temperature limit will be discussed, showing the validity of the arguments at 
arbitrarily high temperatures. 

Finally, this will be concluded by showing that already from the two-point functions 
gauge-invariant quantities like the critical temperature and order parameters can possibly 
be extracted. 

It should be noted that the approach has also been applied to the case of finite density, 
after coupling to fundamental fermions. Since this goes beyond Ya ng-Mills theory, this wil l 



not be described here. Discussions of this topic can be found, e. g., in {2271 . l34ll . l357H360l . |385 



5.1. Finite temperature 

In the following only a thermostatic setting. 



1. e. 



equilibrium, will be discu ssed |373 



386 



is not yet 



The extension of the framework presented here to non-equilibrium states 
complet ely co nstructed, especially beyond perturbation theory, but progress in this direction 
is made 387 . 



5.1.1. Formulation 

An adequate setting for equilibrium calculations is the Matsubara formalism |373j . With- 
out going into detail, it essentially requires to change to Euclidean spac^ and then com- 
pactify the time direction with size = 1/T. Therefore, equilibrium is manifest as no 
genuine time dependency remains. 

A further consequence of this distinction of the time direction is that energy and three 
momentum become independent variable s. T hough this may at first seem like a breaking 
of Lorentz invariance, this is not the case (388 1. In this standard setting, merely a choice 



for the four velocity of the heat bath is made such that the latter is at rest. By explicitly 
including this four velocity in the calculations, Lorentz symmetry can be made manifest once 
more. Despite the gain in symmetry, the practical calculations become rather cumbersome, 
and it is therefore simpler to have only hidden Lorentz symmetry. This will be done here 
as well. However, it should be kept in mind that this is just a choice of frame, without 
importance for the physics. A further consequence is that at four-momenta with ^ 
temperature becomes negligible, and Lorentz symmetry, or in Euclidean space-time 0(4) 
symmetry, becomes restored: If all momenta are large compared to the temperature the 
(perturbative) four- dimensional behavior of correlation functions w ill be rec overed. That 



164j . |389|, and applies 



is essentially a consequence of the Appelquist-Carrazone theorem 
within the same domain of validity. 

Furthermore, the reference direction of the heat bath makes it possible to distinguish 
between a polarization in the direction of the heat bath or transverse with respect to it. As 



In a sense, the finite temperature results are therefore much more closer to the real world than the 
results presented at zero temperature where analytic continuation or the use of the reconstruction theorem 
are necessary to make contact with Minkowski space. 
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a consequence it is necessary to decompose the gluon propagator as 
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D;t{po,p) = P^APo,P)D^T\po,f) + P^Apo,p)Dl\po,f), (115) 
i. e., there are two independent dressing functions. The tensor structures are given by 



iPo,P) 



l-5^o)(l-5, 



lyO) 



PL{Po,P) 



pM-p^Ap)- 



A third possible tensor structure, having a mixture of longitudinal and transverse indices, 
is due to the requirement of Landau gauge that Dui/ is four- dimensional transverse not 



independent, and thus doesnot appear explicitly, 
propagators are given by [216 



In terms of correlation functions, these 



Drip) =^ 



1 




pI 



Dl{p) 



MM. 

„2 



{d-2)Ng 



a';{p)a;{-p)-^a-,{p)a-,{-p) 



(116) 
(117) 



where explicit use has been made of the Landau gauge condition. Of course, the ghost 
propagator as being a scalar is still represented by a single scalar function, though now also 
depending separately on energy and three-momentum. Furthermore, due to CP invariance, 
the correlation functions do not change under the replacement po into —po- Thus, only 
positive (and zero) energy solutions will be discussed below. 

As a consequence of a compactified time direction the energies of the system become 
quantized. Because of the Kugo-Martin-Schwinger condition gluons as bosons have to have 

i )erio dic boundary conditions, while fermions would have anti-periodic boundary conditions 
390|. Thus, the energies gluons can have are 



Po 



27rTra, 



:il8) 



with n an integer. To maintain gauge invariance, it is then necessary that ghosts, despite 
their Grassmannian nature, also have periodic boundary conditions, and thus the same 
energy levels 39l|. Since for n = the energy vanishes there is a natural separation of 
states in soft states n = and hard states n ^ 0. The latter can never have p^ = 0, and thus 
cannot be on the tree-level light cone. In a perturbative setting only the soft modes can 



go on-shell, though non-perturbatively this statement loses its meaning |384j |. Furthermore, 
from the perturbative point of view, there are then three different energy scales. One is given 
by the hard scale T, while there is a separation of soft states g(T)T and ultrasoft states 
g{T)'^T. Here, this further distinction will not be made. A more important distinction is to 
declare a momentum to be infrared if and only if it is not only much smaller than Aym, but 



in addition also much smaller than T, as already given in fl63|l |216 
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With this distinction comes another insight. When the hmit of infinite temperature 
is takeiJ^. which will be discussed in det ail in sect ion I5.2.3[ only the soft modes have 



finite energy and are therefore dynamical |l94l l392| . Thus, it is formally admissible to 



remove all hard modes from the Lagrangian. The result is the Lagrangian of a three- 
dimensional Yang-Mills theory coupled to an adjoint scalar. Formally, the ghost remains 
unaltered in the process, and the spatial components of the gluon field become the three- 
dimensional gauge field, while the zero component of the gluon field becomes the adjoint 
scalar. From the relations f lll6|l and (11171) this is immediately clear, as the soft mode of 
the transverse propagator is just the corresponding spatial gluon propagator and the adjoint 
scalar propagator is given by the soft longitudinal propagator. However, in contrast to an 
arbitrary such theory, the mass and self-coupling of the adjoint scal ar ar e not free parameters. 



but are fixed by their descent from the four-dimensional theory [393|. In particular, since 
the parent theory is only logarithmically divergent the mass of the adjoint scalar is protected 
from large radiative corrections of the order of the cutoff with the four- dimensional gauge 
symmetry acting as a custodial symmetry. The infinite-temperature limit can therefore 
be treated in such a reduced three-dimensional setting, which will be made use of below. 
Since the longitudinal propagator is then entirely given by chromoelectric components of 
the gluon field and the transverse part entirely by chromomagnetic ones, these notations 
are also used at finite temperature. Since the ghost in the limit of infinite temperature 
becomes the one of the three-dimensional Yang-Mills theory, it is usually also associated 
with the chromomagnetic gluon, and both together are referred to as the magnetic sector of 
the theory. 

5.1.2. Modifications of the methods 

After these general consequences of introducing finite temperature, the necessary alter- 
ations to implement the Matsubara formalism in the methods employed should be noted. 

In lattice calculations it is str aightforw ard to implement finite temperature by explicitly 



compactifying the time direction |43l . llSTj . For that purpose, the time extension aNt = 1/T 
must be made much smaller than the spatial extensions aNg = L. The reason is that 
otherwise the time direction cannot be considered compactified. Such a difference can be 
achieved by two possibilities. One is that the discretization in time direction is made different 
than in the spatial direction, <^ a*. The other is that the number of lattice points in 
time direction is made much smaller than in space direction, Nt -C Ng (or a combination of 
both). While the prior possibility provides less violations of rotational symmetry, the latter 
is much cheaper. Given the necessity to reach large volumes, it is the method of choice 
here. Note that in the infinite-volume and continuum limit also at must go to zero and Nt 
to infinity, but while the product agNg has to diverge, in this limit the product atNt must 
be kept fixed at 1/T. Other than this, no further modifications in the lattice methods are 
necessary. 

The implementation inside functional methods is also comparatively simple. Besides 
including the general structure of propagators and arguments explicitly, it amounts to re- 



This limit has to be taken after renormalization 221|. 
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placing all integrals over energies with discrete sums over Matsubara frequencies 
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27r 



(119) 



where n is the integer in f lllSp . This is sufficient. The evaluation can be simplified by 
making use of the fact that all correlation functions remain unchanged under the replacement 
Po -Po- 

Due to the discretization, the correlation functions can be regarded as an infinite tower of 
correlation functions enumerated by n. Thus, in principle, an infinite number of independent 
dressing functions have to be determined. In practice, as shown below, it turns out that the 
approximation 

T{po,f)^T{0,pl+f) (120) 

is even for n = 1 rather well fulfilled ^1^, H lij?) . This is a consequence of the effective 
restoration of Lorentz symmetry at large momenta. Therefore, only a small number of 
Matsubara frequencies have to be treated independently, permitting the development of 
efficient algorithms 226 . 

It should be noted that all of this and the following can also be transferred qualitatively 
unaltered to dimensions different than four. However, du e to t he im mediate practical ap- 
plication to heavy-ion experiments and the early universe |373|, |390|, the presentation here 



is restricted to four dimensions. Some investigations concerning the three-dimensional case 
in the present framework are available in and show no qualitative difference to the 

four- dimensional case. 



5.2. Propagators 

As discussed previously, due to f ll20p . in the ultraviolet at momenta large compared to 
both the temperature and Aym, just the zero temperature perturbative behavior is obtained. 
In case of functional methods the truncation has to respec t this. An explicit calculation of 



how this restoration takes place can be found in [19J, |216 



In case the temperature is comparable or larger than 27rAYM; it is possible to describe 
most of the temperature effects quite successfully with finite-temperature perturbation the- 



ory 
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This is the domain of hard-thermal loop calculations [395], which are included 
in both lattice and functional methods by construction. Both latter approaches, by intro- 
duction of non-perturbative effects, also solve the Linde problem. This problem essentially 
states that from order onwards all orders of pe rturbation theory contribute equally, lead- 
ing to a breakdown of perturbation theory 373| even if the factorial growth leading to its 
ultimate breakdown [33|] would not yet be a problem. The reason is that in perturbation 
theory the gluons acquire a screening mass, which is of order gT in case of the longitudinal 
propagator. By resummation, this produces an inverse dependence on the coupling con- 
stant, yielding by power-counting the breakdown of perturbation theory since powers of g 
in numerator and denominator cancel to one. A further obstacle is that the cancellation 
of infrared divergences of zero temperature perturbation theory breaks down due to the 
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screening of longitudinal gluons, which then cannot cancel infrared divergences generated 
by perturbatively unscreened transverse gluons. 

To overcome these problems constructively, improved perturbative methods have been 
developed 382|, relying essentially on one non-perturba tive input parameter, which can, e. 
g., be determined on the lattice, at least in principle 
introducing an effective transverse gluon screening mass 395 



396| . This essentially amounts to 
yielding the missing contri- 
bution of non-perturbative and higher order contributions. Instrumental in this process is 
the possibility t o ma p the theory at large temperatures on a three-dimensional theory, as 
discussed above [393|. 

Here, the infrared suppression of the transverse gluon propagator will be generated with- 
out an external input. Thus the dangerous infrared divergences no longer arise. This genuine 
non-perturbative effect is obtained in the non-perturbative methods, which therefore do not 
suffer from the breakdown of the perturbative series at order . Thus, the Linde problem 
is absent from the calculations presented here 194 . 



5.2.1. Infrared 

The infrared behavior is most transparent once more in the scaling case. The situation 
in the finite-ghost case will be discussed afterwards. It should be noted that due to (E3]) 
the momenta are infrared only when they are much smaller than the temperature. On the 
other hand, at momenta much larger than the temperature the zero-temperature behavior 
will remain. In particular, if a momentum window exists with momenta satisfying 



< p < A 



YM, 



:i2ii 



i. e., at very small temperatures, the temperature acts as an infrared cut-off, and in the 
same way as at large finite volume 145| within this momentum domain the four-dimensional 
infrared behavior, and in particular a power-law dependence of the correlation functions, will 



be obtained 216 



The derivation of the finite-temperature DSEs can be found elsewhere |l66l . |221| . The 
important point is that by an appropriate projection, two coupl ed eq uations for the dressing 
functions of the transverse and longitudinal part are obtained 166| . Keeping for now only 
the contributions from the ghost loop, the equations for the soft modes are given by 



G{0,k) 



1 

Z{0,k) 
1 

H{0, k) 



dqdO 



(27r)3 

v4t(0, go, k, q)G{qo, q)Z{qo, q-k) + Al{0, go, k, q)G{qo, q)H{qo, q-k) 



(122) 



(2vr) 



90 



^ / dqdeR{0,qo,k,q)G{qo,q)G{qo,q + k) 



<?0 



^ / dqd9P{0,qo,k,q)G{qo,q)G{qo,q + k). 



(123) 
(124) 



The case of the hard modes will be discussed below. The wave-function renormalization 
constants for the longitudinal and transverse propagator can become different at finite 
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temperat ure, depending on the renormalization scheme, and are thus denoted differently 
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3901 1 ■ The i nteg ral kernels At, Al, R, P for the case of a undressed ghost-gluon vertex 



can be found in |22l| . By only keeping the go = contribution in the Matsubara sums these 



equations become equivalent to the ones of the three-dimensional Yang-Mills-adjoint-Higgs 



system in the same truncation |166| . In that limit the longitudinal equation decouples, the 
longitudinal propagator becomes tree-level-like, and the transverse and ghost sector behave 
as a three-dimensional Yang-Mills theory. Thu s, in the scaling case they will exhibit the 



characteristic three-dimensional exponents |166| . On the other hand, in the limit of zero 



temperat ure i t is possible to show that these equations yield again the four- dimensional 
behavior 
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Equations ( I122til24p can also be written as 



1 



G{0,k) 
1 

Z(0, k) 
1 

H{0, k) 



Z:, + UG{0,k) + J2^G{qo,k) 



Z3T + nz(o,fc) + ^nz(go,fc) 

qo^o 

ZsL + ^UHiqo,k), 



(125) 
(126) 
(127) 



where Ilj denotes the various self-energies. As can be gleaned from the infinite-temperature 
limit in the longitudinal equation, the zero-component g'o = explicitly vanishes, since 
P(0,0,fc, = jl66j |. This is a direct consequence of the tensor structure of the ghost- 
gluon vertex, and it would require a highly non-trivial, temperature-dependent dressing to 
alter this behavior. As a consequence, the result for the longitudinal mode H{0, k) is entirely 
determined by the hard modes and depends only implicitly on the soft ones. 

This requires to consider the hard modes, which have an effective tree-level mass of 
27rTn. Since the external momentum is much smaller than this mass, the hard-mode dressing 
functions become essentially constant in the infrared. Thus, for /c — )■ 0, the dressing functions 
are, respectively, given by constants Az{ko), Ah^ko) and Ag{ko), which are assumed to be 
bounded as a function of fco- Actua lly, because of asymptotic freedom, they will decrease 
logarithmically with k^ for k^ ^ oo [216|. 

For the soft modes power-law ansatze 



G{0,k) 
Z(0, k) 
H{0, k) 



Bzk^' 



(128) 
(129) 
(130) 



are made with independent exponents k, t and / and constant coefficients Eg, and Bh- 
The scaling condition requires k < 0. Then, th e sof t self-energies 11^(0, k) and 112(0, k) take 
exactly the same form as in three dimensions |l66l |. i. e. 

2i-4Kp(^2 + 2k) 



nG(o,fc) 



-BgBzQ TCa 



k{?> + A{-2 + k)k)V{2k + 1) 
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y 



nz(0,A:) - -B^gTC^ .^(1 + (2. + |) ^ ' ^'^'^ 

where y = k"^ = k"^. These are exactly the resuhs obtained in the three-dimensional case 
(9^ , and thus will yield a three-dimensional behavior as the leading infrared behavior for the 
transverse gluon and ghost if they are not exactly canceled or dominated by the remaining 
Matsubara sums. Setting the hard- mode dressing functions equal to the constants Ag{ko), 
Azi^ko) and Ah^ko), the self-energies in the ghost equation fll25p can be rewritten as 



IlG{qo^O,k) = q'TCa (132) 




The first term is logarithmically divergent and can be absorbed in the wave-function renor- 
malization constant. The second term is sub-dominant when compared to the 3d-term and 
is fin ite after summation over go- Finally, the function Tr{k) vanishes identically as — )■ 



216| . Thus, the leading part of the ghost equation is the same as in the 3d-case, i.e. it is 
given by UG{0,k). 

For the transverse equation it is necessary to first get rid of the spurious divergences 
in both gluon equations. After subtraction of these divergences, the contributions of th e 



hard-modes in the self-energies vanish as /c — )■ 0, separately for each Matsubara term |221 
Thus, the subtracted part cannot contribute to the IR behavior in the transverse equation. 
On the other hand, the infrared cont ributions from the unsubtracted self-energies cannot 



be neglected in general |174 Il75l . Il90l | . This requires a regularization of the spurious diver- 
gence£j. To achieve this, one can replace the approximately constant dressing functions of 
the hard modes by Ag{qo){q'^ + q^)''^ , which are suppressed when r > 0. This is the pre- 
scription commonly used for regularizing the divergences at zero temperature with r = —k, 



94 1175L II86 



Then, the integrals can be performed and yield in the limit — )■ 

'TCa (C - 3)r (r - \) 
P 327r3/2r(2 + 2r) 



nf(^o,A^^o) = - ^t""' %J^L }1 (^--^)"" E -4.H>r-" (133) 



where n = go/2vrT. These sums diverge for r < 1/2, due to the term ^q"^^. For r = 1/2, 
this exponent becomes equal to —1, i.e. the sums are logarithmically divergent. Finally, 
for r > 1/2, the sums are finite and can be resummed analytically. One should also note 
that the terms in fll33til34|) behave like 1/fc^, i.e. like a mass term in the IR limit. In the 



^^Note that in the transverse case, divergences appear for each hard mod e, while in the longitudinal case 
only the sum over the hard modes is affected by spurious divergences |216l |. 
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transverse case - but not in the longitudinal case |221l | - this term may not be renormalized, 
as this is not allowed by gauge-invariance jsQoj. If either a divergence not stronger than 
logarithmic should be encountered or the same exponents r for the transverse and for the 
longitudinal case should be used, then the only possibilit}|l§ is to have r > 1/2. F or al l such 



values of r, the contribution fll33p is sub-leading in the transverse equation ( 11261) |216l . |397 . 
At the same time, the screening mass in the longitudinal equation would then be solely due 
to the regularized contribution fll34p . This result is a consequence of the truncation scheme. 
In fact, due to the decoupling of the soft modes in the longitudinal equation (see equation 
(11271) ). the result is dominated by the hard modes, which are very sensitive to truncation 
artifacts, since they live on a scale that is effectively mid-momentum. Thus, this truncation 
scheme is not able to yield a consistent description of the chromoelectric screening mass, 
and a determination of its value is not possiblaH 

More sophisticated truncations which eliminate the s purio us divergences already at zero 



temperature, e. g. the one discussed in section [323] or of |223l |. would also eliminate the spu- 
rious mixing of physics and truncation artifacts at finite temperature [s^], and thus permit 
to access the electric screening mass. Indeed, the lattice calculations shown below show that 
the phase transition dynamics is encoded in this screening mass, making its self-consistent 
determination the most important task at finite temperature. The question to which ex- 
tent the modeling of the vertices in turn determines the thermodynamic properties then 
becomes a central question. Thus, taking the vertices self-consistently into account appears 
to become mandatory at finite temperature. Ultimately, this should resolve whether it is a 
change of gluon properties or of bound state properties which drives the phase transition 
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Nonetheless, the primary message is that the longitudinal gluon propagator exhibits a 
non-zero screening mass, giving a qualitatively correct description of the physics involved. It 
is also this screening mass which prevents the longitudinal gluon propagator from modifying 
the infrared behavior of either the transverse gluon propagator or of the ghost propagator. 
Therefore these two propagators should behave at all non-zero temperatures exactly as in 
the three-dimensional case at momenta much smaller than the temperature T and than 
Aym- 

A similar discontinuous change from zero to non-zero temperature is also found in 
background-gauge calculations, usi ng r enormalization group equations and a scaling as- 



sumption for the running coupling [399|. It is thus possibly a more generic feature of the 
scaling case. 

The situation for the finite-ghost case is even more involved. In this case, the screening 
masses obtained at zero temperature now obtain temperature-dependent contributions also 
for the transverse gluon propagator. Similar modifications apply to the finite dressing func- 
tion of the ghost. In close correspondence to the longitudinal gluon propagator in the scaling 



'''^Of course, r = 1/2 just corresponds to the limiting scaling case of an infrared finite gluon propagator 
and an infrared divergent ghost dressing function. This possibility will not be discussed here for simplicity, 
as this does not add anything new compared to the discussion in section [3.2.41 

Another, though somewhat simplified, truncation has been investigated in [398j . and yielded a qualita- 
tively similar result. 
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case the transverse propagator's screening mass will now also mix truncation artifacts and 
temperature effects in a generic truncation. As a consequence, in the finite-ghost case all 
screening masses will receive a temperature-dependent contribution, but it will require a 
better understanding of truncation artifacts to separate the spurious and physical contribu- 
tions. However, this also implies that in the finite-ghost case no qualitative change occurs 
when a non-zero temperature is introduced. 

5.2.2. Intermediate momenta and temperatures 

The investigation using functional methods provides a qualitative description of the gen- 
eral infrared behavior for the propagators. In particular, they provide direct access to the 
properties related to the asymptotic state space in section 15.2.41 However, the fact that the 
longitudinal gluon propagator is sensitive to contributions on the scale of the temperature 
already indicates that a pure infrared analysis is not sufficient to describe the thermody- 
namics. That no sign of the well-known phase transition has been found is another serious 
drawback. 

A treatment of the DSEs over the total momentum range has been done to also investigate 



the intermediate momenta [1661 . |216| . 122 1| . However, since the truncation employed could 
remove the quadratic divergences only up to a finite truncation artifact, this provided another 
source for an electric screening mass: In this case the spurious divergences are formally 



under control [22 1| . since at finite temperature a mass renormalization of the zero mode 
of the chromoelectric gluon propagator is actually admissible. But this renormalization is 
in general arbitrary. Therefore, it is not possible to obtain physically useful results from 
such an approach. It would require at least a truncation scheme, like the one discussed in 
section I3.2.5[ which removes the quadratic divergences already at zero temperature^, to 
give at least the possibility for a correct behavior. However, this is a necessary requirement, 
but not necessarily a sufficient one. In particular, if, e. g., the thermodynamics were driven 
by glueball dynamics, this would possibly not be captured, except for an ingenious choice 
of truncation. It remains to solve the DSEs at finite temperature with such an improved 
truncation to clarify whether more is required. Calculations using FRGs [223i] permitted 
to have a better control over these artifacts, but are not yet in fully satisfactory agreement 
at temperatures around the phase transition. Nevertheless, they constitute a major step 
towards control of the truncation artifacts. An alternative and rather successful approach 
was to use gauges which directly include the Polyakov-loop dynamics, like Polyakov gauges 
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Fortunately, lattice calculations do not have similar problems |43|, Il37l | , and the calcu- 
lations performed at zero temperature can also be performed at finite temperature. The 
resulting zero-modes of the chromoelectric, chromomagnetic, and ghost propagators as a 
function of temperature are shown in figure [33 To show the adequacy of the approximation 
fll20|) . some examples for the hard modes are given in figure 1551 These results have been 
obtained in the minimal Landau gauge. The volumes yet accessed are rather small, so a 



significant gauge-dependence would naively only be expected for the ghost propagator |237 . 



There are also alternative truncation schemes proposed for that purpose [2021 
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Temperature dependence of the magnetic gluon propagator for su(2) Temperature dependence of the magnetic gluon propagator for su(3) 




Temperature dependence of the electric gluon propagator for su(2) Temperature dependence of the electric gluon propagator for su(3) 




Temperature dependence of the ghost dressing function for su(2) Temperature dependence of the ghost dressing function for su(3) 




Figure 37: The transverse (top panels) and longitudinal (middle panels) gluon propagat or and the ghost 



3pi 

dressing function (bottom panels) as a function of temperature and momentum 2371 1394{ . The left panels 
show results for the gauge algebra su(2) and the right panels for su(3). Independent determinations o f the 
critical t emperature Tc ca n be found in 
40i - l406t and for su(3) in [407l| . 



400|,|401|. Further results can be found for su(2) in ^12 & ,2161 . |342 . 



122 



Transverse gluon propagator for su(2) 

> 
(5 



Longitudinal gluon propagator for su(2) 



0.5 



• 




Matsubara frequency 
• P„=0 GeV 
o Po=2t T 


CM 

>1.5 

<u 

C5 


-* 




• 






1 






_ 












o 








• 




• 

o 












• o 



• o 






0.5 
n 


- o 
•o 

o 

• o 

\ 





5 

p [GeV] 



5 

p [GeV] 



Figure 38: Comparison of the soft mode to the first two higher Matsubara frequencies at T = Tc for su(2) 
on a 4 X 46"^ latt ice. In the left panel the transverse gluon propagator is shown, and in the right panel the 
longitudinal one 237[. All results are shown as a function of the four momentum = Pq + [pf' ■ 



However, the comparison to absolute Landau gauge |405l . |406[ | also shows an influence on 
the gluon propagators , though t his may be int ertwined with severe lattice artifacts observed 



at finite temperature [216|, M, |406|, |408|, |409 



It is visible that the only significant temperature dependence is seen for the electric 
gluon propagator. The temperature dependence of the magnetic gluon propagator is rather 
weak, and in particular not qualitatively affected by the phase transition. For the ghost 
propagator, besides effects which can be easily due to the different spatial volumes at the 

Sinc e the relation ( I120p is found to 
none of the hard modes show 
not surprisingly given their 



216, 237 



different temperatures, no effect is seen 
hold rather accurately even at the phase transition [23'i 
a pronounced dependence on the phase transition [62|, l23'i 
effective energy of at least 2tiTc ~ 2 GeV. Thus, the soft and hard magnetic sectors depend 
rather smoot hly o n te mperature. Also, the hard electric sector does not show a significant 



dependence |216l . |237| | . This also implies that the electric- magnetic asymmetry observed 



410l . |411| in dime nsion -two gluon condensates |412l | at finite temperature is driven just by 



the electric sector 237 



Only the soft electric part is then left. In agreement with the argumentation based on the 
infrared analysis in section 15.2. ![ there is a significant imprint of the phase transition visible 
for the chromoelectric propagator. In particular, its screening mass is severely affected, as is 
shown in figure [221 The phase transition point can indeed be determined from the maximum 
of the susceptibility of the electric screening mas^ 



T, 



dT 



(135) 



^^The overall scale and many details of the curve shown in figure [511 are sensitive to lattice artifacts 
394l 408{ . However, only its temperature dependence is relevant, but this has to be studied further. 
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I Electric screening mass for su(2) | | Electric screening mass for su(3) | 




T/T, T/T, 



Figure 39: The chromoelectric screening mass (top and middle panels) and the corresponding sus cept i bility 



J5|) (bottom panel) for the gauge algebra su(2) (left panels) and for su(3) (right panels) 
The line for the su(2) case is a fit of type p36p and for the su(3) case of type (I137p . 



233, Hf- 
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for both su(2) and su(3). It is therefore possible to determine the gauge-invariant phase 
transition temperature from the propagators alonJ^. Note that the minimum of the screen- 
ing mass, though close to the phase transition temperature, cannot be the location of the 
phase transition for a second order transition, since at this point the susceptibility would be 
zero instead of infinite. T he lo cation of the phase transition is thus rather indicated by the 



peak in the susceptibility |394 



For a second order phase transition, as is the case for su(2) |374| . it is somewhat surprising 
that the screening mass remains finite at the phase transition. This is even more peculiar 
since the temperature-dependence of the electric screening mass, after subtracting a zero- 
point mass, shows a critical behavior 394| . This is best illustrated using a fit in the critical 
region close to Tc of type 



DL{0)--2{T)=mo + ai9{T,-T)il 



T 

f: 



+ ahd{T - T, 



T 



1 



(136) 



where i/<> is a critical exponent. This fit is shown also in figure l39l It is a current question 
whether this critical behavior is related to the one expected for su(2) based on the universal- 
ity class of the 3-dimensional Ising model. This would require a value of the c ritic al exponent 
of z/> = z/< ^ 0.63. The fitted values are = 0.15(1) and z/> = 0.65(1) |394 |. These do 
not show the expected equality, though the exponent in the high-temperature phase is close 
to the expected one. In the low-temperature phase, the very small pre-factor a/ makes the 
determination of the exponent rather complicated, given the statistical uncertainty. The 
constant offset uiq may then be linked to the screening mass in the minimal Landau g auge , 
and thus obscuring the critical behavior. A more detailed analysis can be found in 394 . 



Taking tall of his into account, the screening mass indeed shows a behavior consistent with 
the critical scaling of a second order phase transition. 



In the su(3) case, the transition is of first order 374j |. This leads to the expectation of a 
discontinuity of the screening rnass, which is reflected in the lattice results shown in figure 
l39l Indeed, a behavior of type |394 



Dl(0)-5 (T) = mo + e{T - T,) aVS + t 



(137) 



is describing the data acceptable close to Tc, as is visible in figure [391 Thus, the order of 
the phase transition seems to be deducible from the longitudinal gluon propagator. Thus, 
all interesting properties of the phase transition - critical temperature, order, and critical 
behavior - seem to be accessible directl y from the correlation functi ons, though systematic 



effects remain to be better understood j216l 1233, M |406|, M, |409 . 

Of course, the overall scale of the figure [391 thus of the fit fll36p . is renormalization- 



group- dependent, as is the screening mass in general. However, this is a rather generic feature 



^^One should be somewhat wary with this statement. It is known that gauge-dependent qua ntities can 
indicate a (quantum) phase transition, though there is no physical transition associated with it 319l 413 1. 
However, in the pres ent case the free energy and other physical observables confirm the presence of the 
phase transition [37J|. 
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of many order parameters 39J], including such important quantities like the Polyakov loop 
and the chiral condensate discussed in section [5.3.21 Nonetheless, this implies that it cannot 
be measured directly, even if it would be gauge-invariant. 

Another nice result is that, as in the discussion of the scaling case above, the resulting 
propagator here in the minimal Landau gauge also solves both aspects of the Linde problem 
immediately. One is the appearance of infrared divergent perturbative integrals as soon 
as a magnetic propagator is involved. This has been solved previously by a regularization 



procedure and the introduction of a phenomenological magnetic screening mass |395| . This 



screening mass is now recovered here from first principles. The second component is that 
the perturbative expansion breaks down already at order (7^, since from this order on all 
orders of perturbation theory equally contributqij. This problem is also alleviated by the 
present approach, since it provides a non-perturbative expression, which does not require 
any expansion in the coupling constant. It is worthwhile to note that using the Gribov- 
Zwanziger Lagrangian these problems are explicitly o vercome, and permit even perturbative 
calculations in finite-temperature Yang-Mills theory 4151 . |416 . 



5.2.3. Infinite-temperature limit 



At the current time there is quit e a numb er of hints from both experiments |l9l-l2ll. |417 



418l | and theoretical investigations |417H419| that the matter above the phase transition is 
strongly interacting up to temperatures as large or even larger than 2-3 times the critical 
temperature. It is thus far from the naively expected quark-gluon plasma. This expectation 
was based on the idea that the coupling evaluated at the temperature becomes small with 
increasing temperature. Of course, in a renormalizable and interacting quantum field theory, 
this simple argument is not fully adequate. However, it is imaginable that at very large 
temperatures all processes are essentially dominated by leading-order perturbative effects. 

Indeed, when measuring thermodynamic observables using lattice gauge theory at very 
high temperatures, they show a behavior which could be in agreement with a logarithmic 
approach of a Stefan-Boltzmann be havior, ari d thus thermodynamic bulk quantities appear 



as being dominated by free gluons j374j . |420| . However, there are a number of arguments 



against that this is a generic behavior. Already the infrared analysis of section [5. 2. II suggests 
that even at very high temperatures, at lea st in the soft magnetic sector, non-perturbative 



long-range interactions should remain [216 . 

Indeed, there is a very general argument supporting this line of thought 
Take the limit of infinite temperature. This is equivale nt to the static limit of the theory 



194, 392 



as the extension of the time direction shrinks to zero [392| . Thus it becomes effectively 



three-dimensional. T he corresp onding three-dimensional Lagrangian can be obtained as an 



effective field theory |166l . |393| . In this course the original Aq component of the gluon field 
becomes an additional adjoint Higgs field. This preserves the number of degrees of freedom, 
as a three-dimensional gluon field has only one transverse polarization. The constants that 



^^To resolve this problem, resummation techniques have been introduced, like hard-thermal loops (39E 
414| . Though they provide an enormous improvement in the intermediate momentu m regi me, they eventually 



fail in the infrared, disagreeing with the results for the propagators presented here [19J|. 
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appear in the three-dimensional theory, hke the Higgs mass, can be obtained by matching 
with the original theory. The most important quantity in this context is the (now dimen- 
sionful) three-dimensional gauge coupling, which to leading order is given by = g'^T with 
the four- dimensional gauge coupling g. The limit actually contains some subtleties, since 
the starting theory is a renormalizable theory, while the limiting theory is fi nite i n th e sense 



that all renormalization constants are finite, at least in Landau gaugr^^' 
resulting effective Lagrangian is then given by 



190, 2211. The 



C 



1 



h 



D 



ab 



which still requires the usual gauge- fixing techniques. As noted, the coupling constants, g^, 
TTih, and h, are not free, but tied to their finite-temperature expressions, and in particular 
functions of the four-dimensional gauge coupling and the temperature only. 

Such a theory is actually confining, which could have been expected on the basis that 
a three-dimensional Yang-Mill s th eory is qualitatively not much different frorna four- 
dimensional Yang-Mills theory 42l|. The adioint scalar could actually change this 422| . but 
for the value of the parameters relevant this is not the case 3421 . l343l |. In fact, the adjoint 
scalar is a residue of the electric sector of the theory, while the Yang-Mills sector corresponds 
to the magnetic sector. Thus, it should be expected that even the infinite-temperature limit 



is strongly interacting |194 



Motivated by this it is possible to also find non-perturbative quantities in the finite 
temperature system which do not show a b ehav ior compatible with a perturbatively domi- 
nated system, like the spatial string tension 42 3| . The system is therefore non-perturbative. 
Nonetheless, it remains that the thermodynamic bulk quantities are apparently perturbative. 
The resolution of this apparent contradiction will be discussed below in section 15.31 

Of course, such a system can be equally well solved using either lattice or functional 
methods. In case of DSEs, the truncation employed is shown in figure HOI Employing a bare 
scalar-gluon vertex can be justified at least by quenched lattice results, i. e. using only test 
scalars, but not with sea scalars. The result is shown in figure HU As a consequence, the 
adjoint Higgs decouples in the infrared for any non-zero tree-level mas^. Since the tree-level 
mass is essentially the electric screening mass, and thus proportional to the temperature, 
this is an acceptable assumption. 

Also the tadpoles have to be included in a very p articular way in this approach, in 
contrast to the more simpler Yang-Mills case 1661 . Il90| . On the one hand, the finite Higgs 
mass yields a finite shift of its mass due to a Higgs tadpole. On the other hand, in principle 
any three-dimensional Higgs- Yang-Mills theory in Landau gauge features linear divergences. 
However, this particular theory is the infinite-temperature limit of a four- dimensional Yang- 
Mills theory in Landau gauge. Since temperature cannot introduce new divergences 390l |. 



^^From the perturbative point of view it is superrenormalizable, and outside Landau gauge it is actually 
linear divergent. Both points are of no relevance he re. 



80 



This can radically change for a massless scalar 190l 19. 
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Figu re 40: The graphical representation of the DSE truncation in the infinite-temperature hmit used here 
[l66j. Curly lines are gluons, dotted lines are ghosts, and dashed lines are the adjoint scalars. Lines with a 
circle are full propagators, and open circles are full vertices, which are set bare in this truncation. 



I Scalar-gluon vertex, symmetric poinT 



Scalar-gluon vertex, orthogonal momenta 




Figure 41: The quenched scalar-gluon vertex in three dimensions for su(2) in minimal Landau gauge [6211424 1 
Note that the tensor structure of the scalar-gluon vertex is the same as for the ghost-gluon vertex, ((9T|) 
The lattice size is (3.1 fm)^ with a = 0.13 fm. The mass of the quenched scalar is 1 GeV. 



and in the current truncation no new artificial ones appear 166ll216l . l22l| . the effective three- 



dimensional theory cannot have any divergent renormalization constants. This dictates that 
the additionally introduced coupling constants have to arrange themselves so that no such 
divergences arise, which has to be taken into account in the truncation, finally yielding the 



truncation used here and described in detail in 166 



Similarly, the system can be simulated on the lattice |342L l343l | . Again, it is necessary to 



tune the three indepen dent constants in the system such that the theory is finite. This is a 
non-trivial task 342], and guidance can be taken from the perturbative equivalence 
in particular since this problem arises in the deep ultraviolet. It is then possible 



393 



to determine the propagators of all involved fields in both lattice and DSE calculations. 
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Figure 42: The results for the su(2) infinite temper atur e limit from DSEs (in the scaling case) |166l | and 
lattice calculations (in the minimal Landau gauge) 342l |. The top- left panel shows the gluon propagator 
(the lattice data are for the Yang-Mills case [240j, for comparison), the top-right panel the ghost dressing 
function, the bottom-left panel the adjoint Higgs propagator, and the bottom-right panel the running gauge 
coupling. 



The results are shown in figure |42j It turns out that the presence of the Higgs is actually a 
negligible effect for the propagators of the gluon and the ghost, and henc e th e resu lt for them 
is almost identical to the one of a three-dimensional Yang-Mills theory 1661 . l342j . This is in 



line with the argumentation above using functional methods that the chromoelectric sector 
essentially decouples at sufficiently high temperature, being dominated by its screening mass. 
The latter encodes the infiuence of the hard modes on the chromoelectric sector in the infinite 
temperature limit. 
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Figure 43: The soft transverse (left panel) and longitudinal (right panel) gluon Schwinger functions above 
the phase transition temperature. The top panels show the results from lattice calculation in the mi nima l 
Landau gauge [g^I, the middle panels using DSEs in the scaling case at about 1.5 Ge V tem perature 221 1, 
and the bottom panel in the infinite temperature limit using DSEs in the scaling case jl66 ] . Note that the 
'time' in this equilibrium case is actually along a spatial direction, as discussed in section [4.21 
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5.2.4- Schwinger functions and confinement 

Once the propagators are available, it is natural to investigate once more what their ana- 
lytic structure is. From the infrared analysis of section 15.2.11 it can already be expected that 
the transverse gluon propagator will not be an asymptotic state in very much the same way 
as in four dimensions. Indeed, the corresponding Schwinger function shows this behavior, as 
can be seen in figure H3] for the finite temperature case , though the employed lattice volumes 



are not sufficient to make any final decision yet [409| . Only that the Schwinger function is 
not monotonously decaying at these volumes can be observed. Note that for the not shown 
hard modes their effective tree-level mass is so large that it is not yet possible to find any 
significant deviation from a massive particle-like behavior, nor any trace of their thermal 
width. 

Thus, the magnetic sector again appears to be essentially unaffected by temperature. In 
particular, this implies that magnetic gluons are not asymptotic states at all temperatures, 
in contrast to the idea of a plasma of free gluons. However, given the problems encountered 
when trying to construct a gauge-invariant gluon state, this outcome is rather com forti ng. 



The situation is somewhat different for the chromoelectric gluon. First of all [216[, at 
the perturbative level the zero component of the gluon in Landau gauge is a member of an 
elementary (perturbative) BRST quartet j2|,[l2|- such, it should notbe observable, i. e., it 



is effectively confined. This is also not altered at finite temperature 39d |. However, pending 



a final understanding of a non-perturbative BRST construction, this is at the current time 
not fully satisfactory. 

Therefore, again its Schwinger function is quite interesting, which is also shown in figure 
H3| for both the infinite and finite temperature casj^. At infinite temperature, both the 
propagator and the Schw inger function can be well fitted by the expression fllOOl) for 



an unstable particle [166|. Indeed, the propagator is in very good agreement with a double 
complex pole structure. It is even possible, at least in the infinite temperature case, to 
identify a subleading pole structure. The lattice data, however, are not yet sufficiently 
precise to make a final statement, except that the behavior is different from the transverse 
case: The Schwinger function, and a possible zero crossing of it, seems to depend stronger 
on temperature than in the transverse case. Thus, it is in principle possible that there could 
be a change of the analytic structure either when going to non-zero temperature or at the 
phase transition. In particular, the latter option is interesting, as this would be closer to 
the idea of a deconfining transition. 

5.3. Thermodynamic quantities 
5.3.1. Thermodynamic potential 

With all these non-pert urbative interactions around, it is a non-trivial question, why var- 



ious results [374j . l382l . |420| show that for asymptotically large temperatures a Stefan-Boltz- 



mann-like behavior is observed for bulk thermodynamic quantities. It has been shown, em- 



^^It should be noted that the result is qualitatively independent of the interaction strength of the Higgs 
with the gluons in the infinite-temperature case [19Cll | , a nice display of how even arbitrarily weak interactions 
can have an influence. 



131 



ploying the dimensional reduction as a tool, that the thermodynamic potential can be quite 
accurately reproduced using perturbation theory and only one further effective parameter to 
incl ude sub-le ading non-perturbative contributions of the three-dimensional Yang-Mills the- 



ory 



382, 396 



These results show that the non-perturbative, effectively three-dimensional, 
contribution is sub-leading by at least one power of the temperature compared to the hard 
interactions, and thus the bulk thermodynamic quantities at sufficiently high temperatures 
are essentially given by the perturbati ye ex pression. Arguments based on DSEs |l66| and 
the Gribov-Zwanziger effective theory 4151 ] are in-line with this conclusion. 



From the physical point of view this is just the statement that the average energy scale 
of interactions is given by the temperature, and this is then deep in the asymptotic domain. 
However, this is essentially an off-shell process, since, as discussed previously, the gluons 
themselves are likely confined. Therefore, a more appropriate picture may be a dense glueball 



soup, which undergoes rapidly hard collisions dominated by hard partonic processes [2161 . 

Given the manifestation of the phase transition in the correlation functions, it is an 
interesting question whether also the other thermodynamic information can be obtained from 
the correlation functions. Indeed, it is, in principle, possible to determine the t herrnodynami c 
pote ntial by a 2 P I / Luttinger- War d- Cor nwal l- Jackiw- Tomb oulis construction 2271 . 1266L 1425 - 



43l[. The exact expression is given by |431 
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where tl denotes the tree-level propagators. The contribution V2 contains contributions 
from the vertices, which are dropped in the truncation at propagator level. In the infinite- 



temperature limit, this expression reduces to [166 
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where Zq, Gq, and Hq are the tree- level dressing functions of the dimensionally- reduced 
theory ( ITggD 

Zo{p) = Go{p) = 1 
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By dropping the expression V2, although consistent with the truncation scheme of the func- 
tional equa tions, this approximation has the serious drawback of not closing thermodynam- 
ically 43 1| . In particular, this yields spurious divergences in the necessary cut-off, which 

132 



Contributions 



I — 

0%l 



02 



-0.002 



-0.004 - 



-O.OOl 



-A--A- 



. I 



0.95 



1.05 



I Total I 



a 

0.001 



_I I I I I I L_ 



-0.001 



1.1 0.9 

T/T. 





Soft mode ghost 


- ♦ - 


Lowest hard mode ghost 




Soft mode magnetic 


- -■- - 


Lowest hard mode magnetic 




Soft mode eiectric 


- -A- - 


Lowest mode eiectric 




Aii soft modes 


- - 


Aii iowest hard modes 


■o ■ 


Fuii 



■Q 



G... 



■O. 



O. 



0.95 



1.05 



1.1 

T/T, 



Figure 44: Results for the approximation (I139P of the thermodynamic action using the soft mode and the 
lowest hard mode and a three-dimensional cut-off of 1.5 GeV for the lattice data [62.] . 



are very hard to determine at finite temperature. The s purious di vergences can be better 



isolated and subtracted in the infinite-temperature hmit |l66l . Il90| , yielding an expression 



proportional to the temperature cubed, and therefore as expected s ub-l eading compared to 



the Stefan-Boltzmann behavior generated by the hard interactions [3 731 . 

A possibility to obtain a comparison at finite temperature is to use the lattice result in 
the expression (11391) at the same truncation leveQ However, it suffers from discretization 
artifacts, which require to include an artificially low cut-off of about 1 to 2 GeV with the 
available lattice data. The corresponding temperature-dependence of this approximation is 
shown in figure HH It is essentially dominated by the truncation of the Matsubara sum, and 
has little resemblance to the expected behavior. This indicates that possibly contributions 
from higher energy and higher Matsubara frequencies are important. This is to be expected, 
since the cut-off imposed is below 27rTc 1.8 GeV, and thus below the natural scale of 
the hard interactions at the phase transition. Thus interesting in principle, this approach 
requires much better input. 

Of course, other bulk quantiti es, h ke the entropy, can in principle be determined as soon 



as this problem is under control j373 



5. 3. 2. Order parameters 

As already noted, the critical temperature and the order of the phase transition can 
be accessed rather directly from the propagators. It is therefore interesting whether it is 
possible to also determine other order parameters. Of course, this is in principle not of direct 
releva nce, since QCD likely does not exhibit a genuine phase transition at finite temperature 



at all |432| . though this may change for the standard model [43 3| . Moreover, for most gauge 



'I am grateful to Jan M. Pawlowski for motivating this investigation. 
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alge bras the finite temperature phase transition of Yang-Mills theori es is first order 375 



380(1 . Therefore no order parameters need to be associated with them [152]. However, it has 



been quite a striking observation that t he or der parameters employed for su(2) Yang-Mills 



theory, in particular the Polyako v lo op |180| . remain to be rather good indicators for the 



first order ph ase transitions [375H38Cll ]. or even for the cross-overs observed in theories with 



432l . l434j - |438j . It is thus worthwhile to determine them using the present 



matter fields 
setup. 

This will be done here for the one central order parameter in Yang-Mills theory, the 
temporal Polyakov loop [43|] 



trexp ( / dtAo 



:uv. 



where the closed curve C wraps once around the compactified time direction. It is associated 
with the center symmetry Zat of SU(A^) theor\j^. being non-zero if the symmetry is broken 



180 



Unfortunately, the expression for the Polyakov loop (11411) is an exponential in the fields. 
Therefore, as for any Wilson line ||9i], it is an infinite series of contracted correlation functions 
of arbitrary high order. Therefore, an exact evaluation using correlation functions is in 
general not possible. 

As a consequence, a number of approximate evaluation schemes have been developed, 
aside from the direct evaluation using lattice gauge theory 43|. Two of them are closely 
tied to the correlation functions. 

One i s based on the Weisz-potential to determine an upper bound of the Polyakov loop 
This approach involves only the ghost and gluon propagators. 



61, 266, 354, 377 



and 

even with the approximation of using their zero-temperature behavior provides rather good 
qualitative and even quantitative results. This indicates that the dynamics of the Polyakov- 
loop at the phase transition is not tied strongly to the finite-temperature behavior of the 
propagators, but is dominated by the interaction strength mediated by the gluon propagator 
essentially at the scale of the temperature. 

Another possibility is to use dual observables, which denote a class of observables ob- 
tained by performing Fourier transformations of quantities with respect to the (fermionic) 
boundary condition of a test particle. I. e., the boundary condition is taken to be not just 
(anti-Vperiodic, but periodic up to a phase exp The simplest example is using a test 
quarlo- The momentum integral of the trace of the quark propagator S yields then the 



^•^Strictly speaking, this symmetry has to be factored out in the standard model for anomaly reasons [8|. 
However, even for thcories which do not possess a non-trivial center, the Polyakov loop is an acceptable 
order parameter [377 -380). and can therefore be expected to remain so also for the standard model. It is 
then, howe ver, no longer related to a symmetry-breaking transition, and is also not st rictly zero in either 
phase [439| . though to a very good approximation it turns out to be so [379l[380l |439| . There are further 
arguments concerning the product-group structure of the standard model suggesting its usefulness )433j] . 

*'*An alterati on to a sea quark corresponds to the introduction of an imaginary chemical potential |440|, 
as discussed in 12661. 
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Polyakov loop for SU(2) | 



Polyakov loop for SU(3) | 




Figure 45: The Polyakov loop in the quenched approximation obtained from the d ual chira. 1 condensate 
(|142p by Fourier-transformation (|143p for su(2) (left panel) and su(3) (right panel) |237l |442| . The lower 
panels show a magnification around the critical temperature. 



gauge-invariant dual chiral condensate jssol . 44 1| 

(#), = Z2C^t5^ J 0^^^^s{p,2nT{n + ^)), (142) 

where Z2 is the wave-function renormahzation constant of the quark, and the trace is over 
both Dirac and color space. Of course, a.t ip = n this is just the usual, renormalization- 
group-dependent quark condensate. 

It can now be shown that the corresponding Fourier transform 

"^^e-^'^(^V^), (143) 
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equals the Polyakov-loop in the quenched case 44 1| . Thus, it is only necessary to determine 
the trac e of the quenched quark propagator 350|]. This can be done either by lattice calcu- 
lations 443 1 or using functional equation^. The results obtained with the latter method 
are shown in figure HS] for the quenched su(2) and su(3) case. Though the involved approx- 
imations make the Polyakov loop not exactly zero in the low temperature phase, the phase 
transition is clearly discernible. However, since the available input data, as shown in figure 
[39| do not clearly signal the properties of the phase transition yet, the distinction of su(2) 
and su(3) is also not completely evident in the Polyakov loop yet. Thus, again better data for 
the gluon propagator would be necessary for this purpose. Nonetheless, this demonstrates 
how it is possible to access thermodynamic order parameters using correlation functions. 

5.4- Vertices 

As remarked, vertices become tremendously complex at finite temperature due to the 
large nurnber of po ssible vertex structures once Lorentz/Euclidean symmetry is no longer 



manifest [390l . l395l |. As a consequence, at non-zero temperature there have been only few 



investigations of the ver tices beyond perturbative evaluations yet. The ghost-gluon vertex 



has been investigated in |223[, and found to be moderately affected by the temperature, but 
in a way intricately linked to the propagators. 

In the infinite-temperature limit it is possible to also make some statements about the 
vertices, since the tensor structure becomes simple once more. However, even in this case, re- 
liable investigations without quenching the adjoint H iggs are n ot yet available. The quenched 
case has been shown in figures |27l [28], andHD |62ll89l . ll53l . l424| for the ghost-gluon, the three- 
gluon, and the two-adjoint-Higgs-gluon vertices, respectively. The former two are identical to 
the ones of pure Yang-Mills theory, and show a very similar behavior to the zero-temperature 
case. For the latter it is found that the vertex is essentially fiat and close to tree-level, sup- 
porting the approximations made in the DSE treatment of the infinite-temperature limit in 
section I5.2.3[ though this requires confirmation in the unquenched case. 



6. Summary and discussion 

6.1. The state of the art, pitfalls, and odds and ends 

Summarizing, quite a lot of progress has been made in the determination of (gauge- 
dependent) correlation functions since the first investigations using functional 177l | or lat- 



tice methods 95 



At the present time the technical tools for investigations using lattice 
methods in Landau gauge are quite well under control. In functional methods, a wealth of 
understanding of different approximations and truncations has been obtained. As a con- 
sequence, sophisticated truncation schemes have been developed, which on the one hand 



permit to reach rather good agreement with lattice results for the propagators [37|, while 



at t he same time already give access to more complex quantities like three-point vertices 



89l Il53l . 1154 . Il95l | , and to some extent four-point vertices |l55l . II6OL Il82| , as well as genuine 



85i 



Where currently actually lattice result s for the gluon propagator are used as input, and an assumption 



has to be made on the quark-gluon vertex [2371 l350l l371| . 
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bound-states [13|, l333[ | and thermodynamic properties |237l . l341l . l354| . In particular the 
latter two points bring the functional approach, though quite less intensively developed than 
lattice methods, to the same general level: Being a tool for investigating physical processes. 

This progress has been made possible by a synergistic development of lattice and contin- 
uum methods, which also includes non-functional approaches, like effective theories. Only 
this combination has made it possible to have as good a control over and understanding of 
the correlation functions as has been presently achieved. 

With this a turning-point has been reached. The investigation of correlation functions has 
so far mainly been conducted for its own sake, and for the exploratory investigation of how to 
access physical observable^ and to understand more about fundamental non-perturbative 
phenomena like confinemen t, chiral symmetry breaking and topology directly on the level 



of correlation functions [13|, l289j . Now the technology developed is sufficiently advanced 
to start accessing completely new fields. And there are two avenues of investigations, which 
are currently followed. One is the application of the combined framework of lattice and 
functional meth o ds to new sys tems. First successes have been achieved in this program 



M \M m uM I342l - l344j . 1362 



[444] . but a lot of work remains to bring to bear the full 
possibilities of the synergistic effects. The other is to actually use the functional methods 
to extrapolate the lattice results into domai ns where the l attice cannot venture currently 



for p ractical reasons, like cold, dens e matter |227l l34ll . l349l . 13571 . l359l . |444| . non-equilibrium 



445| . and scattering processes |275| . It can therefore be expected that the combination of 
methods will provide access to a number of unresolved problems in the future. 

This said, there are still topics to be understood, even on a fundamental level. Though 
quite a number of hints have been presented here for the idea of how the various solutions 
available fit together, this is far from a proof. Being skeptical, it could easily turn out that 
indeed all possible ways to select Gribov copies inside the first (or any fixed) Gribov region 
could coincide in the thermodynamic limit. In particular, the discussion given in section 
14.3.21 implies that it is very well possible that the scaling case of functional calculations can 
never be connected with a feasible numerical simulation, though this does not imply that it 
does not exist. 

However, as has been repeatedly emphasized, this is actually not a hmitation , since th e 



physical results turn out to be independent of the choice of solution |333l l339l . |370| - |372 



including questions like confinement [37|, |354| . Thus, in principle, it is possible to sit back 



relaxed, and wait what is realized, and thus available, in Yang-Mills theory, and all the while 
keep on working with the already confirmed set of correlation functions. Their properties 
are anyhow equally well under control. 

One has to be a bit less relaxed when thinking about their interpretation. Though by 



now a non-perturbative realization of BRST symmetry |57l . l59l . l60l . Il30| . and with it an 



algebraic construction of the state space |l2l . l37l |. is available, it is by far not clear how 



*^0f cour se, in particular concerning the hadronic spectrum these have aheady been quite successful 
HQim, but so far with somewhat ad-hoc approximations, mainly guided only by chiral symmetry and 



neglecting most of the Yang-Mills sector. The reasons why these drastic approximations can deliv er a rathe r 
good description of physical phenomena can only with the full solutions now be fully understood 217 , 333l | . 



137 



to connect it to any of the proposed gauges on the lattice, as emphasized above. On the 
other hand, in the finite- ghost cases, it is still not entirely clear how a Hilbert space has 
to be constructed 
localities 



71, 298 



371 1X1 
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and in particular this may involve explicit non- 
Though this is not something with direct relevance to 
the determination of observables, it is a pain in the back of the head, and conceptually not 
yet satisfactory. In particular, though manifest quark and possibly gluon confinement has 
been obtained from correlation functions 37|, IGJ, l214l 12771 . |284J . l354| . its full understanding 
remains a challenge. 

Another unclear question is yet how the string-like properties of Yang-Mills theory are 
realized. This entails two questions. One is, whether it is in prin c iple poss i ble t o encode the 



string-like features in low-order correlation functions a t all |180L Il95l . 1 1961 . |289| . Since they 



are already manifest in, e. g., the Regge trajectories |180| of mesons encoding the string 
tension, the answer is an affir mative yes, since the meson properties can be determined from 



the four-point functions |273| . But this implies that possibly full self-consistency at the level 
of four-point functions may be required, and this is definitely not yet achieved. It is also 
not clear whether and how a corresponding truncation can be determined. The situation is 
even worse for string-breaking, which requires at least six-point functions to observe. 

Even more complicated is the Wilso n potenti al. Determining it may require rather 



been obtained also in Landau gauge |195 



different gauges, like the Coulomb ga uge 14461 l447l| . though some approximate results have 



355| . But even with the approximate results in 



place, the understanding of how to realize the difference between Yang-Mills theory with 
different gauge groups or for Y ang-Mills theory coupled to matter in different representations 



remains a subject of research [1281 . Il94j . Il97l . Il98l . l448 



Finally, though again exploratory results are available, the connection of the correlation 
function approach to other views of low-energy Yang-Mills theory and QCD still requires a 
better unders tand i ng. I n particu l ar, only how t o relate correlation functions to topological 



scenarios 



44J. I255I . |285| . I324J . l328H330l . l449H45l| as well as to stochastic quantization ideas 



65L I94J . I122I . Il29| has been investigated so far. Still a direct relation to many other views, like 



the stochastic vacuum, analytic perturbation theory, and many others remains unknown. 

6.2. Summary and outlook 

In total, in this manuscript an overview has been given on the determination of (gauge- 
dependent) correlation functions with both lattice and functional methods at zero and finite 
temperature in gauge theories. Particular attention has been paid to the problem intro- 
duced by the Gribov-Singer ambiguity. The emphasis has been put on Yang-Mills theory, 
with a simple Lie group, though connections to other theories and/or setups have been indi- 
cated. The results here have concentrated on gauges fulfilling th e perturbative Landau gauge 



condition. Th erefore, rnany resul ts in Coulomb gauge [17 



13, 61 
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lating gauges [I7ll. l342 
Abelian gau g es |l51 



483 



linear covariant gauges 



4931- 



0495 



47g-'483, 

485-49^1, direct Laplacian gauges [53l 
2191 . 1220 . 49617 (49 7| have not been covered, though results at a similar 



446 .1447.1452447811 . interpo- 



d. I13.4 14831 . 1484 - maximal 



, and other gauges 



level have been obtained in these cases as well. In particular, a similar situation with the 
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finite-ghost case and scaling appears to be realized in Coulomb gauge |45J, l460l . l464l |473| , 
though this has to be investigated further. 

Of course, the current important topic is the full inclusion of matter fields to investigate 
QCD, the standard model, and beyond-the-standard-model theories. In the end, however, 
all this will rest on the concepts presented here, and will be achieved with similar methods. 
It is thus the creation of this technological framework and the understanding of its results 
within the last fifteen years, which is now paving the way for a complete new generation 
of applications. Thus, with the next decade it will be decided whether the promise of the 
possibilities these methods and approaches hold will be realizable, and perhaps even open 
up the way to new frontiers in non-perturbative physics. 
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